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ABSTRACT 

Gamma Ray Burst (GRB) prompt emission spectra in the keV-MeV energy range are usually 
considered as adequately fitted with the empirical Band function. Recent observations with the 
Fermi Gamma-Ray Space Telescope {Fermi) revealed deviations from the Band function, sometimes 
in the form of an additional black-body (BB) component, while on other occasions in the form of an 
additional power law (PL) component extending to high energies. In this article we investigate the 
possibility that the three components may be present simultaneously in the prompt emission spectra 
of two very bright GRBs (080916C and 090926A) observed with Fermi, and how the three components 
may affect the overall shape of the spectra. While the two GRBs are very different when fitted with 
a single Band function, they look like “twins” in the three-component scenario. Through fine-time 
spectroscopy down to the 100 ms time scale, we follow the evolution of the various components. We 
succeed in reducing the number of free parameters in the three-component model, which results in a 
new semi-empirical model -but with physical motivations—to be competitive with the Band function 
in terms of number of degrees of freedom. From this analysis using multiple components, the Band 
function is globally the most intense component, although the additional PL can overpower the others 
in sharp time structures. The Band function and the BB component are the most intense at early 
times and globally fade across the burst duration. The additional PL is the most intense component at 
late time and may be correlated with the extended high-energy emission observed thousands of seconds 
after the burst with Fermi /Large Area Telescope (LAT). Unexpectedly, this analysis also shows that 
the additional PL may be present from the very beginning of the burst, where it may even overpower 
the other components at low energy. We investigate the effect of the three components on the new 
time-resolved luminosity-hardness relation in both the observer and rest frames and show that a strong 
correlation exists between the flux of the non-thermal Band function and its E pea k only when the three 
components are fitted simultaneously to the data (i.e., F^ T -Epg^ ki relation). In addition, this result 
points toward a universal relation between those two quantities when transposed to the central engine 
rest frame for all GRBs (i.e., L^ T -Ep® s */| T relation). We discuss theoretical implications of the three 
spectral components within this new empirical model. The BB component can be interpreted as the 
photosphere emission of a magnetized relativistic outflow. The Band component can be interpreted as 
synchrotron radiation in an optically thin region above the photosphere, either from internal shocks 
or magnetic field dissipation. The extra power law component extending to high energies likely has 
an inverse Compton origin of some sort, even though its extension to a much lower energy remains a 
mystery. 

Subject headings: Gamma-ray burst: individual: GRB 080916C- Gamma-ray burst: individual: 

GRB 090926A Radiation mechanisms: thermal Radiation mechanisms: non- 
thermal - Acceleration of particles 
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1. INTRODUCTION 

The fireball model remains the most popular scen ario 


& Rees) 1978 

Paczynski 1986; Goodman 

1986; Sherm & 

Piran 

1990[ 

Rees & Meszaros||1992( Meszaros & Rees 

1993 

iees & Meszaros| 1994). In this m 

odel, the GRB 

le or a rapidly 
star formed by 
star (collapsar; 

centra 

spinni 

either 

engine is a stellar-mass black ho 
ng and highly-magnetized neutron 
the collapse of a supermassive 

Woosleyl 

1993; IMacFadyen & Woosleyl 19991 IWoosley 

& Heger 

2006) or the merger of twe 

compact ob- 


In both cases, the original explosion creates a bipolar 
collimated jet composed mainly of photons, electrons, 
positrons and a small fraction of baryons. The relativis¬ 
tic explosion ejecta within the jet are not homogeneous— 
they form multiple high density layers, which propagate 
at various velocities. When the fastest layers catch up 
with the slowest, the charged particles contained in the 
layers are accelerated through mildly relativistic colli¬ 
sionless shocks (inter nal shocks; Rees & Mcszaros]|1994 


Kobayashi et al. 1997 Daigne & iVlochkovitch|I998|). The 
particles subsequently cool via emission processes such as 
synchrotron, Synchrotron Self Compton (SSC), and In¬ 
verse Compton (IC). The internal shock phase is usually 
associated to the so-called CRB prompt emission, mainly 
observed in the ke V—MeV energy range (see e.g., the 
spectral catalogs by |von Kienlin et al7||2011[ |Gruber et| 
al.|2014 ) and usually lasting from a few ms up to several 
tens to hundreds of seconds. As the ejecta interact with 
the interstellar medium they slow down via re lativistic 


collis i onless shocks (externa l shocks; Rees & Meszaros 
1992 Meszaros fe Rees|1993l accelerating charged parti¬ 


cles, which then emit non-thermal synchrotron photons. 
This external shock phase is usually associated to the so- 
called CRB afterglow emission observed at radio wave¬ 
lengths up to X-rays and in some cases even up to the 
GeV regime hours after the prompt phase, and days and 
even years for the lowest frequencies. The detailed origin 
of the gamma-ray emission, however, is not fully under¬ 
stood and many theoretical difficulties remain, such as 
the composition of the jet, the energy dissipation mech¬ 


anism s. as well as the radiation mechanisms (e.g., Zhang 
20111 . 


"Another prediction of the fireball model is the exis¬ 
tence of intense, thermal-like emission from the jet photo¬ 
sphere, expected to be observe d simultaneously with the 
non-thermal prompt emission (Goodman|1986| Meszaros 
2002 |Daigne fe Mochkovitch |2002l Rees fe Meszaros 
20051. Indeed, the high density of the outflow makes it 
optically thick to Thomson scattering close to the source, 
but as the jet expands its density decreases and radia¬ 
tion can escape resulting in a thermal-like component 
more or less affected by sub-photospheric processes and 
jet-curvature effects. 

GRB prompt emission spectra in the keV-MeV energy 
regime we re predominantly f i tted by the so-called Band 
function (Band et al.|[l993 Greiner et al.|fl995) prior 
to the launch of the Fermi Gamma-ray Space Telescope 
(hereafter Fermi ) in 2008. This empirical function is a 
smoothly broken power law whose shape is defined with 
four parameters: two indices a and j3 corresponding to 
the spectral slopes of the low- and high-energy power 


laws (PL), respectively; a break energy parametrized 
to cor respond to the maximum of the spectrum 
Ep ea k (|Gehrels| 1997]); and a normalization factor. The 


derived Band spectra usually suggested a non-thermal 
origin of the radiation leading to the natural proposition 
that they were produced by synchrotron mechanisms. 
However, the high values obtained for a were often in 
conflict with the synchrotron scenario predi ctions in both 


the s l ow and fast electron cooling regimes ( Crid er et al., 
1997 Preece et al. 1998 Ghisellini et al. 2(J(Jl)| ) With' 


the aim to id entity emission from the fireb all model’s jet 


photosphere, Ghirlanda et al. (2003) and Ryde (20041 
fitted a thermal spectral shape—using a pure black- 
body component (BB), and a combination of a BB and 
a PL, respectively—to the GRB prompt emission ob¬ 
served with the Burst And Transient Source Experiment 
(BATSE) on board the Compton Gamma Ray Observa¬ 
tory ( CGRO ). Despite the good fits obtained in a few 
cases, it was often difficult to assess whether thermal 
spectra were be tter than the non-therm al Band ones. At 


the same time, Gonzalez et al. (2003) fitted simultane¬ 
ously BATSE and 6 'G"^G/Energetic Gamma Ray Exper¬ 
iment Telescope (EGRET) data of GRB 941017, finding 
a significant deviation from the Band function at high en¬ 
ergies (> 1 MeV); these fits improved with the addition 
of a PL to account for the high-energy deviation. 

The Gamma-ray Burst Monitor (GBM) and the Large 
Area Telescope (LAT) on board Fermi have significantly 
improved GRB prompt emission observations after 2008. 


Although the Band function remains a good model to 

describe GBM spectra (e.g., Gruber et al. 

2014), the 

joint spectral analysis of some GRBs using 
and LAT confirmed that in some cases an ad 

Doth GBM 
ditional PL 

is required to improve the fit quality (e.g., 

Abdo et al. 

2009bj Guiriec et al. 112010; Ackermann et al. 

■iOioMMil 

2013p. Moreover, Fermi results show that this additional 


tion b elow a few tens of keV. Indeed, jAckermann et al~ 


(20111 reported the existence of an intense peak in the 


prompt emission light curves (LCs) of GRB 090926A, 
from the lowest to the highest energies, associated to this 
additional PL. Similar associations of sharp structures 
in GRB LCs with an a dditi onal PL were al s o rep orted 


Guiriec et al. 


(2011b) and Gonzalez et al. (2012), the 


latter using GURU data. In GRB U90926A, the addi¬ 
tional PL exhibited a spectral break around 1.4 GeV, 
which was interpreted as resulting from 7-7 opacity mak¬ 
ing possible an estim ate o f the bulk Lorentz factor, T, 
between 200 and 700. Guiriec et al. (2010) reported that 
a similar PL deviation from tEe" Band function could 
also be identified in GBM-only data of some GRBs, sug¬ 
gesting a spectral turnover well below the 1.4 GeV of 
GRB 090926A, since little or no emission was observed 
in the LAT data for these events. The origin of this 
spectral component remains challenging. 

The quest for GRB prompt photospheric emission con¬ 
tinued in the Fermi Era. Within the framework of 
the fireball model, the expected photosphere compo- 
nent should outshine the non-thermal Band component 
(|Zhang & Pe’er||2009|, but this photosphere component 
remained undetected. Such a surprising observational 
result triggered a heated debate on the GRB prompt 
emission mechanism. A highly-magnetized jet was sug¬ 
gested to suppress the photosphere emission component 
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(|Daigiie_fc_^lochkovitch 2002 Nakar et al. 2005 Zhang 


Pe’er 2009|). Alternatively, it was suggested that the 


sphere (e.g., 

Beloborodov 2010 Lazzati et al. 2011 

Rees 

& Meszaros 

2005 

Thompson 2006). Ryde et al. ( 

2010 ) 

and P^r et 

al. 

2 U12) fitted simultaneously a thermal 


and suggested that the main keV-MeV prompt emission 
could be of thermal origin and, the refor e , the signat ure 


of the jet’s photosphere. Then, Guiriec et al. (2011a) re¬ 
ported for the first time that the prompt emission vF v 
spectra of GRB 100724B were best fitted with a dou¬ 
ble curvature model (BB+Band) rather than the Band 
functi o n alo ne. Contrary to the co nclusions of |Ryde| 


et al. (2010) and Pe’er et al. 


(2012 ), the photosp heric 


component (BB) identified in Guiriec et al. (|2011a|) was 
subdominant compared to the non-thermal component 
(Band function). Such a subdominant BB is not ex¬ 
pected from the pure fireball model, but is more con- 



nent, the outflow must be highly-magnetized close to the 
source but the jet must have a low magnetization at large 
radii for t he in ternal s hocks to be efficient. Alternatively, 


Zhang & Yan (2011b) suggested that efficient energy dis¬ 
sipation is still possible, even if the magnetization pa¬ 
rameter a is still greater than unity. They envisaged a 
collision-induced magnetic reconnection and turbulence 
(ICMART) to effecti vely dissipate the magn etic energy in 
the outflow (see also Zhang & Zhang||2014[ for a simula¬ 


tion of GRB lightcurves within such a model.) A similar 


GRBs ((Guiriec et al. 2011a|b 

2013 

Burgess et al.J|2011 

Axeisson et ai.|20I2) observed with i 

^'ermi. Interestingly 


firmed in Guiriec et al. (2013), fits using a Band function 


with a BB result fn Band function shapes that are much 
more compatible with the predictions of the sy nchrotro n 
emissi on origin. Moreover, as shown by |Guiriec et al.| 


(2013), in those GRBs in which an intense but still sub- 
dominant thermal-like emission is detected, a strong cor¬ 
relation appears between the energy fluxes and the 
spectral peak energy of the non-thermal component only 
(hereafter Ff T -E^ k)i relatiorj^J) For GRBs in which 
the thermal-like contribution affects very little the shape 


the data leads to a similar F-E pea k relation (Golenetskii 

et al. 1983 

Borgonovo & Ryde||2001; Liang et al. 

200 TF 

Guiriec et a 

1. 2010[ Ghirlanda et al 2010, 2011a|b[ 

Lu et 

al. 

2012 

). The Fj N1 -Fp e ; ak;i relations have similar slopes 


the non-thermal prompt emission, and, when corrected 
for the redshift, the fits between the luminosity and Ep® s * k 

of the non-thermal component (hereafter Lf JT -E(^ k I 'f T ) 
align perfectly for all tested GRBs (Guiriec et al. 2013). 


In summary, until now GRB prompt emission spectra 


16 NT stands for non-thermal component; in the context of the 
multi-component model that includes thermal and non-thermal 
components, this non-thermal component is usually adequately fit¬ 
ted with a Band function or a cutoff power law. 


have been found to be composed of: 1) a Band function 
alone; 2) a combination of a Band function and a PL; 3) 


a combination of a BB and a PL; 

or 4) a combination of a 

Band function and a BB (Abdo et al.|2009a|b Guiriec et 

al.,2010 

Ackermann et al.|2010, 2011; Guiriec et al. 2010, 

20 lTaTb 

2013; Pe’er et al. 2012 

Rydeetai. 2010; Zhang! 

et al. 201 lcf. The non-therma 

Band function usually 


component usually carries less than a few tens of percent 
of the total energy, and the additional PL component 
extends from low to high energies. 

Here we reanalyze two of the brightest Fermi GRBs, 
GRB 080916C and GRB 090926A, in the context of the 
multiple components described above and discuss the 
possible simultaneous existence of such separate com¬ 
ponents by comparing spectral analyses at various time 
scales from time-integrated to very fine-time intervals. 
Through the evolution of these components with time, 
we trace back the emission mechanisms which produced 
them, following the jet physical mechanisms back to the 
central engine energy reservoir. In Section [2j we de¬ 
scribe the data selection as well as the observations, and 
in Section [3] we present the procedures we followed for 
the analysis. Sections [4] and [5] are dedicated to the re¬ 
sults of the time-integrated and the coarse-time analy¬ 
ses, respectively. Section [6] reports the fine-time spectral 
analysis results, which are the main topic of this arti¬ 
cle. In Section [TJ we investigate the impact of the mul¬ 
tiple spectral components on the relations between the 
energy-flux of the non-thermal component and its ob¬ 
served vFv spectral-peak energy (namely Fb T -Ep e J ik i ) 
and between the luminosity of the non-thermal compo¬ 
nent and its intrinsic vF v spectral-peak energy (namely 

L[ fT -E™ s * k ^ T .) Section|8]introduces our suggested model 
for GRB prompt emission. Section |9] discusses an inter¬ 
pretation of our observational results. Finally, in Sec¬ 
tion lOj we summarize the most important observational 
results and comment on their future impact. 

2. DATA SELECTION AND OBSERVATIONS 

The GBM consists of 12 Sodium Iodide (Nal) detec¬ 
tors covering energies between 8 keV and 1 MeV. GBM 
also includes two bismuth germanate (BGO) detectors 
(bO, bl) covering energies between 200 keV and 40 MeV, 
located on opposite sides of the spacecraft. The LAT 
is a pair-conversion telescope sensitive to photons with 
energies from 20 MeV to >300 GeV. When a photon en¬ 
ters the LAT, it is converted into an electron-positron 
pair when interacting into the conversion tungsten foils 
located between the tracker silicon strip detector planes. 
Detailed information about the GBM and the LAT can 


be found in Meegan et al. (2009) and Atwood et al. 


(2009), respectivelj 17 
GBM detected GKB 080916C and GRB 090926A on 


2008 September 8 a t T o =00:12:45 UT (Goldstein & 
van der Horst| 2008 ) and on 20 09 September 26 at 
l'o=04:20:26.99 UT ( Bissaldi|[2009 ), respecti vely.. These 
bursts were also detected with the LAT (Tajima et 

~ (12 009 a j) for GRB 080916C, 


al. 


Abdo et al. 


and Uchara et al. 


Ackermann et al. (2011) for 


17 See also http://fermi.gsfc.nasa.gov for additional and up-to- 
date information about t'ermi. 
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Time since GBM trigger (s) 


Fig. 1.— GRB 080916C : Light curves as observed with GBM Nal (8-500 keV) and BGO (500 keV-20 MeV) detectors and with 
LAT/LLE (>20 MeV). (al—6) 0.1 s time-resolved light curves, (bl—6) Light curves with the binning used for the fine-time analysis 
presented in Section [6] 


TABLE 1 

List of GRB properties. 


GRB Namt 

a 

Location 

Redshift z 

t 9 o(s b 

Observed Highest Energy Photon (GeV) 

GRB 080916C 
GRB 090926A 

(RA,Dec) = (119.84717°,-56.63833°) (±0.5” 
(RA,Dec)=(353.40070°,-66.32390°) (±1.5’” 

C 

r 

4.24±0.2(P] 

2.106ft- 1 

63±1 

20±2 

33 

19.6 


a See GBM GRB catalog at http://fermi.gsfc.n asa.gov/ssc _ 

“Duration computed between 50 and 300 keV (Ivouveliotou et al. 
The Gamma-Ray bu rst Opti cal j Xear-Ii 1 frarect I Sel ect or (CiltU.N 


iGreiner et al. 

(2009 


c Switt/'U VUT 

Vetere et al. 2009 


f Very Large Telescope (VLT) (JVI 

alesani et al. 2009) 


GRB 090926A.) The properties of the two GRBs are 
reported in Table ITJ We show the light curves of 
GRB 080916C ana GRB 090926A in several energy 
ranges in Figures [T] and [2j respectively. 

For the analysis we only selected the optimal Nal de¬ 
tectors with angles to the source smaller than 30° based 
on the best source locations and with no blockage by an¬ 
other part of the spacecraft (e.g., LAT, radiators, solar 
panels) as well as with no shadowing by another GBM 
module. According to these criteria, we selected detec¬ 
tors nO, n3 and n4 for GRB 080916C and n3, n6 and n7 
for GRB O9O926A0 


1 19931_ 

) - Clemens et al.|2~008|l 


Usually, the source is only in the field of view of one 
BGO detector, and only this BGO detector is retained 
for the analyses. Therefore, for GRB 080916C, we used 
bO. However, in the case of GRB 090926A, the burst 
occurred in the median plane of the spacecraft, which 
is parallel to the sides to which the BGO detectors are 
attached. Therefore, the two BGO detectors can be si¬ 
multaneously used in this specific case. Because of the 
peculiar location of GRB 090926A in the spacecraft coor¬ 
dinates, the reconstructed signal may be distorted due to 
the inaccurate modeling of the absorption in the back of 
the BGO modules (photomultipliers for instance), espe- 


18 The Nals are named nx with x varying from 0 to 9 for the first 


10 detectors, and “a” and “b” for detectors 11 and 12, respectively 
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Fig. 2.— GRB 090926A : Light curves as observed with GBM Nal (8 keV—1 MeV) and BGO (1-20 MeV) detectors and with LAT/LLE 
(>20 MeV). (al—6) 0.1 s time-resolved light curves, (bl-6) Light curves with the binning used for the Hne-time analysis presented in 
Section |6] 


daily at low energy. Our analysis revealed that although 
the signals reconstructed from the various Nal detectors 
match perfectly, discrepancies between the Nal and BGO 
modules are observed in the overlapping energy region 
between 200 and —750 keV for bO and between 200 and 
-600 keV for bl for GRB 090926A. Beyond -750 keV 
and —600 keV for bO and bl, respectively, Nal and BGO 
data are in agreement. 

We used in our analysis time tagged event (TTE) data, 
which have the finest time (i.e., 2 /is) and energy resolu¬ 
tion that can be achieved with the GBM. We used the 
Nal data from 8 keV to —900 keV cutting out the over¬ 
flow high-energy channels as well as the K-edge from —30 
to —40 keV. For GRB 080916C, we used BGO data from 
200 keV to —39 MeV cutting out the overflow energy 
channels. For GRB 090926A, we also excluded the low- 
energy channels <750 keV and <600 keV for bO and bl, 
respectively, since they may be affected by calibration 
problems. We then generated the response files for each 
GBM detector based on the best source location. 

To optimize the analysis according to the instrument 
performance, we study the photons converted in the front 
and in the back of the LAT tracker as two separate data 
sets (hereafter FRONT and BACK, respectively) using 
the appropriate instrument response functions. We an¬ 
alyzed the FRO NT and BACK data pas sing the Pass 
7 transient cuts (Ackermann et al. 2012b) within a 10° 
region of interest (KOI) centered on the best source loca¬ 


tion and with energy above 100 MeV. To perform spec¬ 
tral analysis of the LAT data between 20 and 100 MeV, 
we used the LAT low-energy (LLE) data, which are de¬ 
signed to increase the LAT sensitivity below 100 MeV 
for transient sources as well as to make possible spectral 
anal ysis belo w 100 MeV (for instance, see Appendix of 
Ajelio et al.||2014 |. 


1’he background in each of the GBM detectors as 
well as in the LLE data was estimated by fitting poly¬ 
nomial functions to the light curves in various energy 
ranges before and after the source active time period. 
The background was then measured by interpolating the 
functions to the source active time period. For GBM 
data, the background was fitted to the CSPEC data 
files, which have the same energy resolution as the TTE 
data but with a coarser time resolution. CSPEC data 
cover a much longer time range, making the estimation of 
the background more reliable especially for long GRBs. 
We then exported the background estimated from the 
CSPEC data to our analysis of the TTE data. The back¬ 
ground for the LAT FRONT and BACK data was esti¬ 
mated by averaging the data over several orbits (during 
which the GRB was not present) when Fermi was lo¬ 
cated at the same position in the orbit as at the trigger 
time and when the LAT w as pointin g in the same direc¬ 
tion as at the trigger time (Vasilelou] 2013). 
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Energy (keV) Energy (keV) Energy (keV) Energy (keV) 

Fig. 3.— Sketches of the various tested models. 


3. SPECTRAL ANALYSIS METHODOLOGY 

To perform the spectral analysis, we used the XSpec 
and Rmfit tool kits. Starting from a cutoff power law (C) 
or a Band function (Band), we built twelve increasingly 
complex spectral models with two or three components 
adding a power law (PL), a black body (BB), C or Band. 
Sketches of the various models are presented in Figure [3] 
The best spectral parameter values and their l-a un¬ 
certainties were estimated by optimizing the Castor C- 
statistic (hereafter Cstat), which is a likelihood technique 
that converges to y 2 for a specific data set when there 
are enough counts. We performed the spectral analysis 
over three time scales—we refer to these later in the text 
as: time-integrated, coarse-tinre bins and fine-time bins. 

In order to check the consistency of the Nal and BGO 
detectors through all tested models, we fitted both detec¬ 
tor types simultaneously and we applied a free effective 
area correction (EAC) factor between bO and each of the 
Nal detectors selected for the analysis. For each model 
fitted we then fixed these EAC values. To evaluate the 
consistency of the GBM data (Nal and bO) with the LAT- 
LLE, LAT-FRONT and LAT-BACK data, we also fitted 
them simultaneously adding a new, free EAC factor be¬ 
tween bO and each of the LAT data types. When EAC 
values were compatible with no correction needed, they 
were set to unity. We only applied the EAC corrections 
in the time-integrated and coarse-time bin analysis, be¬ 
cause in the fine-time bin analysis the number of counts 
is much smaller, leading to unconstrained EAC factors. 
When significant EAC factors (i.e., significant deviation 
from unity) between the GBM and LAT data sets are 
required with the simplest models, those data sets are 


in much better agreement with the more complex ones. 
The strong corrections required with the simplest models 
seem to be unreasonable based on our current knowledge 
of the instruments; therefore, this reinforces the scenarios 
involving the more complex spectral shapes. The results 
of cross-calibration analysis are reported in greater detail 
in Appendix [A] 

Due to the very limited number of counts at high en¬ 
ergies, we only used the LAT-LLE data in combination 
with GBM in the coarse-time analysis, and we did not 
use any LAT data in the fine-time analysis. We then 
estimated the probability that the Cstat improvements 
obtained with the most complex models compared to the 
simplest ones are not merely due to signal and/or back¬ 
ground statistical fluctuations, by performing multiple 
sets of Monte Carlo (MC) s imulations following th e same 
procedure as discussed in |Guiriec et al. |2011aJ 2013| 
(see Appendix |B|) In Sections pT] and [5] we compare the 
results of joint GBM and LAT data fits with those of the 
GBM-only ones for the time-integrated and the coarse¬ 
time analysis, respectively; the consistency of the joint 
GBM and LAT analysis with the GBM-only one is im¬ 
portant to support the results of the fine-time analysis 
presented in Section | 6 j for which we only use GBM data. 
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Fig. 4.— Time-integrated GBM and LAT spectra of GRB 080916C (from Tq- 0.1 s to To+71 s) when fitted (a) with a Band function 
alone and (b) with B+BB+C. The count rate spectra are presented in panels (al) and (bl)—the solid lines corresponding to the fitted 
model—and the deconvolved z/F„ spectra—the dashed lines correspond to the individual components of the fitted model and the solid ones 
to the total emission—in panels (a2) and (b2). Panels (a3) and (b3) correspond to the residuals of the fits. 



Fig. 5. — Time-integrated GBM and LAT spectra of GRB 090926A (from To to To+20 s) when fitted (a) with a Band function alone, 
(b) with B+C and (c) with B+BB+C. The count rate spectra are presented in panels (al), (bl) and (cl)—the solid lines corresponding 
to the fitted model—and the deconvolved i/p„ spectra—the dashed lines correspond to the individual components of the fitted model and 
the solid ones to the total emission—in panels (a2), (b2) and (c2). Panels (a3), (b3) and (c3) correspond to the residuals of the fits. 


4. TIME-INTEGRATED SPECTRAL ANALYSIS 

We performed a time-integrated spectral analysis over 
a period corresponding to the most intense part of the 
prompt emission of the two GRBs in the keV-MeV en¬ 
ergy range (i.e., from T o -0.10 s to To+71.00 s and from 
T 0 to To+20.00 s for GRBs 080916C and 090926A, re¬ 
spectively.) The fit results using the v ario us mod els dis¬ 
cussed above are reported in Tables |A1| & A2 of Ap¬ 
pendix [C] We discuss below the fits and present the most 
relevant ones in Figures [4]&[5j 


4.1. GBM-only 

For GRB 080916C, we find that the B+BB and B+PL 
fits improve the Band-only fit by 74 and 54 units of Cstat 
for 2 additional degrees of freedom (dof), respectively, 
while for GRB 090926A the improvement is 139 and 81 
units of Cstat, respectively 


tivelyP 1 ] 
d in [Hu 


As already reported in IGuiriec et all (12011 al 120131), 


19 In the case of nested models and in the gaussian regime an 
improvement of ~10 units of Cstat per additional degree of freedom 
corresponds to a ~5 a level improvement. 
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the addition of the BB component to the Band func¬ 
tion results in shifting systematically the E pea k values 
toward higher energies and making both the low- and 
high-energy Band PL indices steeper. The values of 
a obtained with the B+BB fits are, therefore, more 
compatible with the pure slow-cooling synchrotron emis¬ 
sion scenario for the two GRBs and even with the pure 
fast-cooling synchrotron process for GRB 080916C. The 
impact of the parameter changes on the F^-Ep^ ; 

and L^ T -Ep® s ^^ T relations as well as on the interpre¬ 
tation of the prompt emission mechanisms will be dis¬ 
cussed in great detail in Sections [7] and [9} The temper¬ 
atures of the BB are ~40 keV for GRB 080916C and 
GRB 090926A corresponding to the Planck function 
peaking at ~120 keV (right within the sensitivity range 
of the Nal detectors.) 

It must be noted that when the high-energy PL of the 
Band function has an index <-3, the Band function can 
be replaced with C without affectin g m uch the Cstat 
value of the fit (see for instance Table A1 of Appe ndix [C] 
or time interval from To-0.1 s to To+4.3 s in Table [A3| of 
Appendix [C] ) This does not mean that the high-energy 
slope is as steep as an exponential cutoff, but that the 
slope is at least as steep as ~-3 and that our data may 
not allow us to better measure /3. In the time-resolved 
analysis and in the context of the multi-component mod¬ 
els, we will not make any distinction between Band and 
C for the main keV-MeV spectral contribution. 

Adding a PL to the Band function (i.e., B+PL) has 
th e opposite effect on the Band parameters, as reported 


Guiriec et al. (2010): E pea k is shifted toward lower en¬ 
ergies and both low- and high-energy Band PLs become 
less steep. The index of the additional PL is ~-2.00 for 
both GRBs. 

Although we performed simulations to compare all 
models, here we only describe the results obtained when 
comparing Band-only fits with C+BB+PL because they 
are the most relevant ones for the rest of the analy¬ 
sis. Indeed, Band-only is traditionally considered as 
a good fit to the data and can, therefore, be consid¬ 
ered as a reference —for instance, it ha s been the case 
for GRB 080916C (fAbdo et al.l 2009a])- but we show 
throughout this article that C+BB+PL is a globally bet¬ 
ter description of the data for the two GRBs by following 
the statistical test procedure presented in Appendix |B| 
For GRB 080916C, C+BB+PL improves the Band-only 
fit by 77 units of Cstat for 3 additional dof, while for 
GRB 090926A, C+BB+PL leads to an improvement of 
162 units of Cstat for 3 additional dof. For both GRBs, 
none of the 10 5 synthetic spectra using the Band-only 
fits as null hypothesis give a ACstat value as high as 
the observed ones, thus supporting the statement that 
the probability that C+BB+PL better fits the data than 
Band is due to statistical fluctuation of signal and back¬ 
ground around the Band function is <10 -5 . Moreover, 
the Cstat value resulting from the fit of the real data 
with Band-only is much higher than the Cstat values ob¬ 
tained when fitting the synthetic spectra, which would 
not be expected if the true spectrum were a Band func¬ 
tion. We performed the same exercise, this time using 
the fit of C+BB+PL to the real data as the null hypoth¬ 
esis and recovered Cstat values for the synthetic spectra 
that were compatible with the ones obtained when fit¬ 


ting the real data with either Band or C+BB+PL. All 
the above results reinforce the hypothesis that the Band 
function is not the best description of the data and that 
C+BB+PL is a significantly better model. 

4.2. GBM+LAT 

Panels (a3) of Figures E] & [5] show strong wavy pat¬ 
terns in the residuals of the fits when a Band function 
alone is fitted to the GBM+LAT time-integrated spectra 
of GRBs 080916C and 090926A.The spectral parameter 
values resulting from the simultaneous fit of GBM and 
LAT data with the various models—and especially the 
most complex ones—are very well compatibl e wi th tho se 
reported from the GBM-only fits (see Tables Al & A2 of 
Appendix [Cj ). While the high-energy power law of the 
Band function is statistically compatible with an expo¬ 
nential cutoff when fitting B+BB+PL and C+BB+PL 
to the time-integrated GBM data alone, B+BB+PL is 
preferred over C+BB+PL in the simultaneous fit of the 
time-integrated GBM and LAT data. 

We investigated the effect of a possible cutoff in the 
additional PL of the B+BB+PL model by replacing the 
PL with a PL with an exponential cutoff (i.e., B+BB+C) 
or with a second Band function (i.e., B+BB+B). We 
find that for GRB 080916C, B+BB+C does not re¬ 
duce the Cstat value obtained with B+BB+PL. For 
GRB 090926A, both B+BB+C and B+BB+B are sig¬ 
nificantly better than B+BB+PL. 

5. COARSE-TIME SPECTROSCOPY 

We concluded above that both GRB 080916C and 
GRB 090926A were better fitted with a three-component 
model. However, it is impossible to conclude from the 
time-integrated spectra whether these separate compo¬ 
nents are artifacts due to a strong spectral evolution 
during the event duration, and whether they can be as¬ 
sociated directly to physical processes. To address these 
questions, we performed a coarse-time spectral analysis 
selecting time intervals including the main stru ctur es o b- 
served in each burst light curve (see Table A3 and A4 of 
Appendix |C] for GRB 080916C and GRB U90926A, re¬ 
spectively. (This selection resulted in 3 and 5 time inter¬ 
vals for GRB 080916C and GRB 090926A, respectively. 

We fitted either the GBM data alone or in combination 
with the LLE data. We did not include the LAT transient 
data in this analysis to avoid complications and possible 
biases due to the low number of photons and possible 
calibration inconsistencies. We also tested the need of 
EAC either among the GBM detectors or between the 
GBM and the LLE data; the results are reported in Ap¬ 
pendix [A] 


5.1. GRB 080916(7 
5.1.1. GBM-only 

We find that C+BB is a significantly better fit than 
Band only in the first and third time intervals, improv¬ 
ing the Cstat values by 28 and 45 units, respecti vely , for 
only one additional dof. As reported in Section [4+| the 
addition of a BB component to the Band function sys¬ 
tematically shifts Ep ea k to higher energy and decreases 
both values of a and /3. In fact, the slope of the high- 
energy PL of the Band function becomes so steep that it 
can be replaced with C. That is why C+BB and B+BB 
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have similar Cstat values in these two time intervals. In 
the second time interval, the addition of the BB compo¬ 
nent to a Band function does not improve the fit signifi¬ 
cantly (i.e., 11 units of Cstat for two additional dof), and 
it impacts much less the Band parameters. In particular, 
the value of /3 remains pretty high and inconsistent with 
an exponential cutoff. In all intervals, the temperature 
of the BB is similar to the temperature measured in the 
burst time-integrated spectrum. 

Adding a PL component to the Band function does not 
work for the first time interval (i.e., no convergence is 
obtained with a PL compatible with a positive flux). Al¬ 
though no improvement of the Cstat value is observed in 
the second time interval either, it is possible to fit B+PL 
to this time interval, and the fit results in PL with a 
positive flux (i.e., the flux of the additional PL at 100 
keV is 3.7 er above 0.) In the latter case, the fitting en¬ 
gine prefers a solution for which most of the high-energy 
emission is associated with the PL and not with the Band 
function. To enable this solution, the high-energy slope 
of the Band function becomes very steep and compatible 
with an exponential cutoff, becoming identical to slopes 
of the two other intervals when adding a BB to the Band 
function (i.e., B+BB is equivalent to C+BB in term of 
Cstat values.) Therefore, C+PL is equivalent to B+PL 
(i.e., similar Cstat values) in the second time interval 
with one less dof. In the third time interval, a B+PL 
fit significantly improves the Band only fit (i.e., 39 units 
of Cstat for 2 additional dof); however, C+BB remains 
the best two-component model. In the second and third 
time intervals, the indices of the additional PL are com¬ 
patible with the values resulting from the time-integrated 
analysis. 

Since for the three time intervals, C+BB (or B+BB) is 
significantly better than B+PL (or C+PL), we investi¬ 
gated the effect of an additional PL to the C+BB model 
(i.e., C+BB+PL). In the first time interval, we could 
not obtain an additional PL with a positive flux. In the 
second time interval, adding a PL to the B+BB model re¬ 
duces the Cstat value by 13 more units and impacts the 
parameters of the Band function. With the additional 
PL, Epeak is shifted toward lower energy and while the 
value of the low-energy PL index of the Band function 
increases, its high-energy PL becomes compatible with 
an exponential cutoff. The changes in both slopes mark 
the contribution of the additional PL not only at high- 
energy but also at low-energy. C+BB+PL is equivalent 
to B+BB+PL (i.e., similar Cstat values) and an improve¬ 
ment of 13 units of Cstat is obtained compared to B+BB 
for one additional dof. In the third time interval, the 
Cstat values obtained with C+BB and C+BB+PL are 
very similar; however, the addition of the PL affects the 
spectral parameters of the other components in a way 
that the PL has a positive contribution to the total flux 
(i.e., the flux of the additional PL at 100 keV is ~13 er 
above 0 for the 2 nd and the 3 rd time intervals, respec¬ 
tively.) Finally, we investigate the possibility of a break 
in the additional PL of the C+PL scenario in the third 
time interval by replacing the PL with a second cutoff PL 
(i.e., C+C2). With this additional dof, the fitting engine 
prefers to have C2 mimicking a thermal shape instead of 
the non-thermal shape of the additional PL. Indeed, the 
index of the C2 is +0.95+0.22 and Eq is 64+7 keV which 


matches very well with the BB component measured in 
the B+BB and C+BB scenarios. 

5.1.2. GBM+LLE 

The spectral results obtained when fitting simultane¬ 
ously GBM and LLE data remain globally similar to 
those derived from the fits of GBM data alone, both for 
the parameters and for the statistics. However, the LLE 
data allow for a better constraint of the high-energy PL 
of the Band function in the most complex models such as 
B+BB, B+PL and B+BB+PL, and confirm steep slopes 
for the high-energy PL of the Band function with indices 
below -2.5. In addition, possible breaks have been iden¬ 
tified in the additional PL when using LLE data with 
Eo between ~100 and ~200 MeV in the second time in¬ 
terval, and between ~250 and ^700 MeV in the third 
time interval—because Eo=E pea k/(2-apL) with E pea k be¬ 
ing the energy at the break in the photon PL spectrum 
and q;pl being the PL index, the actual break energy 
E P eak in the photon PL spectrum is between ~50 and 
~100 MeV in the second time interval and between ^75 
and ~210 MeV in the third one. 

5.2. GRB 090926A 

5.2.1. GBM-only 

A B+BB fit improves significantly compared to the 
Band-alone fit in all the time intervals but the last, in 
which the intensity of the signal is weaker. A B+PL 
leads to similar improvements although the B+BB re¬ 
mains slightly better. Before T 0 +5s, the Cstat value 
decreases by 67 units for two dof when adding a BB to 
the Band function; a B+PL model, however, cannot be 
constrained. An additional PL is identified only after 
T 0 +5s in the coarse-time analysis. In that case, both 
E p e a k and P values remain similar to those measured 
in the Band-only scenario; however, the value of a in¬ 
creases. Both the temperature of the BB and the index 
of the additional PL are similar to those measured in the 
time-integrated analysis. This reinforces the idea that 
the three components identified in the time-integrated 
spectrum are not artifacts due to strong spectral evolu¬ 
tion as they are recovered in the time-resolved analysis 
and from one time-interval to the other. 

From 5s to 9s after the trigger time, B+BB+PL im¬ 
proves B+BB by 13 units of Cstat. The resulting high- 
energy PL of the Band function is very steep and compat¬ 
ible with an exponential cutoff. Therefore, C+BB+PL 
is as good as B+BB+PL and an improvement of 9 units 
of Cstat compared to B+BB is measured for only one 
additional dof. After To+9s, a limited improvement is 
obtained when adding a PL to B+BB; however, the fit 
results indicate that in the C+BB+PL model, the fluxes 
of the three components are positive (for instance the 
fluxes of the additional PL at 100 keV are 9.0, 4.4 and 
1.2 a above 0 for the 3 rd , the 4 th and the 5 th time inter¬ 
vals, respectively.) The reason the Cstat values do not 
change much is that the global shape of the whole spec¬ 
trum changes, especially the low- and high-energy PLs 
of the Band function. 

5.2.2. GBM+LLE 

The simultaneous fits of GBM and LLE data give sim¬ 
ilar spectral results to those obtained when fitting GBM 
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data alone with the various tested models; however, LLE 
data enhance the significance of the additional BB and 
PL components. Indeed, before T 0 +5s, B+BB is better 
than Band alone by 97 units of Cstat while an improve¬ 
ment of 67 units of Cstat was obtained from GBM only 
data. In the second and third time intervals, B+PL is 
clearly a better improvement than B+BB, while the two 
models were competitive when using only GBM data. 
This is expected since only the additional PL contributes 
to the emission detected at high-energy with the LAT. Fi¬ 
nally, we included an exponential cutoff in the additional 
PL to look for a high-energy break, but the significance 
of this result remains weak due to the low data counts in 
these energies. 

5.3. Partial Conclusion 

The coarse-time spectral analysis supports the complex 
spectral shape identified in the time-integrated spectra of 
GRBs 080916C and 090926A as well as the possible pres¬ 
ence of multiple spectral components. The various com¬ 
ponents identified in the time-integrated spectra of these 
two GRBs cannot be explained as artifacts due to strong 
spectral evolution during the event durations. Indeed, we 
recovered the multiple components in the smaller time 
intervals with values of their spectral parameters consis¬ 
tent with those obtained in the time-integrated analy¬ 
sis. Furthermore, this refined analysis suggests that the 
three components do not have the same relative inten¬ 
sities during the burst durations. The BB component 
and the Band function seem to be the most intense at 
early times, while the contribution of the additional PL 
is stronger at later times. 

6. FINE-TIME SPECTROSCOPY 

To better understand the behavior of the various spec¬ 
tral components as well as their evolution with time, we 
performed a very fine-time spectral analysis. At very 
fine time scales, GBM and LAT data are in very dif¬ 
ferent count regimes; therefore, to avoid any complica¬ 
tion related to this difference as well as to possible cross¬ 
calibration issue we only consider GBM data throughout 
this Section. 

With fine time intervals it is often difficult to unam¬ 
biguously conclude that the most complex models are 
significantly better than the simpler ones based on the 
pure statistical considerations of the likelihood ratio test 
solely. Instead, our goal here is to verify that the com¬ 
plex shapes identified as significantly better than the 
Band function in the time-integrated and the coarse-time 
spectral analyses are also valid options at very fine time- 
scales. It is important to keep in mind that even if the 
most complex models are not statistically better than 
the Band function based on the likelihood ratio test re¬ 
sults, it does not necessarily mean that the contributions 
of the additional components (i.e., BB and PL) to the 
total flux are negligible and compatible with a null flux. 
As presented in the previous Sections, the addition of a 
BB and a PL to the Band function significantly modi¬ 
fies the parameters of the Band function. Therefore, it is 
possible to have complex models that are statistically not 
favored compared to a Band function alone—based solely 
on the likelihood ratio test- but whose additional com¬ 
ponent fluxes are clearly positive. Moreover, we show 
in Appendix [5] that the results of the likelihood ratio 



Fig. 6.— Snapshot of a i/F v spectrum resulting from the fine-time 
analysis of GRB 090926A. This shows the new three-component 
model (i.e., (1) a non-thermal Band function with a high-energy 
power law index <-3 and which is statistically equivalent to a cut¬ 
off power law, C, (2) a thermal-like component adequately approx¬ 
imated with a black body component, BB, and (3) a non-thermal 
additional cutoff power law, C) overplotted with the Band-only fit 
for GRB 090926A. The solid yellow and black lines correspond to 
the best Band-only and C+BB+C fits, respectively. The dashed 
yellow and black lines correspond to the 1-cr confidence regions 
of the Band-only and C+BB+C fits, respectively. The solid blue, 
red and green lines correspond to the cutoff power law, to the BB 
component and to the additional cutoff power-law resulting from 
the best fit with the C+BB+C model (i.e., solid black line) to the 
data, respectively. The C+BB+C model exhibits a break in the 
second C component below 100 MeV. 


test are not always sufficient to decide whether the sim¬ 
pler model is truly a better description of the data than a 
more complex one and that additional statistical tests are 
sometime necessary; following the procedure described in 
Appendix [B] we conclude that C+BB+PL is globally a 
better description of the data than all the other simpler 
models. In this Section, we will focus mostly on the con¬ 
sistency and continuity of the spectral results from time 
interval to time interval and on how the complex models 
modify the values of the spectral parameters of the Band 
function. 

We have shown in panels (al-6) of Figure 111 and [2] 
the two burst light curves in various energy bands witn 
a 0.1 s time resolution; we show in panels (bl-6) the 
same light curves binned at the intervals we used for the 
fine-time spectral analysis (for the bin size determina¬ 
tion see Appendix 0; For this analysis we started by 
fitting to every time interval both Band-only and the 
model that we identified in the previous sections as the 
best choice, namely C+BB+PL. When all three com¬ 
ponents could not be fitted simultaneously, we kept as 
many of them as possible (i.e., if C+BB+PL could be fit¬ 
ted but only C+BB, we report the results from C+BB.) 
An example of the three-component fit using C+BB+C 
for one time interval of GRB 090926A is presented in 
Figure 6l The co mpl ete fit results are reported in Ap¬ 
pendix C]in Table [Ab| for GRB 080916C, and in Table ptb| 
for GRB 090926A(See also Appendix [Ej Figures [7] and 
[8] show the reconstructed photon light curves in multiple 
energy bands using the Band-only and the C+BB+PL 
analyses of GRBs 080916C and 090926A, respectively. 





A New Model for GRB Prompt Emission Spectra 


11 



Fig. 7. — GRB 080916C : (al-4) Count light curves as observed with GBM between 8 keV and 1 MeV, and with LAT/LLE between 20 
and 100 MeV. The time bins correspond to those used for the fine-time analysis; (bl-4) Reconstructed photon light curves resulting from 
the fine-time analysis with a Band function alone. The energy bands and time intervals are the same as for the count light curves (i.e., 
panels (a)) for comparison. The 20-100 MeV light curve is an extrapolation of the model (i.e., Band) fitted to the GBM energy range 
into the LAT/LLE energy domain, (cl-4) Reconstructed photon light curves resulting from the fine-time analysis with the C+BB+PL 
model. The blue, green and red lines correspond to the cutoff power law, to the black-body and to the additional power law, respectively. 
The black dashed line corresponds to the sum of the 3 components. The energy bands and time intervals are the same as for the count 
light curves (i.e., panels (a)) for comparison. The 20-100 MeV light curve is an extrapolation of the model (i.e., Band) fitted to the GBM 
energy range into the LAT/LLE energy domain; (dl-4) Reconstructed photon light curves resulting from the fine-time analysis with the 
C-bBB+PL 5 p arams model. The blue, green and red lines correspond to the cutoff power law, to the black-body and to the additional power 
law, respectively. The black dashed line corresponds to the sum of the 3 components. The energy bands and time intervals are the same 
as for the count light curves (i.e., panels (a)) for comparison. The 20-100 MeV light curve is an extrapolation of the model (i.e., Band) 
fitted to the GBM energy range into the LAT/LLE energy domain. For clarity, the uncertainties on the count light curves and on the 
reconstructed photon light curves are not shown. Therefore, those figures are for qualitative purpose. 


These figures show the time evolution of the individ¬ 
ual spectral components as well as the time evolution 
of the total emission, to be compared to the observed 
count light curves in the same energy bands. The evo¬ 
lution of the spectral parameters of the Band-on ly a nd 
C+BB+PL models are presented in Figure [9] and [TO] for 


GRBs 080916C and 090926A, respectively, and the pa¬ 
rameter distributions of the two fits are presented in Fig¬ 
ure HU to HU 

Note that in the following we sometimes use the ter¬ 
minology Band when we talk about the C component 
of the C+BB+PL scenario; indeed, the C component 
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Fig. 8. — GRB 090926A : (al-4) Count light curves as observed with GBM between 8 keV and 1 MeV, and with LAT/LLE between 20 
and 100 MeV. The time bins correspond to those used for the fine-time analysis; (bl-4) Reconstructed photon light curves resulting from 
the fine-time analysis with a Band function alone. The energy bands and time intervals are the same as for the count light curves (i.e., 
panels (a)) for comparison. The 20-100 MeV light curve is an extrapolation of the model (i.e., Band) fitted to the GBM energy range 
into the LAT/LLE energy domain, (cl-4) Reconstructed photon light curves resulting from the fine-time analysis with the C+BB+PL 
model. The blue, green and red lines correspond to the cutoff power law, to the black-body and to the additional power law, respectively. 
The black dashed line corresponds to the sum of the 3 components. The energy bands and time intervals are the same as for the count 
light curves (i.e., panels (a)) for comparison. The 20-100 MeV light curve is an extrapolation of the model (i.e., Band) fitted to the GBM 
energy range into the LAT/LLE energy domain; (dl-4) Reconstructed photon light curves resulting from the fine-time analysis with the 
C+BB+PL 5 params model. The blue, green and red lines correspond to the cutoff power law, to the black-body and to the additional power 
law, respectively. The black dashed line corresponds to the sum of the 3 components. The energy bands and time intervals are the same 
as for the count light curves (i.e., panels (a)) for comparison. The 20-100 MeV light curve is an extrapolation of the model (i.e., Band) 
fitted to the GBM energy range into the LAT/LLE energy domain. For clarity, the uncertainties on the count light curves and on the 
reconstructed photon light curves are not shown. Therefore, those figures are for qualitative purpose. 


of C+BB+PL corresponds globally to the traditional 
Band function and it is sometimes easier to refer to it 
as Band. The C and Band component are usually statis¬ 
tically equivalent as long as Band j3 is <-3. 
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Time since GBM trigger (s) 


Fig. 9.— GRB 080916C : Evolution with time of the spectral pa¬ 
rameters resulting from the fits of a Band function alone (blue) and 
of a C+BB+PL model (red) to the fine time intervals. The dashed 
regions >-2/3 and >-3/2 correspond to the domains in which the 
values of a are incompatible with pure synchrotron emission from 
electrons in both the slow and fast cooling regimes, and with syn¬ 
chrotron emission from electrons in the fast cooling regime only, 
respectively. 



Time since GBM trigger (s) 

Fig. 10. — GRB 090926A : Evolution with time of the spectral pa¬ 
rameters resulting from the fits of a Band function alone (blue) and 
of a C+BB+PL model (red) to the fine time intervals. The dashed 
regions >-2/3 and >-3/2 correspond to the domains in which the 
values of a are incompatible with pure synchrotron emission from 
electrons in both the slow and fast cooling regimes, and with syn¬ 
chrotron emission from electrons in the fast cooling regime only, 
respectively. 


6 .1. Evolution of the Various Spectral Components in 
the C+BB+PL Scenario 

The Additional PL Component 

At early times we were able to fit only C+BB for 
both GRBs, i.e., until T 0 +1.2 s and T o +4.0 s for GRBs 
080916C and 090926A, respectively, while keeping all pa¬ 
rameters of the C+BB+PL model free; no additional PL 
was required. This is consistent with the results reported 
in Section[5] Interestingly, while here we only fitted GBM 
data (i.e., <40 MeV), the additional PLs appear contem¬ 
poraneously with the rise of the LAT-LLE emission above 
20 MeV for both GRBs as shown in Figure [7] and [8] This 
is consistent with the fact that the mechanisms respon¬ 
sible for the additional PL extends from low energies in 
GBM (i.e., <1 MeV) to more than several tens of MeV 
in the LAT. 

After it kicks off, the additional PL remains present 
until the end of each burst. At first, its intensity rises 
until it overpowers the other components and reaches 
its maximum at s and ^10 s for GRBs 080916C 
and 090926A, respectively. The peak intensity of the 
PL is perfectly correlated in time with a very intense 
and sharp structure observed in both light curves (see 
Panels (a) of Figures [7] and [+ For GRB 090926A, the 
PL is so intense at its maximum that it overpowers the 
two other components at all energies: this results in a 


sharp spike lasting ^0.1 s visible at all energies in the 
count light curves. Since the relative intensity of the 
PL to the other components—between 100 keV and 1 
MeV—is lower at the peak for GRB 080916C, the sharp 
structure—also lasting for ^0.1 s—is only clearly visible 
<20 keV and >20 MeV in the count light curves. In 
GRB 080916C the PL remains after the peak the most 
intense component at late time especially below 20 keV. 
In the case of GRB 090926A, the PL seems to fade faster 
than the C component in the C+BB+PL model, but it 
remains until late times. 

While the intensities of the additional PLs evolve with 
time, their indices remain between — L0 and —2.0 ac ros s 
the burst durations (see Figure [9fe &[l0fe and Figure Fl2jl 
&[l4ji). For GRB 080916C the values of the PL index 
are ~-1.5 during the whole burst (see Figure [9ji fc[T+.) 
For GRB 090926A, the PL index values decline steadily 
and also cluster around -1.5 (see Figure pd}' &jl4ji). 

We note that the index of the additional PL during 
the prompt phase is similar with the PL indices reported 
for the extended LAT emission seen up to thousands of 
seconds after the prompt phase and contemporaneou s 
with the afterglow emission (e.g., Ackermann et al. 2014). 
This may suggest that the latter is the result of the same 
physical process that produced the PL in the prompt 


emission. 
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Fig. 11.— GRB 080916C : Distribution of the spectral param¬ 
eters resulting from the fine-time analysis using a Band function 
alone. The filled blue histograms correspond to the distribution of 
the best parameters (i.e., central values) and the dashed red and 
green histograms corresponds to the distributions of the 1 —<7 lower 
and upper limits, respectively. The vertical black dashed line at- 
2/3 and at -3/2 correspond to the limits above which the values of 
a are incompatible with pure synchrotron emission from electrons 
in both the slow and fast cooling regimes, and with synchrotron 
emission from electrons in the fast cooling regime only, respectively. 


Fig. 13. — GRB 090926A : Distribution of the spectral param¬ 
eters resulting from the fine-time analysis using a Band function 
alone. The filled blue histograms correspond to the distribution of 
the best parameters (i.e., central values) and the dashed red and 
green histograms corresponds to the distributions of the 1-cr lower 
and upper limits, respectively. The vertical black dashed line at 
-2/3 and at -3/2 correspond to the limits above which the values of 
a are incompatible with pure synchrotron emission from electrons 
in both the slow and fast cooling regimes, and with synchrotron 
emission from electrons in the fast cooling regime only, respectively. 
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Fig. 12.— GRB 080916C : Distribution of the spectral param¬ 
eters resulting from the fine-time analysis using a the C+BB+PL 
model. The filled blue histograms correspond to the distribution 
of the best parameters (i.e., central values) and the dashed red and 
green histograms corresponds to the distributions of the 1-cr lower 
and upper limits, respectively. The vertical black dashed line at 
-2/3 and at -3/2 correspond to the limits above which the values of 
a are incompatible with pure synchrotron emission from electrons 
in both the slow and fast cooling regimes, and with synchrotron 
emission from electrons in the fast cooling regime only, respectively. 


Fig. 14. — GRB 090926A : Distribution of the spectral param¬ 
eters resulting from the fine-time analysis using a the C+BB+PL 
model. The filled blue histograms correspond to the distribution 
of the best parameters (i.e., central values) and the dashed red and 
green histograms corresponds to the distributions of the 1-cr lower 
and upper limits, respectively. The vertical black dashed line at 
-2/3 and at -3/2 correspond to the limits above which the values of 
a are incompatible with pure synchrotron emission from electrons 
in both the slow and fast cooling regimes, and with synchrotron 
emission from electrons in the fast cooling regime only, respectively. 
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The BB Component 

For both GRBs, the BB is more intense at early time 
(i.e., during the first 4 to 5 s of each burst) than at late 
time (see red curves in Figure a & [8j;). Although its 
intensity decreases with time, it is possible to identify 
this component during the entire prompt phase. For 
both GRBs, the BB temperature peaks around 30-40 
keV with a larger spread of values for GRB 090926A (see 
Figure \l2\: fc [lT^) . While a global cooling trend may be 
observed tor GRB 090926A, the temperature seems to re¬ 
main constant for GRB 080916C (see Figure |9ji &[Ioji). 
Although the BB is a subdominant component over the 
entire energy range, it approaches and even sometimes 
overpowers the C component at the peak en erg y of the 
Planck spectrum around 100 keV (see Figure [7p & §)■ 

The C Component 

The C component carries most of the prompt emission 
energy and it clearly overpowers the other components 
between 100 keV and 1 MeV. After a fast rise at early 
times, the intensity of the C component decays steadily 
with time (see Figure [7b &[8b). This simple shape struc¬ 
ture is much less evident when a Band function alone 
is fitted to the data. For both GRBs, the a index of 
C, remains more or less constant with time, especially 
for GRB 080916C (see Figure J 9 |d &JI0J3). Although the 
values are more spread for GRB 090926A, both distribu¬ 
tions are centered around —0.7 (see Figure 12 r & 14r). 
In addition, the l-er lower limits of a are between 1.5 
and ~—0.7 —beside a few outliers—, which are the limits 
for synchrotron fast and slow cooling, respectively. 

It is very important to note that the values of a are 
much lower (i.e., <-l) in some time intervals where only 
C+BB can be fitted to the data (see first time interval of 
GRB 080916C as well as the two first and last time inter¬ 
vals of GRB 090926A in Figure [ 9)3 &[10]d, respectively.) 
In these intervals, the values of a are closer to the mean 
value of the additional PL index than in the other time 
intervals. As we will see in Section [8j this most likely 
indicates that what we think is the contribution of the 
C component at early times might actually be due to 
the additional PL contribution, which is very intense at 
early and late times. Therefore, the additional PL may 
be the first visible component in some bursts, and this 
component may last much longer and even extend into 
the afterglow phase. 

Finally, the E pea k globally decreases with time from 
> several MeV to hundreds and even tens of keV at late 
times (see Figure [9^, &fl0^i). The E pea k also tracks the 
Band function energy flux as we will see in Section [7] 

6.2. Band vs C+BB+PL 

If the prompt emission is indeed composed of sepa¬ 
rate spectral components as proposed in the C+BB+PL 
scenario, the fit of any single component (i.e., a Band 
function) to the data should correspond to an average of 
the complex model shape. It is indeed what we o bser ve 
for GRB 080916C and GRB 090926A (see Figure jAlj &; 
A2 of Appendix |E| : the yellow lines corresponding to 


the C+BB+PL scenario, the systematic wavy pattern 
observed in the fit residuals when a Band function alone 
is fitted to the data (see panel (a3) of Figure 0&§ is 
naturally expected: taking into account the various com¬ 
ponents not only reduces the Cstat value, but it globally 
flattens the systematic wavy pattern of the fit residuals. 

Generally, the fit of a Band function alone to a spec¬ 
trum composed of a Band function and a subdominant 
BB component results in spectral parameters that are 
systematically biased: both a and ft are greater than 
their ‘real’ values, and the value of E pea k is underesti¬ 
mated. On the other hand, a Band-only fit on a spectrum 
which is truly Band+PL overestimates fi and underesti¬ 
mates E pea k and a. This is in perfect agreement with 
the results presented in Sec tions [ 4j an d [5l as well as with 


those already published in Guiriec et al. (2010 2011a 


the Band function fits are averages of the black ones, 
which correspond to the sum of the three components 
of the C+BB+PL model. Therefore, in the context of 


2013). 

We discuss below two cases, one where the additional 
PL is intense compared to C but the BB contribution is 
limited (e. g., ti me interval 26.0-32.0 s of GRB 080916C— 
see Figure Al of Appendix |E|), and one where the BB 
component is intense relative to C but the additional PL 
is limited ( e.g., time interval 1.5-1.8 s of GRB 090926A 
see Figure |A2| of Appendix |E|) 

(i) time interval “26.0-32T) s” of GRB 080916C: The 
Epeak of the Band-only fit is located between the two 
humps of the C+BB+PL model; the low-energy hump 
of the latter corresponds to the energy at the maximum 
of the BB component and the high-energy one to the 
Epeak of C. At the same time, the high-energy slope 
of the Band-only fit (i.e., ft) is much harder than the 
high-energy slope of C (i.e., exponential cutoff) of the 
C+BB+PL fit, because the high-energy power law of the 
Band-only fit is an average between the high-energy emis¬ 
sion of C and of the additional PL component. Finally, 
and more importantly, the value of the Band-only a in¬ 
dex is mostly an average of C and PL, since these are the 
most intense components in this time interval. Indeed, 
a of the Band-only fit is ~ —1.1, which is about the 
average of the C and the PL component slopes, whose 
indices are ~ —0.7 and ~ —1.5, respectively. The re¬ 
sults derived here can be extrapolated to all intervals for 
GRB 080916C. Indeed, if we closely examine the evolu¬ 
tion of a over time for GRB 080916C (see Figure]^)) and 
compare it to the intensity evolution of the various spec¬ 
tral components in the C+BB+PL scenario (see panels 
(c) of Figure[7]), we can easily explain the evolution of the 
parameters of the Band function fits. Figure [9b shows 
that the Band a values (blue line) exhibit a dramatic 
change around 4-5 s: a drops from ~ —0.7 before the dis¬ 
continuity down to < —1.0 after it. This discontinuity in 
the values of a matches perfectly with the time at which 
the additional PL starts to overpower the other compo¬ 
nents, especially at low energies that have the strongest 
impact on the Band a (see panels (c) of Figure [7]). The 
values of Band-only a are the lowest (i.e., <-1.3) around 
5-6 s which corresponds to the intensity peak of the ad¬ 
ditional PL in the C+BB+PL model. However, the a of 
the C+BB+PL fit remains constant around —0.7 during 
the entire burst, except for the first time interval where 
it is < —1.1 (red line). We will see in Section [8] that this 
low value of a at early times in the C+BB+PL model 
may be an indication that the additional PL is already 
present at very early times in some GRBs. Finally we 
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note the strong changes in the values of /?, which vary 
from ~ —1.7 to < —5.0 (see Figure |9^). The values of /3 
are the highest from ~5 to ~10 s, which is when the addi¬ 
tional PL is the most intense in the C+BB+PL scenario 
(see panels (c) of Figure[7j) Also, values of /J above —2.0 
are clearly not physical as these imply an infinite amount 
of energy, so if the Band function were indeed the correct 
fit, then a cutoff in its high-energy PL would be required 
to adequately explain the data when /3>-2. Therefore, it 
is perfectly natural to account for such strong variations 
of the Band /? index with a model including an additional 
PL at high energies as explained earlier in this Section. 

(ii) time interval “1.5-1.8 s” of GRB 090926A: The 
Epeak of the Band-only fit is also located in between the 
two humps of the C+BB+PL model. The high value of 
/3 (i.e., > —2) in the Band-only scenario is too hard to be 
physical and a spectral cutoff at high-energy should be 
present. On the other hand, such a high f) value is nat¬ 
urally accommodated in the C+BB+PL model, where f3 
would be an average of the decaying parts of BB and C. 
The value of Band a would also be an average of the BB 
and C components. Since a pure Planck function is ad¬ 
equately approximated with a C with an index of +1.0, 
and since the C component of the C+BB+PL model has 
an index of ~ —1.1, we expect to recover a value of a 
greater than ~ —1.1 when fitting Band-alone to the data 
if the true model were a combination of a C and BB. 
Indeed, when fitting Band-alone to the data, we measure 
a value of a > —0.6. We will see in Section [8] that the 
low value of a (i.e., ~-l.l) resulting from the C+BB+PL 
fit may indicate the presence of an underlying additional 
PL component at early times in GRB 090926A. With the 
same considerations as for GRB 080916C in the previous 
paragraph, we can extend the results obtained in this 
time interval to the entire burst by comparing the evo¬ 
lution of a between the Band-only and the C+BB+PL 
models (see Figure [lOjj) with the intensity evolution of 
the various components with time (see panels (c) of Fig¬ 
ure pi) . The BB component is the most intense during 
the first 4-5 s of the burst, which also corresponds to the 
time intervals where the highest values of a are measured 
when fitting Band-alone to the data. The BB component 
biases Band-only a towards higher values in this time pe¬ 
riod. From ~ 5 to ~ 10 s, C clearly overpowers the other 
components, which explains the values of a around —0.7 
obtained during this time. Finally, after 10 s, the ad¬ 
ditional PL becomes intense and biases the values of a 
down to values as low as ~ —1.5. 

While a values for GRB 080916C in the Band- 
only scenario are mostly between —1.5 and —0.7 
namely the fast and slow synchrotron cooling regimes, 
respectively—most of the m a re greater than —0.7 for 
GRB 090926A (see Figure [TIJd 

In the C+BB+PL scenario, the distributions of 
a values peak around —0.7 for both GRBs. For 
GRB 080916C the a distribution is very sharp, but it 
is wider for GRB 090926A (see Figure 12 d & 14:>); how¬ 
ever, the 1 -a lower limits are nearly all below —0.7. Con¬ 
versely, the values of Epeak are spread over hundreds of 
keV with C+BB+PL, while they are all clustered around 
300-400 keV with Band-alone fits. 

An interesting result appears when we extrapolate the 
reconstructed photon light curves obtained by fitting ei¬ 
ther Band or C+BB+PL to the GBM data (i.e., between 


8 keV and 40 MeV) into the LLE energy band (i.e., 20 to 
100 MeV). While both Band and C+BB+PL well repro¬ 
duce the count light curves in the energy bands ranging 
from 8 keV up to 1 MeV, only C+BB+PL adequately 
mimics the LLE light curve between 20 and 100 MeV 
(see panels (c) of Figure [7]& [8]). This is particularly true 
for GRB 080916C. Indeea, while C+BB+PL adequately 
reproduces the observed position of the peak intensity of 
the 20-100 MeV light curve as well as its relative inten¬ 
sity after 8 s, the Band-only fits result in a peak position 
which is shifted to earlier times and in an extended high- 
energy emission which is too intense before 20 s and too 
weak afterwards. 

Finally, we performed MC simulations to investigate 
the validity of the C+BB+PL model based on the 
method described in Appendix[B] For each time interval, 
we generated 10 5 synthetic spectra choosing the Band- 
only fit as the null hypothesis; these were then fit either 
with Band or C+BB+PL. When fitting a Band function 
alone to the synthetic spectra, we adequately recovered 
the input parameters. However, the Cstat values were 
much lower than those resulting from the fit to the real 
data. It was often not possible to fit C+BB+PL to syn¬ 
thetic spectra corresponding to time intervals where the 
input /? values were high (i.e., ~ —2) because the fitting 
engine did not converge. In the cases it did converge and 
the resulting parameter values were “compatible” with 
the observed ones, the resulting Cstat values were much 
higher than those obtained when fitting C+BB+PL to 
the real data. Therefore, Band was always a better de¬ 
scription of the synthetic spectra. C+BB+PL was usu¬ 
ally an adequate model to fit the synthetic spectra cor¬ 
responding to time intervals with steep Band-function /3 
(< —5); in these intervals C+BB+PL and Band-alone 
fits led to similar Cstat values. However, C+BB+PL 
fits resulted in null fluxes for BB and PL, leaving only 
the C component. In addition, the C parameters were 
corresponding to the input parameters of the Band func¬ 
tion and not to those observed when fitting C+BB+PL 
to the real data. We then performed the same exercise 
but choosing C+BB+PL as the null hypothesis. When 
fitting the synthetic spectra with C+BB+PL, the fitting 
engine converged >80% of the time, allowing us to re¬ 
cover the input parameters. The resulting Cstat values 
were also similar to those obtained when fitting the real 
data though a little bit lower. When fitting a Band func¬ 
tion alone to these synthetic spectra, we recovered the 
parameters obtained when fitting the same function to 
the real data. In addition, the Cstat values were also 
similar to those resulting from the real data fits. This 
reinforces the idea that C+BB+PL is an overall better 
description of the time-resolved spectra for both GRBs. 
We note that this exercise could not be performed in 
all time intervals while leaving all parameters from the 
C+BB+PL model free; it was, however, possible to sim¬ 
ulate all time intervals where the additional PL clearly 
overpowers the two other components. 
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Fig. 15.— (a) Flux of the Band function (or of the cutoff power law) versus its E pea k when fitting Band-only (blue) or C+BB+PL (red) 
to the fine time intervals of GRB 080916C. The dashed red line corresponds to the best fit of the red data points with a power law; (b) 
Flux of the Band function (or of the cutoff power law) versus its E pea k when fitting Band-only (blue) or C+BB+PL (red) to the fine time 
intervals of GRB 090926A. The dashed red line corresponds to the best fit of the red data points with a power law; (c) Flux of the cutoff 
power law (i.e., non-thermal component, NT) versus its E pea k when C+BB or C+BB+PL are fitted to the time resolved data of four GRBs 
including short (green) and long ones (red, yellow and green). This is a refined version of the figure published in [Guiriec et all] ( |2013| ) 
showing the similarity of the F[^ T —E^^ k . relation for both short and long GRBs; (d) Luminosity of the cutoff power law versus its E pea k 

when corrected for both the redshift and the k-correction (i.e., L T— Ep^^l T relation) for the same sample of GRBs as in panel (c). The 
color dashed lines correspond to the best power law fit to the data of each individual GRB. The solid black line corresponds to the best 
power law fit to the whole sample simultaneously. 


7. Ff T -ENT AKii AND Lf T -ESf RELATIONS 

Figure pHTfi and b show the energy fluxes of the Band 
function (or C in the case of C+BB+PL) as a func¬ 
tion of Ep ea k when fitting either Band-alone (blue) 
or C+BB+PL (red) to the time-resolved spectra of 
GRB 080916C (Figure (1+) and GRB 090926A (Fig- 
We note that me correlation between the two 


ure 


15 r). 


quantities is much str onger with C+BB+PL . As reported 


for the first time in Guiriec et al. (2013), if this new 


Ff^-Epg^k ; relation is a specific property of the non- 
thermal component only (i.e., Band or C in the new 
multi-component model), it is expected that the fits of a 
Band function alone to GRBs that have additional spec¬ 
tral components will lead to large scatter in the data and 
in a weaker and biased relation or even in no correlation 
at all if the addi tional components are very intense. 

?NT_pNT 

peak,i 


In Figure 


15 


the Ff 


relations of several 


GRBs are overplotted. The results on short GRB 
12 0323A and long GR B 110721A were already reported 


Guiriec et ai. 1(2013) and preliminary results of GRBs 


080916C and 090926A using C+BB+PL were also shown 
in the same article. Here, the detailed analysis of GRBs 
080916C and 090926A confirms the results reported in 


Guiriec et al. (2013): in the observer frame, the Fj 
Epeak,i relations have similar slope for all GRBs whatever 
is the duration of the burst. In addition, when corrected 
for redshift, all GRBs lie along the same L^ T -E™ S *^ T 
relation—with a very limited scatter—suggesting a uni¬ 
versality for the phenomenon (see Figure 151). This new 
detailed a nalysis ti ghtened even more the relation pre¬ 


sented in Guiriec et al. (2013). 


The PL fit of the Lf T -Ep e e s *^ T relation for each indi¬ 
vidual GRB result is shown in colored dashed lines: 
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Fig. 16.— Distribution of the spectral parameters resulting from 
the fine-time analysis of both GRBs 080916C and 090926A using a 
Band function alone. The filled blue histograms correspond to the 
distribution of the best parameters (i.e., central values) and the 
dashed red and green histograms corresponds to the distributions 
of the l-cr lower and upper limits, respectively. The vertical black 
dashed line at -2/3 and at -3/2 correspond to the limits above which 
the values of a are incompatible with pure synchrotron emission 
from electrons in both the slow and fast cooling regimes, and with 
synchrotron emission from electrons in the fast cooling regime only, 
respectively. 


Fig. 17.— Distribution of the spectral parameters resulting from 
the fine-time analysis of both GRBs 080916C and 090926A using 
a the C+BB+PL model. The filled blue histograms correspond to 
the distribution of the best parameters (i.e., central values) and the 
dashed red and green histograms corresponds to the distributions 
of the l-cr lower and upper limits, respectively. The vertical black 
dashed line at -2/3 and at -3/2 correspond to the limits above which 
the values of a are incompatible with pure synchrotron emission 
from electrons in both the slow and fast cooling regimes, and with 
synchrotron emission from electrons in the fast cooling regime only, 
respectively. 
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The simultaneous fit of all the GRBs with a PL results 
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The strong compatibility of the results for all the tested 
GRBs reinforces the possibility to eventually use this re¬ 
lation as a redshift estimator. 

We must note that the initial rising part of the first 
pulse of a burst usually does not follow the typical F+ - 
®£oak i relation as already reported in 


Guiriec et al. 


(2013), but an anti-correlation is observed between F{ N 1 


and E pdk,i instead. 


8 . A NEW MODEL FOR GRB PROMPT EMISSION 

In Section [6j we reported the constancy in the values 
of the a and the indices of the additional PL when fit¬ 
ting C+BB+PL to the time-resolved data of both GRBs. 
Not only did these values remain mostly constant dur¬ 
ing each burst (see panels (b) and (e) of Figure [9] and 
10 ), but they are also nearly the same for both events 
{see panels (b) and (d) of Figure [l2| and fl4|. Figures[l6] 
and [TT] show the Band-only and (J+BB+PL combinea 
parameter distributions, respectively, for GRBs 080916C 
and 090 926A . This reinforces the conclusions reported in 
Section |6.2[ the Band-only fits result in a narrow dis- 


tribution of the E pe ak values around 300 keV and in a 
broader distribution for a, with a large fraction of the 
latter values > —0.7. In contrast, the E peak values for 
C+BB+PL are much more spread with a maximum of 
the distribution around 500 keV, the values of a are nar¬ 
rowly peaked around —0.7 and the a l-cr lower limits 
are all bellow —0.7 besides a few cases, and the values 
of the additional PL index are clustered between —1 and 
—2 with a distribution peak around —1.5. 

Since the slopes of the low-energy PL of C and of the 
additional PL in the C+BB+PL model do not change 
much with time and remain similar in the two bursts, 
those two parameters can be frozen to their typical val¬ 
ues in the fits, reducing the complexity of the C+BB+PL 
model by two dof. This new model, therefore, which we 
will call C+BB+PL 5 params hereafter, has only 5 param¬ 
eters . C+BB+PL 5 params has only one parameter more 
than Band-alone, which makes the two models very com¬ 
petitive from a statistical point of view. 
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Fig. 18.— Energy flux evolution between 8 keV and 100 MeV for (al) GRB 080916C and (bl) GRB 090926A in the context of the 
C+BB+PLsparams model. Panels (a2) and (b2) of each figure show the contribution of each component to the total energy flux of 
GRBs 080916C and 090926A, respectively. For clarity, no uncertainty is displayed for the total energy flux (black line). Because of its 
narrow shape, the contribution of the BB component to the total 8 keV-100 MeV flux is ver y lim it ed; h owever, it has an important relative 
contribution to the total emission at the BB spectral peak around 100-200 keV (see Figure |A3| & |A4| of Appendix | e|). 


We fitted C+BB+PL 5 params to the time-re solve d spec¬ 
tra and the results are reported in Table |A5| of Ap¬ 
pendix 0 & Figure |A3| of Appendix [El and Table |A 6 | 
of Appendix [C] & Figure |A4| of Appendix [E] for GRBs 
080916C and 090926A, respectively. As expected, 


C+BB+PL 5 params results in similar Cstat values as 
C+BB+PL except for a few cases where the additional 
PL slope appears to be a bit steeper than —1.5. However, 
this can be easily explained by the presence of a cutoff in 
the additional PL. Indeed, the C+BB+C 6 par ams fits (i.e., 
C+BB+PL 5 params with a exponential cutoff in the addi¬ 
tional PL) seem to indicate the presence of a cutoff whose 
energy evolves smoothly with time. The cutoff could be 
hidden by the C component at times, and it could be at 
much higher energies at others. This result may ques¬ 
tion the use of the additional PL cutoff to estimate th e 
Lorentz factor as proposed in Ackermann et al. (2011). 
Instead, such an evolution of the cutoff energy may be 
a signature of the emission mechanism(s) producing this 
component and not an effect of the 7-7 opacity. 

It is very important to note that the number of time 
intervals in which the three components cannot be si¬ 
multaneously fitted is very limited for C+BB+PL 5 params 
compared to C+BB+PL. Even more interestingly, 
C+BB+PL 5params results indicate that the additional PL 
is present from the very beginning of the prompt emission 


in both GRBs. Moreover, the additional PL is the most 
intense component at very early times. This is supported 
by the very low values of a (i.e., < — 1 ) obtained in the 
first time intervals when fitting C+BB+PL with all the 
par amet ers free (see Figure [9jo fe |1T)|) ). Indeed, in Sec- 
tion |6.1| we reported that we could not fit the C+BB+PL 
model at early times, but only C+BB. Intriguingly, in 
these time intervals, the values of a were significantly 
lower than average and closer to the typical index val¬ 
ues of the additional PL. This can be explained if we 
consider the fit resulting from C+BB+PL 5params as the 
“true” model. With C+BB+PL 5params , the additional 
PL overpowers the other two components, and when fit¬ 
ting C+BB to the data, the value of a is a hybrid (i.e., 
~ —1.1) between the real value of the C component (i.e., 
~ —0.7) and the index of the additional PL that is not 
included in the fit (i.e., ~ —1.5). The same conclusions 
can be drawn for GRB 090926A to explain the low val¬ 
ues of a obtained when fitting C+BB at late times (see 
Figure [lOb). Indeed, in the C+BB+PL 5params scenario, 
the contribution coming from the additional PL is the 
most intense at late times (see panels (d) of Figure [ 8 ]). 

The reconstructed photon light curves resulting from 
C+BB+PL 5params are presented in panels (d) of Figure]?] 
and [ 8 |for GRBs 080916C and 090926A, respectively. For 
both bursts, these light curves are in very good agree- 
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ment with the observed count light curves, in particular 
the LLE light curve between 20 and 100 MeV. Indeed, 
although we do not use the LAT-LLE data in our fits, 
the extrapolations of our model to higher energies repro¬ 
duce very well the observations with a smooth evolution 
of the additional PL flux with time. The resemblance of 
the observed 20-100 MeV light curves with the recon¬ 
structed one is particularly striking for GRB 080916C: 
not only the intensity peak is well recovered, but the 
long tails before and after it are also very adequately 
reproduced. Figure [18] shows the evolution of the en¬ 
ergy content of each component with time (panels (al) 
and (bl)) as well as their relative contributions to the 
total energy flux (panels (a2) and (b2)) between 8 keV 
and 100 MeV in the observer frame. The Band function 
contribution to the total energy flux dominates at early 
time while the additional PL contribution is dominant 
at later time; however, the contribution of the additional 
PL is also important at very early times during the first 
instant of the bursts. Because of the narrow shape of 
the thermal-like component, its contribution to the total 
energy flux is limited although it may reach few tens of 
percent in some time intervals in the best case scenario. 

A larger sample of bursts will help to better estimate 
the typical values of a and of the additional PL index. 
While for the two GRBs studied in this article, as well as 
for many others, a value of a around —0.7 seems to be 
a good estimate, f or others the value of a seems to peak 
around —1.2 (e.g., Guiriec et al. 2013) ; both values seem 
to be typical ones for a. 

The F^ T -Ep^ k; relation can be used to better 


constrain the C+BB+PL 5 p arams model and to even 
further decrease its complexity. Indeed, in the 
C+BB+PL 5 params model the values of a, “/3” (i.e., an ex¬ 
ponential cutoff for C), and the index of the additional 
PL are frozen, therefore, the Ff T H5pg(( kii relation can 
be reduced to a correlation between the amplitude of 
the Band function (or C), A, and its peak energy E pea fc. 
This relation between A and Epeakj can then be deter¬ 
mined with an initial time-resolved analysis of the data 
for each burst with the C+BB-|-PL 5 p arams model, and 
subsequently be used as input to the C+BB-|-PL 5 p arams 
model to tighten even more the fit results. 

Moreover, if the Lf T -Ep® s £ |( k r ^ T relation is confirmed and 
shown to hold for a large sample of bursts, then the re¬ 
lation between A and E pea k may be firmly determined 
beforehand for GRBs with known redshift and would al¬ 
low us to constrain the C+BB+PL 5 paiams model from 
the beginning; therefore, the C+BB-|-PL 5 p arams model 
would then have only 4 independent parameters (i.e., 
C+BB+PL 4 params ) like the Band function alone. For 
GRBs with unknown redshift, z would be an additional 
parameter (i.e., 5 independent parameters and a con¬ 
straint between A and E pea k) and it may be possible 
to evaluate the redshift of the burst from the fits of 
C+BB+PL 4 paramS!Z to each time-interval. Finally, it 
could also be possible to fit directly the cosmological pa¬ 
rameters to the time-resolved spectra in the case of GRBs 
with known redshifts (i.e., C+BB+PL 4 p aramSi n M ,n A ). Al¬ 
though preliminary results look promising, a comprehen¬ 
sive study is required to validate or not those relations 
prior to any possible use; those results will be presented 
in a follow up article. 


9. INTERPRETATION 

There are several theoretical possibilities for the phys¬ 
ical origin of each of the three components identified in 
the spectra of GRB 080916C and GRB 090926A. How¬ 
ever the general framework of the discussion is well es¬ 
tablished: (i) GRBs are produced by ultra-relativistic 
outflows and (ii) the prompt emission has an internal 
origin. More precisel y, (i) is required to avoid a strong 
7-7 annihila tion (e.g. |Baring & Hardingj[1997| |Lithwick 
fe SaripOOT ). A detailed calculation leads to constraints 
for the Lorentz factor and the radius of the emission site 


to allow the escape of photons at a given energy E (Gra- 
not et al.|2008 Hascoet et al.|20f2 ). For the Band (or G) 
and BB components, this calculation does not provide a 
strong constraint on the emission site, that can be any¬ 
where from the photosphere to the deceleration radius. 
On the other hand, as the PL component is extending 
at least up to 100 MeV, it must be produced above the 
photosphere. Condition (ii) is obtained from the study 
of the variability of the prompt emission, which cannot 
be reproduced by the external shock during the phase 
when the_outflow is decelerated by the ambient medium 
(e.g. |Sari fe P ir an|199 7). This implies that the Band (or 
C) and BB components must be of internal origin. For 
the PL, our analysis clearly shows that variability below 
the one-second timescale is also present (see panels (c) 
and (d) of Figure [7]&[ 8 j red curves), indicating that the 
origin is at least partially internal. On the other hand, 
the evolution of this component, which becomes brighter 
at the end of the burst suggests that a link with the 
deceleration of the outflow should also be considered. 


9.1. Origin of the Band (or C) and BB components 

In the first class of models—dissipative photospheres— 
the bulk of the keV-MeV prompt emission is released 


at the photosphere of the relativistic outflow (Goodman 
1986 Paczynski 1986), i.e. at a low radius compared to 


other models. The spectrum can have a complex shape, 
far from a thermal spectrum, due to a sub-photospheric 
dissipation process that accelerates electrons, allowing 
Comptonization at high-energies, and, if the outflow is 
magnetized, synchrotron radiation at low-energies (e.g. 
Rees fe Meszdros 120051 [Pe’er et al. 2006| |Beloborodov| 
2010| ?urm et al.|2011 ). In this scenario, the Band (or G) 


component and the BB component should be interpreted 
as a single component with a complex shape. The physics 
governing the evolution of this shape is related not only 
to the well-understood physics of the photosphere, but 
also to the nature of the dissipative sub-photospheric 
process, which is not well constrained. It is then too 
early to conclude if this scenario can reproduce the ob- 


served spectral evolution (see the discussions in Vurm 
et al. ||2013l |Asano & Meszaros||2013[ |Gill & Thompson 

20T4F. 


In the second class of models the photosphere is not 
strongly affected by sub-photospheric dissipation. The 
spectrum of the photosphere then remains thermal, with 
a Planck spectrum slightly modified at low-energy due to 
the pec uliar geometry of the photosphere in a relativistic 


outflow ( Bel oborodov 2011 Lundman et al. 2013 Deng & 
Zhang |2014 ). This provides an obvious interpretation of 
the BB component in our a nalysis. However, as pointe d 
out in previous cases (e.g. Guiriec et al. 2011a 2013), 
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the photospheric emission in our two bursts appears to 
be less hot and bright than the prediction of the stan¬ 
dard fireball model where the acceleration of the ultra- 


the dynamical ti mescale but comparable to the electron 


relativistic outflow is purely therma l ( Meszaros fe Ree s 
2000 Daigne & Mochkovitch 2002 iNakar et al. 2005). 


This implies that the o utflow is strongly magnetize d close 
to the c entral source (Zhang & Pe’er 2009[ Hascoet et 


al. 2013). This important conclusion of previous analy- 


sis of bright Fermi bursts remains unchanged with the 
three-components spectral analysis of GRBs 080916C 
and 090926A. On the other hand, our analysis tends to 
show that the BB component is always detected in bright 
bursts where an accurate spectral analysis is possible, 
whereas GRB 080916C was often presented in previous 
studies as a g ood example of burs t without photospheric 
emission (e.g |Abd o et ah 2009a Zhang & Pe’er 2009 
Zhang et al.|20TTc |. it this is also true tor the whole GR1 


population where a multi-component spectral analysis is 
not possible, this implies that all GRBs are associated 
to magnetized outflows, with initia lly only 1 to 10% of 
the energy in thermal form (Hascoet et al.|[2013l, allow¬ 


ing for a residual photospheric emission at the observed 
level, as in the two cases presented here. 

In this scenario, the non-tliermal emission correspond¬ 
ing to the Band (or C) component is produced in the 
optically thin regime, above the photosphere. It is as¬ 
sociated to the emission of relativistic electrons, accel- 


erated either in internal shocks (IRees & Meszaros|1994 

Kobayashi et al. 

1997 

Daigne & Mochkovitch 1998 

) or 

due to magnetic reconnection (e.g. Lyutikov &; BlandJ 

ord 


occur if the conversion of the magnetic energy into kinetic 
energy is efficient during the acceleration of the outflow, 
leading to a low magnetization a < 0.1 at large radius. 
The second possibility corresponds to the opposite case, 
where the acceleration is inefficient and the magnetiza¬ 
tion is still high at large radius, a > 1. In both cases, 
the dominant radiative process should be synchrotron, 
as the synchrotron self-Compton (SSC) predicts an ad¬ 
ditional strong component either in the optical-UV-soft 
X-ray range, or above 100 MeV, in contradiction with 
observations (Piran et al. I2009[ Bosnjak et al.|[2009 ). 

Our analysis confirms that when considering only the 
Band (or C) and BB components, this scenario is in 
good agreement with observations. In the case of GRB 
080916C, the low-energy photon index a of the non- 
thermal spectrum remains below -1 for almost the whole 
evolution, which is compatible with the expected fast¬ 
cooling regime, when taking into account the effect on 


the Klein-Nishina regime 

iDerishev et al. 2001 Bosnjak 

[etaL|[2009 

Wang et al. 

2009 

Daigne et al. 

20IT). In 


26A, a is larger, usually in the -0.7 - -1 range. 
This is of course fully compatible with the synchrotron 
radiation in slow cooling regime. However, this inter¬ 
pretation is not satisfactory as the radiative efficiency is 
very low in this regime, leading to an energy crisis for 
a bright burst like GRB 090926A. Large photon indices 
close to -0.7 can however be obtained if electrons are only 
marginally fast cooling (the cooling break being close to 
the peak energy, the intermediate sl ope -1.5 is masked 


and the slope —2/3 is ob served (e.g. Daigne et al. 2011 
Beniamini & Piran|2013) or if the magnetic held in the 
shocked region is decaying on a timescale smaller than 


cooling timescale (Derishev 

2007 

Lemoine 2013 

Uhm & 

Zhang; 2014 

Zhao et al. 2014|). 


the reconnection scenario can reproduce the observed 
lightcurves. In the case of internal shocks, detailed sim¬ 
ulations show in addition that the spectral evolution can 
also be reproduced, inclu ding the correlation betwe en the 
peak energy and the flux ( Bosnjak fe Daigne|2014 ). Such 


spectral calculations are not available yet tor reconnec- 
tion mo dels. The generi c calculation of synchrotron radi¬ 


ation by Uhm fe Zhang (2014) shows that the continuous 
injection of electrons by reconnection events as expected 
in the ICMART model could lead to the c orrect spec¬ 
tral slope (see also Beniamini & Piran [2014 ). It remains 
to test that the addition of the contributions from sev¬ 
eral simultaneous emission sites in the outflow would not 
broaden the spectrum too much. Detailed calculations 
of the photospheric emission in a variable outflow pre¬ 
dicts that evolution is also expected for the BB compo¬ 
nent. However, due to the weakness of this component, 
we probably detect only the brightest peaks, which does 
not allow u s to eva luate the real variations of the tem¬ 
perature (Hascoet et al.j|2013). 


Overall, both the spectral parameters and the spec¬ 
tral evolution observed in GRBs 080916C and 090926A 
are compatible with the scenario where these two bursts 
are associated with an ultra-relativistic outflow, which 
is initially highly-magnetized, producing a weak thermal 
emission at the photosphere (BB component) and where 
dissipation in the optically thin regime (shocks or re¬ 
connection) leads to synchrotron radiation of accelerated 
electrons (Band component). 

Unfortunately, the situation regarding the synchrotron 
origin of the Band (or C) component becomes more com¬ 
plex in the three-component analysis, as the addition 
of the PL leads to larger photon index a, usually close 
to —2/3, which could be explained by a modified syn¬ 
chrotron radiation as discussed above. However, we point 
out that our analysis shows that a third component, dom¬ 
inant at low (< 10 keV) and high (> 100 MeV), is re¬ 
quired to better fit the observed GRB spectrum, but that 
a power-law is only the simplest possible assumption for 
the shape of this component. The real shape is difficult 
to constrain, especially because this third component is 
masked in the intermediate MeV range. From a theo¬ 
retical point of view, it is difficult to understand which 
process could produce a single power-law extending over 
at least five decades in energy. It is more natural (see 
below) to expect that the spectral shape of the third 
component is in reality more complex. A small change 
of the shape at low-energy would allow us to reconcile 
the three-components analysis with the synchrotron ra¬ 
diation for the dominant component in the MeV range. 

9.2. Origin of the PL component 

The third, PL, component is clearly the most puzzling. 
It has been suggested that it cou ld be entirely of external 
origin (e.g. Abdo et al. 2009b|) However, our analysis 
shows that it is present inPboth GRBs since the very 
beginning, when the deceleration of the outflow is still 
negligible. In addition, it also shows clearly that time- 
variability is present, on the 1 s timescale and proba¬ 
bly below. This suggests an internal origin, at least at 
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early times. One possibility discussed in the literature 
is a hadronic component (see e.g. Abdo et al. 2009b). 
However its low efficiency leads to an energy cr isis tor 


bright bursts as GRBs 080916C and 090926A (Asano 


et al. 2009; Asano & Meszaros 


2012). This leaves only 


electrons as a reasonable candidate tor this emission. If 
the Band (or C) component is due to synchrotron ra¬ 
diation, then the radiating electrons must have initially 
very large Lorentz factors (> 10 4 ) to produce MeV pho¬ 
tons, which limits the inverse Compton scatterings due 
to Klein-Nishina corrections. However, the evolution of 
the Band (or C) component shows large variations of the 
spectral peak, allowing it to be less deep in the Klein- 
Nishina regime when the peak energy is lower, and even 
to enter the much more efficient Thomson regime. There¬ 
fore, inverse Compton emission is expected at high en¬ 
ergy and is not expected to peak exactly at the same 
time as the synchrotron radiation, in agreement with the 
observed lightcurves of the Band (or C) and PL compo¬ 
nents in GRB 080916C and 090926A. Calculations made 
in the context of internal shocks predict a fiat spectrum 
in agreement with the observed index of the PL close to 
-2 (|Bosnjak et al.||2009| |Bosnjak & Daigne||2014 1. 

The inverse Compton scenario does not explain why 
the PL is extending down to the keV range. There are 
however some mechanisms that could correlate a compo¬ 
nent in the keV range and one above 100 MeV, leading to 
a more complex shape than a simple power-law: (i) one 
possibility is associated to the 7-7 annihilation at high 
energies. The signature of this process is a time-evolving 
attenuation leading to a broken power-law at high-energy 
followed by a strong, close to exponential, cutoff. Our 
study cannot confirm this signature in GRBs 080916C 
and 090926A but the fact that a C in replacement of the 
PL is not rejected by the analysis may indicate that 7 - 
7 annihilation is indeed present in the 100 MeV-10 GeV 
range. Then a cascade of pairs is produced, which rad iate 


at lower energy b y synchrotron radiation (see e.g. Asano 


& Meszaros 2011). It should be tested if such an emission 


could emerge in the keV range ; (ii) a second possibility is 
to assume that less relativistic electrons are also present 
in the flow, either because the acceleration process does 
not accelerate all of them, or because the dissipation pro¬ 
cess is weaker in some regions (residual internal shocks 
with small relative Lorentz factors for instance.) Then, 
the synchrotron radiation of these electrons peaks in the 
keV range, and inverse Compton scatterings occur in the 
Thomson regime, producing high-energy photons. This 
could explain the third component fitted by a PL, but 
with a more complex spectrum. A detailed modeling, 
which is out the scope of this paper, will be needed to 
test these ideas. Note that in the inverse Compton sce¬ 
nario for the PL component, the possible cutoff in the 
PL found in some time bins could be an evidence of the 
intrinsic curvature of this component close to its peak 
energy, rather than a signature of the 7-7 annihilation. 

Finally, the fact that the PL component becomes dom¬ 
inant at the end of both bursts suggests that a new 
component, of external origin, is emerging at late times, 
when deceleration starts. The physical origin would then 
be the same as for the long-lasting emission revealed by 
Fermi- LAT, which is observed both in GRBs 080916C 
and GRB 090926A. This origin is still debated, but syn¬ 
chrotron and inverse Compton radiation in the exter- 


nal shock are strong candidates (Zou et al. 2009 Kumar 

& Barniol Duran 2009) Ghiseiiini et al. 

2010 ; [Kumar & 

Barniol Duran 20101; 

Lemoine 2013f 

Be 

[oborodov et al. 

2014| Vurm et al 

201 

~i'l . Interesting 

y, an early tempo- 


of GRB 090926A (Ackermann et al. 2011), which rein¬ 
forces the scenario where the high-energy emission is ini¬ 
tially dominated by a component of internal origin, with 
a transition to an external origin once the deceleration 
starts. 


10. CONCLUSION 

We summarize our results as follows: 


• GRB prompt emission spectra are more complex 
than the shape resulting from the fit to a Band 
function alone. An adequate option to capture the 
complex shape of GRB prompt emission is a com¬ 
bination of at least three separate spectral com¬ 
ponents whose intensity and relative contribution 
to the total energy flux evolve with time, (i) The 
three components are not always present in all 
GRBs; or (ii) in some GRBs certain components 
strongly overpower the others making their identifi¬ 
cation more difficult. Figure [ 6 ] shows an example of 
the three-component model (T7e., ( 1 ) a non-thermal 
Band function with a high-energy power law index 
<-3 and which is statistically equivalent to a cut¬ 
off power law, ( 2 ) a thermal-like component ade¬ 
quately approximated with a black body compo¬ 
nent, and (3) a non-thermal additional power law 
with or without cutoff) in GRB 090926A. 


• The new three-component model is composed of (i) 
a smoothly broken PL with a constant low-energy 
ind ex of ~-0.7 (o r ~-1.2 as for GRB 120323A— 


Guiriec et al. ( 2013|)) and a high-energy slope 


see 

compatible with an exponential cutoff based on the 
quality of our current data (i.e., C), (ii) a thermal- 
like component adequately approximated with a 
BB component, and (iii) another PL component 
with a fixed index of ~-1.5, whose cutoff energy is 
often challenging to identify. However, the cutoff 
would be most of the time below 100 MeV. When 
a cutoff is used, the index of the additional PL gets 
a bit closer to -1. It is important to note that this 
paper aims to establish the existence of at least 
three distinct spectral components, but the spec¬ 
tral shapes of the various components will be ex¬ 
plored in greater detail in a following article. 


• In the C+BB+PL model, (i) the C component can 
be interpreted as synchrotron radiation in an op¬ 
tically thin region above the photosphere, either 
from internal shocks or magnetic field dissipation 
regions, (ii) the BB component can be interpreted 
as the photosphere emission of a magnetized rela¬ 
tivistic outflow and (iii) the extra PL extending to 
high energies likely has an inverse Compton origin 
of some sort even though its extension to a much 
lower energy remains a mystery. 


• We succeed in reducing the number of parameters 
in the three-component model (i.e., C+BB+PL) 
from 7 down to 5 (i.e., C+BB+PL 5 params ) still 
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keeping the same quality for the fits; with 5 free 
parameters C+BB+PL 5 params is statistically com¬ 
petitive with the 4 parameters of the Band func¬ 
tion. Additional constraints to the model will be 
added if the FBand,i~Bpeak.i and Lg anc i.i~Fpeak,i re¬ 
lations are confirmed on a large sample of GRBs; 
This will result in an even simpler model. 


• The C and BB components are intense at early 
times in a burst. The additional PL may kick off 
at very early times and even before the two other 
components and last much longer at low energies. 
The intensity peak of the additional PL seems to be 
often related to sharp and bright structures present 
in the light curves at all energies from 8 keV to tens 
of MeV and even GeV. The additional PL is the 
most intense component at late times at low ener¬ 
gies (tens of keV). This additional PL may connect 
smoothly with the GeV LAT emission observed for 
thousands of seconds after the prompt phase and 
contemporaneously with the X-ray and optical af¬ 
terglow emissions. Because the BB component is 
narrow, its contribution to the total energy flux 
between 8 keV and 100 MeV is very limited com¬ 
pared to the two other components (see figure |18[) ; 
however, the BB has an important relative contri¬ 
bution to the total emission at its s p ectra l peak 
around 100-200 keV (see Figure |A3| & |A4| of Ap¬ 
pendix [e]). 


• From the data set we used, it is difficult to con¬ 
clude if the additional PL is a single component. 
However, it is clear that there is a simultaneous 
evolution of the low- and high-energy fluxes which 
seems to propagate all across the spectrum. 


• There is a strong correlation between the energy 
flux of the non-thermal component—adequately 
fitted with either Band or C—and its E pea k within 
each burst. Taking into account the various com¬ 
ponents in the fit process reduces the scatter of this 
relation. When fitted with a PL, the Ff^-Ep^j 
relations have similar slopes for all GRBs. This re¬ 
lation may be useful to understand the mechanism 
responsible for the non-thermal prompt emission. 

• In the central engine rest frame, the F^-E^k . 
relation points towards a universal L^ T -Ep e e s *^ T 
relation. If confirmed, this relation could be used 
to estimate GRB redshifts as well as to constrain 
the cosmological parameters. 

Beyond the purpose of the current article, this analysis 
using multiple spectral components may have a broader 
impact on other analyses: 

• Since the cutoff in the additional PL seems to 
evolve with time—and most of the time bellow 100 
MeV—and since it is sometimes hidden by other 
components, it may indicate that this cutoff is a 
property of the emission mechanism(s) producing 
the additional PL and not the result of 7-7 opacity. 
If it is the case, it is not correct to estimate the bulk 
Lorentz factor of the jet based on the additional PL 
cutoff energy. 


• The presence of multiple components evolving in¬ 
dependently with time and exhibiting highly vari¬ 
able relative fluxes may strongly bias the conclu¬ 
sion resulting from spectral lag analysis. Indeed, 
lags may be identified between high intensity struc¬ 
tures belonging to different components. 

• Such multi-component spectral analysis may help 
in the context of analysis conducted to constrain 
Lorentz invariance. Indeed, such analysis considers 
that photons are emitted from the same region and 
at the same time. By identifying the various com¬ 
ponents we may improve the association between 
low- and high-energy photons, as it would be the 
case e.g., for GRB 080916C. Based on the count 
light curves of Figure [ZK we may be tempted to 
consider the association of the peak above 20 MeV 
around 6 s with peaks at lower energy happening 
at earlier (or later) times, especially if we do not 
have the knowledge of the low-energy light curve 
below 20 keV. With the multiple component anal¬ 
ysis, we see that the peak at high-energy is directly 
correlated to the one at low-energy. Therefore we 
can increase the accuracy of the Lorentz invariance 
measurement by reducing the uncertainties on the 
time dispersion. 
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APPENDIX 

A. CROSS-CALIBRATION 
A.l. Time-Integrated Spectral Analysis 

We find that a normalization correction is not required between the selected Nal and BGO detectors in both GRBs 
for all fitted models, i.e., all EAC factors are compatible with unit y wit hi n 1- q uncertainties corresponding to ~5% 
and ~10% for GRBs 080916C and 090926A, respectively (see Tables A1 & A2). When we fit the GBM and LAT data 
of GRB 080916C simultaneously, the three LAT data sets (i.e., LLE , LATVFront and LAT-Back) are consistent with 


the GBM ones for all the tested models except B+PL (see Table Al |. For GRB 090926C, the GBM and LAT data are 
only comp atib le for the B+C, B+BB+B2 or C+BB+B2 fits, but large corrections are required with all other models 
(see Table A2). In all cases, however, where EACs are implemented, their uncertainties are large, making it difficult 
to precisely quantify the intercalibration of the two instruments. In addition we note that the inclusion of EACs does 
n ot af fect significantly the values of the spectral parameters of both GRBs in any of the models used (see Tables Al 
& A2|. We decided, therefore, to not consider any EAC corrections in the time-integrated spectral analysis. 


A. 2 . Coarse-Time Spectral Analysis 


GRB 080916C: 


We find that all GBM detec tors are consistent with each other within 5 to 10% uncertainties for all tested models 
and in all intervals (see Table A3). Simultaneous fits of GBM and LLE data are consistent with fits of GBM data 


alone even in the second time interval, where the high-energy emission is most intense (see Table A3). However, since 
the intensities of the BGO and LLE signals are too weak to constrain with enough accuracy the EAC factors between 
GBM and LLE data, a detailed study of EAC factors here is meaningless (we also note that we found a small effect 
of EAC factors on the Cstat values). 

We conclude that all GBM detectors and LLE data have a consistent relative calibration and we note that the spectral 
parameters are similar with or without EACs. Therefore, in the following we will only describe spectra resulting from 
fits that do not include any EAC. 

GRB 090926A: 

In the case of GRB 090926A, we used both BGO detectors and determined that they are consistent within ~10% to 
~20% uncertainties (see Table A4). Beyond T 0 +lls, however, the signal in the BGO detectors becomes weak and it 
is difficult to obtain good constraints on th eir EAC factors. All the Nal detectors are also consistent with each other 
within 10-20% uncertainties (see Table A4). We identified discrepancies of the order of a few percent between the Nals 
and the BGOs supported by the change in Cstat values when EAC factors were included in the fits, although we could 
not quantify these precisely (see Table A4|. Despite the effect on the Cstat values, EAC factors do not significantly 
impact the spectral parameters resulting from the fits using the various models nor the Cstat variations (i.e., ACstat) 
from model to model. Therefore, we conclude that, similarly to GRB 080916C, we do not need to include EAC factors 
between the various GBM detectors. 

Significant EAC corrections are required be twee n the GBM and the LLE data for models that do not include any 
additional PL to the Band function (see Table A4). This is easily explained by the fact that the high-energy emission 
observed with the LAT is clearly incompatible with an extrapolation of the high-energy PL of the Band function above 
several tens of MeV. Therefore, to account for these discrepancies, we can either include an EAC factor or an additional 
component at high-energy; the latter option looks more satisfactory because such strong discrepancies between GBM 
and LAT are not expected. In addition, since the additional PL to the Band function can be identified from GBM only 
data, it is quite natural to consider that the discrepancies observed between GBM and LLE data—when no additional 
PL is used—are not only due to calibration issues but also depend significantly on the accuracy/goodness of the model 
fitted to the data. For the models that include an additional PL, GBM and LLE data look consistent with EAC factors 
close to unity; however, the EAC uncertainties are very large making any correction unreliable. We conclude from this 
analysis that the GBM and LLE data are in agreement but that correction factors between the two data sets cannot 
be accurately determined; however, the results are not significantly affected, so in the following, we only present the 
spectra obtained without EAC between the GBM and the LLE data. 

B. MODEL COMPARISON PROCEDURE 

In every time intervals of the time-integrated, coarse-time and fine-time spectral analysis, we performed multiple 
sets of Monte Carlo simulations (i) to compare all the tested models one to one, (ii) to estimate the goodness of fit 
of every tested model, and (iii) to estimate the trustworthiness of the measured spectral parameters for all the tested 
models. To generate the synthetic spectra of our simulati on sets, we followed exactly the same procedure as already 
described in Guiriec et al. (2011a) and in the Appendix of Guiriec et al. (2013). 

In the article, we will use the statistical terminology and call ‘‘best modeE and “best fit” the model and the fit, 
respectively, that have the lowest Cstat value. 
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To assess which model is a “better description of the real data” set when comparing models Mi and M 2 , we performed 
ten statistical tests addressing the ten following questions (shortened for simplicity). In all tests we assume that both 
signal and background fluctuations are following Poisson statistics. 

• Test 1: estimate the odds to obtain a ACstat value between Mj and M 2 higher than the observed one if the true 
underlying model is Mj — with real fit parameters” => This is the typical “likelihood-ratio test” or “log-likelihood 
statistic test.” 

• Test 2: estimate the odds to obtain the observed parameters of M 2 when fitting M 2 to the data if the true 
underlying model is Mi—with real fit parameters 

• Test 3: estimate the odds to obtain the observed parameters of Mi when fitting Mi to the data if the true 
underlying model is Mi—with real fit parameters => This also serves to estimate the trustworthiness of the 
measured parameters. 

• Test 4: estimate the odds to obtain the specific Cstat value when fitting Mi to the data if the true underlying 
model is Mi—with real fit parameters => This is the typical goodness of fit test. 

• Test 5: estimate the odds to obtain the specific Cstat value when fitting M 2 to the data if the true underlying 
model is Mi—with real fit parameters 

• Test 6: estimate the odds to obtain a ACstat value between Mi and M 2 higher than the observed one if the true 
underlying model is M 2 -with real fit => This is the typical “likelihood-ratio test” or “log-likelihood statistic 
test.” 

• Test 7: estimate the odds to obtain the observed parameters of M 2 when fitting M 2 to the data if the true 
underlying model is M 2 —with real fit parameters => This also estimates the trustworthiness of the measured 
parameters. 

• Test 8: estimate the odds to obtain the observed parameters of Mi when fitting Mi to the data if the true 
underlying model is M 2 —with real fit parameters 

• Test 9: estimate the odds to obtain the specific Cstat value when fitting Mi to the data if the true underlying 
model is M 2 —with real fit parameters 

• Test 10: estimate the odds to obtain the specific Cstat value when fitting M 2 to the data if the true underlying 
model is M 2 —with real fit parameters => This is the typical goodness of fit test. 

We consider that M 2 is a better description of the data than Mi if: 

1. Test 1 is conclusive (i.e., p-value<~10~ 5 —see footnot^Jfor the definition of p-value) and if Test 6 is inconclusive 
(i.e., p-value>>10 -5 .) 

2. Tests 3, 7, 8, 9 and 10 are inconclusive (i.e., p-value>>10 -5 ) and Tests 4 and 5 are conclusive (i.e., p-value<<<.) 
In this case, Test 2 is also usually conclusive. 

In addition, if M 2 is a “better description of the data” than Mi because of criterion 1, we will also say that M 2 is 
“statistically significantly better” than Mi. 

We note that as of today, the GBM instrument team does not anticipate any evidence for global strong instrumental 
systematic effects that will not be accounted in the instrument response function yet; however, we cannot discard the 
fact that possible calibration problems discovered at later time may impact the results of the current analysis as well 
as all the analysis performed with the traditional Band function alone. 


Hereafter we introduce two typical cases: 


1. Comparison of the B+BB and C+BB+PL fits to the GBM data alone in the time interval lasting from To+4.3 s 

to To+7.5 s of the coarse-time analysis of CRB 080916C: 


20 For Tests 1, 4, 5, 6, 9 and 10, p-value is the probability to 
obtain a value for a specific quantity beyond the actually measured 
value of the same quantity from the real data (i.e., one-tail proba¬ 
bility.) For Tests 2, 3, 7 and 8, p-value is the probability to obtain 


an absolute difference between the value of a parameter resulting 
from the fit and the most probable value of the distribution for 
the same parameter beyond the l—a confidence region around the 
most probable value of the distribution (i.e., two-tail probability.) 
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In the following, we will consider that B+BB and C+BB+PL correspond, respectively, to the Mi and M 2 models 
from the description of the various tests presented earlier in this Appendix. B+BB and C+BB+PL have 6 and 7 free 
parameters, respectively, but they are not nested models. 

The B+BB and C+BB+PL fits to the real data result in Cstat values of 508 and 495, respectively. In this time 
interval, C+BB is clearly a very bad fit to the data (i.e., extremely high Cstat value) and it is, therefore, not considered 
here. Indeed, when fitting B+BB to this time interval /3 is >-2.28 and clearly incompatible with the exponential cutoff 
of a C component. 

We first present the results of the typical log-likelihood statistic test using B+BB with the parameter values resulting 
from the B+BB fit to the real data as the null-hypothesis (i.e., Test 1.) Despite the fact that C+BB+PL has one 

more free parameter than B+BB, the ACstatMi _ m 2 distribution (i.e., CstatMi-CstatM 2 ) resulting from the B+BB 

and C+BB+PL fits to the 10 -5 synthetic spectra has only negative values; this is possible because those two models 

are not nested. Since the two compared models are not nested, the ACstatMi_ m 2 distribution does not respect the 

Wilk’s theorem (i.e., the ACstatMi_ m 2 distribution does not converge to a y 2 distribution with 7-6=1 degree of 

freedom.) Among the 10 5 synthetic spectra, we did not encounter a single case for which ACstatMi _ m 2 is larger 

than or equal to the ACstatMi_ m 2 value resulting from the B+BB and C+BB+PL fits to the real data (i.e., the 

maximum value of the ACstatMi _ m 2 distribution resulting from the simulation is << 508-495=13.) Moreover, the 

ACstatMi _m 2 distribution has a very steep slope so it is extremely unlikely to get a ACstatMi_ m 2 value even close 

to the observed one (i.e., 13). We can, therefore, confidently state that Test 1 is conclusive. 

We now present the results of the log-likelihood statistic test but using C+BB+PL with the parameter values 

resulting from the C+BB+PL fit to the real data as the null-hypothesis (i.e., Test 6.) The resulting ACstatMi_ m 2 

distribution has only positive values and the observed ACstatMi_ m 2 (be., 13) is well within the 1-tr of the distribution 

and very close to its most probable value. We can, therefore, confidently state that Test 6 is inconclusive. 

Based on our previous descriptions of the test results, we can confidently conclude in this case that C+BB+PL is 
a “better description of the data” than B+BB and that in addition C+BB+PL is “statistically significantly better” 
than B+BB. 


2. Comparison of the C+BB and C+BB+PL fits to the GBM data alone in the fine-time spectra of GRB 080916C: 


Here we describe the overall results, but the significance depends on the time interval. However, we always observed 
the same trend, which increases the significance of the overall results compared to the significance obtained for each 
time-interval individually. 

C+BB and C+BB+PL are nested-like models and C+BB+PL has two more free parameters than C+BB; in the 
following, we will consider that C+BB and C+BB+PL correspond, respectively, to the Mi and M 2 models from the 
description of the various tests presented earlier in this Appendix. 

In the fine time intervals, the Cstat values resulting from the C+BB and C+BB+PL fits to the real data are 
usually very similar and Test 1 is typically inconclusive. Similarly, if the true underlying model is C+BB+PL with the 

parameters resulting from the C+BB+PL fit to the real data then we would also expect the value of ACstatMi_ m 2 

computed from the C+BB and C+BB+PL fits to the real data; therefore, Test 6 is also typically inconclusive. 

We can often conclude that (i) if the true underlying model was C+BB with the parameter values resulting from 
the C+BB fit to the real data, then the measured Cstat values obtained when fitting C+BB (i.e., Test 4—goodness of 
fit) and C+BB+PL (i.e., Test 5) to the real data are much higher than expected, and (ii) if the true underlying model 
was C+BB+PL with the parameter values resulting from the C+BB+PL fit to the real data, then the measured Cstat 
values obtained when fitting C+BB (i.e., Test 9) and C+BB+PL (i.e., Test 10—goodness of fit) to the real data are 
typical to what we would expect (i.e., within the 1-sigma confidence region of the Cstat distributions.) We can then 
conclude that Tests 4 and 5 are conclusive while Tests 9 and 10 are inconclusive. 

Finally, we also usually observe that (i) if the true underlying model is C+BB, it is likely that we obtain the measured 
parameters resulting from the C+BB fit to the real data (i.e., Test 3 inconclusive) but it is unlikely that we obtain the 
measured parameters resulting from the C+BB+PL fit to the real data (i.e., Test 2 conclusive.) Conversely, if the true 
underlying model is C+BB+PL, it is likely that we obtain the measured parameters resulting from the C+BB+PL 
and C+BB fits to the real data (i.e., Tests 7 and 8 inconclusive, respectively.) 

Therefore, Tests 1 and 6 are inconclusive, but Tests 3, 7, 8, 9 and 10 are inconclusive and Tests 4 and 5 are conclusive 
so we can conclude that C+BB+PL is overall a ’’better description of the data.” 

_ Let us consider the specific case of the time interval lasting from To+8.5 s to To+9.5 s in GRB 080916C (see Table 

A5 of AppendixO) The C+BB and C+BB+PL fits to the real data result in the same Cstat values (i.e., 526 for 472 
and 470 dof for (J+BB and C+BB+PL, respectively.) The result of Tests 1 and 6 are clearly inconclusive because 
C+BB and C+BB+PL are nested models. It is, therefore, not possible to conclude from Tests 1 and 6 solely that 
C+BB+PL is a better description of the data than C+BB. However, Test 4 indicates that if the true underlying model 
was C+BB, it would be extremely unlikely (i.e., <10” 5 ) to get a Cstat value beyond 526 when fitting C+BB to the 
data; indeed, the Cstat distribution peaks narrowly around 485 when fitting C+BB to the synthetic spectra. The 
same conclusion is drawn for Test 5. In addition, if C+BB was the true underlying model, it would be very likely that 
we measure the spectral parameter values resulting from the C+BB fit to the real data (i.e., Test 3 inconclusive), but 
it would be extremely unlikely that we measure the spectral parameter values resulting from the C+BB+PL fit to the 
real data (i.e., Test 2 conclusive.) This is particularly true for the value of a of C and the value of the amplitude of the 
additional PL of the C+BB+PL model; indeed, the values obtained when fitting C+BB+PL to the real data are ~-0.60 
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and >>0 for a and the amplitude respectively, the values of a and of the amplitude resulting from the C+BB+PL fits 
to the C+BB synthetic spectra peak ~-1.4 and ~0, respectively. None of the 10 5 synthetic spectra based on the C+BB 
null-hypothesis result in parameters values compatible with the observed ones when fitting to C+BB+PL. Conversely, 
if the true underlying model was C+BB+PL, it would be very likely to obtain the observed Cstat values resulting 
from the C+BB and C+BB+PL fits to the real data (i.e., 526—Tests 6, 9 and 10 are inconclusive.) It would also be 
very likely to measured the spectral parameters resulting from the C+BB and C+BB+PL fits to the real data and 
this is especially true for a and the amplitude of the additional PL (i.e., Tests 7 and 8 are conclusive). We can then 
conclude that C+BB+PL is a ’’better description of the data than C+BB.” 


C. SPECTRAL ANALYSIS RESULTS 

The values of the parameters resulting fro m the t ime- integrated spectral analysis of GRBs 080916C and 090926A as 
presented in Section [4] are reported in Tables |A1| & |A2[ respectively. 

The values of the parameters res ultin g from the coarse-time analysis of GRBs 080916C and 090926A as presented 
in Section [5] are reported in Tables |A3| & |A4[ respectively. 

The values of the parameters resu lti ng f rom the fine-time analysis of GRBs 080916C and 090926A as presented in 
Section [6] are reported in Tables X5 & |A6[ respectively. 


D. TIME INTERVAL DEFINITION FOR FINE TIME-RESOLVED SPECTROSCOPY 

Although several techniques have been proposed in the past to automate the definition of time intervals, none seems 
to be really adequate to perform the study that we propose here: tracking the spectral evolution of multiple components 
on the shortest possible time scale—>0.1 s—without including light curve regions with strong spectral evolution in a 
single time-interval as well as in avoiding any possible artificial correlation between energy flux and spectral hardness. 

The simplest automated technique consists of using time intervals of constant duration. If the time intervals are 
chosen to be too large, they may include emission periods with very different spectral properties, making the spectral 
fits unusable for good interpretation. Conversely, if the defined time intervals are too short, we may not have enough 
counts in some time intervals to reconstruct adequately the spectra. 

Another technique consists of defining the time intervals based on the signal to noise (S/N) ratio. The signal is 
integrated until it reaches a certain threshold over the background. When this threshold is reached, another time 
interval is created and so forth. Similarly to the previous technique, this one may also create time intervals including 
emission periods with very different spectral properties. In addition, the S/N ratio is a quantity that should depend 
on the spectral shape. Indeed, the number of counts required to reconstruct a soft spectrum may be very different 
from the number of counts required to reconstruct a hard spectrum and the spectral distribution of the counts is also 
very impo rtant . This is not taken into account in the S/N ratio technique. 

[Scargle (19981 proposed to define time intervals as Bayesian blocks. A new time interval is created when the signal, 
which is integrated over a certain period of time and in a specific energy range, varies significantly based on a criteria 
provided by the user. This technique has two major issues. The first one concerns the energy band that we are 
considering to define significant variations of the signal. Indeed, if we choose an energy band that only includes low 
energy data, it will not be possible to detect signal variations due to a component that will only be intense at high 
energy. Conversely, if we select a high-energy band, we will not be sensitive to variations in the low-energy regime. 
Finally, if we select an energy domain that includes both low- and high-energy counts, any variation of the signal at 
high energy will be hidden by the low-energy variations since the number of counts is usually much larger at low energy 
than at high energy (i.e., non-tliermal spectral shape). The second issue is intrinsic to the technique itself. With the 
Bayesian blocks, we define time intervals in which the intensity of the signal does not change significantly; however, 
the spectral distribution may change a lot. Since one of our studies consists of investigating the correlation between 
the instantaneous flux and the instantaneous hardness (i.e., E pea i;) within each burst, then this technique appears 
to be inappropriate. Indeed, it is not possible to know if two contiguous intervals with similar fluxes have different 
hardness since contiguous time intervals are devised to have significantly different fluxes.Therefore, the Bayesian block 
technique may pro duce a rtificial correlations. 


In Guiriec et al. (20131, we applied an original technique to define the time intervals. To devise intervals with the 


least spectral evolution, we applied a Bayesian block technique not to the count light curve, but to the hardness ratio 
evolution computed from light curves in two different energy ranges. This technique appears to be appropriate for the 
study of GRB 120323A since only two spectral components were identified and that they were not intense in the same 
energy domains. Therefore, we computed the time-history hardness ratio in using two energy bands in which each 
component was the most intense. Such a technique is difficult to apply for the current study since we are considering 
three spectral components whose energy domains in which they are intense overlap. 

All the techniques presented here use criteria defined by the user such as energy bands and thresholds. In addition, 
they are applied to observed count light curves, which are tainted with energy dispersion effects. Indeed, a photon 
has only a certain probability to be measured as a count at its true energy and the measured energy can be more or 
less spread around the true photon energy. While the spectral analysis corrects for this effect, a direct analysis of the 
count light curve may be significantly biased. Finally, the user usually defines the energy bands in the observer frame 
while GRBs are known to have large spread in distances. Defining the energy bands in the rest frame is not always 
possible since very few GRBs have redshift measurements. In addition, the energy in which each spectral component 
is the most intense does not only depend on the redshift but also on the micro-physics. For instance, the energy of 
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the emitted radiation depends on the magnetic field profile and intensity as well as on the velocity profile of the solid 
layers within the jet. 

Our procedure to define the fine-time intervals of our analysis is described in the following and it is done using 
Rrnfit. We wanted to define time intervals as short as possible—0.1 s being the shortest allowed time duration -to 
avoid strong spectral evolution within a time interval and to have as many time intervals as possible to track the 
spectral evolution in the both the rising and decaying parts of the pulses. To do so, we combined the base 0.1 s time 
bins until 


• we could adequately fit a Band function to the data (i.e., convergence and amplitude at least 2-cr above 0), 

• the parameter values of the amplitude, a and E pea k resulting from the Band fit to the real data were within 
the 1—<7 confidence region of the distributions of the amplitude, a and E pea i; values, respectively, resulting from 
Band fit to the Band synthetic spectra (i.e., see Appendix |B| 

• the biases on the parameter values for the amplitude, a and E pea it when fitting a Band function alone to the data 
was <l-cr (i.e., the distribution of 1 1 | xUn | ! —-—where X; best is the best value of the parameter resulting from the 
Band fit to the real data, x+ ax is the value of x, at maximum of the distribution resulting from the Band fit to 
the Band synthetic spectra and x“ n is the 1-cr uncertainty on x best —is compatible with 0 within 1-cr), 

• the 1-cr widths of the distributions of the amplitude, a and E pea k were smaller than the 1-cr uncertainties of the 
amplitude, a and E pea k parameters, respectively, resulting from the Band function fit to the real data. 

• the best parameter values for amplitude, a and E pea k did not vary by more than 10% when adding a 0.1 s time 
bin either toward the positive or negative values of time, or by sliding the time interval window by 0.1 s either 
toward the positive or negative values of time. 


E. FINE TIME-RESOLVED ANALYSIS 

Figure [AT| shows the results of the fine time-resolved analysis of GRB 080916C with a Band function alone and the 
C+BB+ PL m odel as presented in Section [6] 

Figure [A~2] shows the results of the fine time-resolved analysis of GRB 090926A with a Band function alone and the 
C+BB+ FL m odel as presented in Section [6] 

Figure A3 shows the results of the fine time-resolved analysis of GRB 080916C with the C+BB+PL 5 pQr . ams model 
as prese nted in Section [8| 

Figure A4 shows the results of the fine time-resolved analysis of GRB 090926A with the C+BB+PL 5 pQI . ams model 
as presented in Section [8] 
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Energy (keV) Energy (keV) 

Fig. Al.— GRB 080916C : i/F u spectra resulting from the fine-time analysis. The solid yellow and black lines correspond to the best 
Band-only and C+BB+PL fits, respectively. The dashed yellow and black lines correspond to the 1-cr confidence regions of the Band-only 
and C+BB+PL fits, respectively. The solid blue, red and green lines correspond to the cutoff power law, to the BB component and to the 
additional power law resulting from the best fit with the C+BB+PL model (i.e., solid black line) to the data, respectively. 
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Energy (keV) Energy (keV) Energy (keV) 

Fig. A2.— GRB 090926A : spectra resulting from the fine-time analysis. The solid yellow and black lines correspond to the best 

Band-only and C+BB+C2 fits, respectively. The dashed yellow and black lines correspond to the 1 -a confidence region of the Band-only 
and C+BB+C2 fits, respectively. The solid blue, red and green lines correspond to the cutoff power law, to the BB component and to the 
additional cutoff power law resulting from the best fit with the C+BB+C2 model (i.e., solid black line) to the data, respectively. 
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Fig. A3. — GRB 080916C : vF u spectra resulting from the fine-time analysis using the C+BB+PL 5 params model (i.e., a=-0.7 and 
index=-1.5). The solid black lines correspond to the best C+BB+PL fits and the dashed ones to the 1-cr confidence regions. The solid 
blue, red and green lines correspond to the cutoff power laws, to the BB components and to the additional power laws resulting from the 
best fits with the C+BB+PL model. 
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Fig. A4.— GRB 090926A : uF u spectra resulting from the fine-time analysis using the C+BB+PL 5 params model (i.e., a=-0.7 and 
index=-1.5). The solid black lines correspond to the best C+BB+PL fits and the dashed ones to the 1 -cr confidence regions. The solid 
blue, red and green lines correspond to the cutoff power laws, to the BB components and to the additional power laws resulting from the 
best fits with the C+BB+PL model. 
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TABLE A1 

Parameters of the various tested models resulting from the time-integrated spectral analysis of GRB 080916C from 

Tq-0.1 S to Tq + 71 S WITH THEIR ASYMMETRICAL 1 -cr UNCERTAINTIES. 


Models 

Sta 

ndard Me 

idel 

Additional Model 

EAC 

Cstat/dof 


C or Band 

BB 

PL, 

C or 

Band 



Parameters 

^peak 

a 

& 

kT 

a 

/3 

E 0 

nO/bO 

n3/b0 

n4/b0 

L/bO 

F/bO 

B/bO 



(keV) 



(keV) 



(MeV) 














b0+n0+n3+n4 







Band 

451 

-1.03 

-2.18 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

666/473 


±25 

±0.02 

±0.08 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 


B + BB 

1321 

-1.30 

-2.41 

45.0 

- 

- 

- 

- 

- 

- 

- 

- 

- 

592/471 


+366 

-243 

±0.03 

+0.17 

-0.30 

±1.6 

- 

- 

- 

- 

- 

- 

- 

- 

- 


B + PL 

269 

-0.35 

-2.03 


-1.84 

- 

- 

- 

- 

- 

- 

- 

- 

612/471 


+ 24 
-20 

±0.16 

±0.06 


+0.08 

-0.13 

- 

- 

- 

- 

- 

- 

- 

- 


B+BB+PL 

866 

-1.06 

-2.47 

11.3 

-1.77 

- 

- 

- 

- 

- 

- 

- 

- 

587/469 


+261 

-178 

±0.15 

+0.25 

-0.98 

±2.3 

+0.10 

-0.51 

- 

- 

- 

- 

- 

- 

- 

- 


C + BB 

1858 

-1.33 

- 

44.9 


- 

- 

- 

- 

- 

- 

- 

- 

597/472 


+330 

-295 

±0.02 

- 

±1.4 


- 

- 

- 

- 

- 

- 

- 

- 


C + PL 

362 

-0.63 

- 


-1.63 

- 

- 

- 

- 

- 

- 

- 

- 

650/472 


+23 

-21 

±0.09 

- 


±0.02 

- 

- 

- 

- 

- 

- 

- 

- 


C+BB+PL 

912 

-1.04 

- 

11.7 

-1.67 

- 

- 

- 

- 

- 

- 

- 

- 

589/470 


+307 

-179 

±0.15 

- 

±2.4 

±0.05 

- 

- 

- 

- 

- 

- 

- 

- 








b0+n0+n3+n4 

+ EAC 







Band 

458 

-1.03 

-2.18 

- 

- 

- 

- 

0.99 

0.98 

0.98 

- 

- 

- 

665/470 


±27 

±0.02 

±0.08 

- 

- 

- 

- 

±0.03 

±0.03 

±0.03 

- 

- 

- 


B+BB 

1311 

-1.31 

-2.43 

44.8 

- 

- 

- 

1.07 

1.05 

1.06 

- 

- 

- 

588/468 


+357 

-238 

±0.03 

+0.17 

-0.32 

±1.5 

- 

- 

- 

±0.04 

±0.03 

±0.03 

- 

- 

- 


B + PL 

259 

-0.30 

-2.03 


-1.84 

- 

- 

1.05 

1.04 

1.04 

- 

- 

- 

609/468 


+23 

-20 

±0.17 

±0.06 


+0.08 

-0.11 

- 

- 

±0.03 

±0.03 

±0.03 

- 

- 

- 


B+BB+PL 

898 

-1.09 

-2.47 

11.7 

-1.78 

- 

- 

1.06 

1.05 

1.05 

- 

- 

- 

584/466 


+290 

-190 

±0.15 

+0.23 

-1.13 

±2.3 

+0.11 

-0.65 

- 

- 

±0.04 

±0.04 

±0.04 

- 

- 

- 


C + BB 

1817 

-1.34 

- 

44.7 

- 

- 

- 

1.07 

1.06 

1.06 

- 

- 

- 

594/469 


+319 

-289 

±0.02 

- 

±1.3 

- 

- 

- 

±0.03 

±0.03 

±0.03 

- 

- 

- 


C + PL 

366 

-0.64 

- 

- 

-1.63 

- 

- 

1.00 

0.98 

0.99 

- 

- 

- 

649/472 


±24 

±0.09 

- 

- 

±0.02 

- 

- 

±0.03 

±0.03 

±0.03 

- 

- 

- 


C+BB+PL 

945 

-1.08 

- 

42.0 

-1.67 

- 

- 

1.06 

1.05 

1.05 

- 

- 

- 

586/467 


+323 

-194 

±0.15 

- 

±2.3 

±0.05 

- 

- 

±0.04 

±0.03 

±0.03 

- 

- 

- 







bO+n 

0+n3+n4 + LLE+Front+Back 






Band 

459 

-1.04 

-2.19 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

728/541 


+23 

-21 

±0.01 

±0.01 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 


B+BB 

1289 

-1.30 

-2.27 

45.0 

- 

- 

- 

- 

- 

- 

- 

- 

- 

650/539 


+ 228 
-194 

±0.02 

±0.03 

±1.6 

- 

- 

- 

- 

- 

- 

- 

- 

- 


B + PL 

333 

-0.68 

-2.17 

- 

-2.12 

- 

- 

- 

- 

- 

- 

- 

- 

685/539 


±20 

±0.08 

±0.01 

- 

+ 0.11 
-0.24 

- 

- 

- 

- 

- 

- 

- 

- 


B+BB+PL 

1026 

-1.18 

-2.25 

42.2 

-2.04 

- 


- 

- 

- 

- 

- 

- 

646/537 


+280 

-187 

±0.08 

+0.03 

-0.06 

±2.2 

+ 0.14 
-0.85 

- 


- 

- 

- 

- 

- 

- 


B+BB+C 

1042 

-1.18 

-2.27 

42.2 

-1.97 

- 

26640 

- 

- 

- 


- 

- 

646/536 


+308 

-200 

±0.09 

±0.13 

±2.2 

+0.11 

-0.82 

- 

+ 28200 
-6450 

- 

- 

- 


- 

- 







bO+nO+n 

3 + n4+LLE+Fr< 

ant + Back + EAC 






Band 

452 

-1.03 

-2.19 


- 

- 

- 

- 

- 

- 

0.92 

1.16 

1.16 

722/538 


±23 

±0.01 

±0.05 


- 

- 

- 

- 

- 

- 

+0.23 

-0.20 

+0.49 

-0.36 

+0.46 

-0.35 


B + BB 

1188 

-1.29 

-2.24 

15.2 

- 

- 

- 

- 

- 

- 

1.02 

1.18 

1.18 

649/536 


+ 241 
-184 

±0.02 

±0.06 

±1.6 

- 

- 

- 

- 

- 

- 

+0.23 

-0.19 

+0.47 

-0.35 

+0.45 

-0.34 


B + PL 

302 

-0.55 

-2.08 


-1.99 

- 

- 

- 

- 

- 

1.34 

2.17 

2.13 

673/536 


±21 

±0.11 

±0.05 


+0.09 

-0.16 

- 

- 

- 

- 

- 

+0.30 

-0.26 

+0.46 

-0.62 

+0.78 

-0.46 


B+BB+PL 

914 

-1.12 

-2.34 

11. G 

-1.91 

- 


- 

- 

- 

0.77 

0.92 

0.91 

643/534 


+ 245 

±0.10 

+0.14 

±2.3 

+0.05 






+0.30 

+0.52 

+0.50 



-184 

-0.17 

-0.13 






-0.26 

-0.37 

-0.30 


B+BB+C 

975 

-1.13 

-2.45 

41.8 

-1.87 

- 

12100.0 

- 

- 

- 

0.67 

0.85 

0.84 

641/533 


+286 

±0.10 

+0.16 

±2.3 

±0.06 


+ 28200 




+0.29 

+0.53 

+0.50 



-200 

-0.22 


-6450 




-0.24 

-0.36 

-0.35 
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TABLE A2 

Parameters of the various tested models resulting from the time-integrated spectral analysis of GRB 090926A from To s 

TO Tq+20 s with their asymmetrical 1 -a uncertainties. 


Models 

Sta 

ndard Model 


Addit 


anal Modi 

31 

EAC 

Cstat/dof 


C or Band 

BB 

PL, C or Band 



Parameters 

■^peak 

a 

£ 

kT 

a 


13 

E 0 

bl/bO 

n3/b0 

n6/b0 

n7/b0 

L/bO 

F/bO 

B/bO 



(keV) 



(keV) 




(MeV) 













bO+bl+n 

3+n6+n7 

with b0>756.48 keV a 

nd bl>609.51 keV 






Band 

294 

-0.75 

-2.44 

- 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 

1138/566 


±5 

±0.01 

±0.03 

- 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 


B + BB 

450 

-1.01 

-2.58 

39.4 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 

999/564 


±16 

±0.02 

±0.06 

±0.8 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 


B + PL 

259 

-0.51 

-2.37 


-2.02 


- 

- 

- 

- 

- 

- 

- 

- 

- 

1057/564 


±5 

±0.04 

±0.03 


+0.13 

-0.24 



- 

- 

- 

- 

- 

- 

- 

- 


B+BB+PL 

436 

-0.96 

-3.37 

37.3 

-1.62 


- 

- 

- 

- 

- 

- 

- 

- 

- 

974/562 


±15 

±0.05 

+0.29 

-0.68 

±1.1 

+0.06 

-0.04 



- 

- 

- 

- 

- 

- 

- 

- 


C+BB+PL 

441 

-0.90 

- 

37.2 

-1.61 


- 


- 

- 

- 

- 

- 

- 

- 

976/563 


±13 

±0.05 

- 

±1.0 

+0.05 

-0.03 


~ 


- 

- 

- 

- 

- 

- 

- 


C+BB+C 

425 

-0.90 

- 

37.6 

-1.51 


- 

21.050 

- 

- 

- 

- 

- 

- 

- 

967/562 


±14 

±0.06 

- 

±1.2 

+0.09 

-0.05 



+ 12.800 
-6.440 

- 

- 

- 

- 

- 

- 

- 





b0+bl+n3+n 

i6 + n7 with b0>756.48 keV and bl>609.51 keV + EAC 





Band 

294 

-0.75 

-2.44 

- 

- 


- 

- 

1.02 

1.03 

0.99 

0.99 

- 

- 

- 

1078/562 


±5 

±0.01 

±0.03 

- 

- 


- 

- 

±0.07 

±0.07 

±0.07 

±0.07 

- 

- 

- 


B + BB 

475 

-1.02 

-2.52 

38.9 

- 


- 

- 

1.05 

0.96 

0.92 

0.92 

- 

- 

- 

937/560 


±22 

±0.02 

±0.07 

±0.7 

- 


- 

- 

±0.07 

±0.07 

±0.07 

±0.07 

- 

- 

- 


B + PL 

257 

-0.51 

-2.33 

- 

-2.00 


- 

- 

0.99 

0.96 

0.92 

0.92 

- 

- 

- 

1000/560 


±6 

±0.04 

±0.04 

- 

+0.11 

-0.24 


~ 

- 

±0.07 

±0.07 

±0.07 

±0.07 

- 

- 

- 


B+BB+PL 

443 

-0.92 

-3.40 

37.3 

-1.60 


- 

- 

1.10 

1.08 

1.04 

1.03 

- 

- 

- 

916/558 


±20 

±0.05 

+0.33 

-0.75 

±1.1 

+0.07 

-0.05 


~ 

- 

±0.07 

±0.08 

±0.08 

±0.08 

- 

- 

- 


C+BB+PL 

446 

-0.91 

- 

37.1 

-1.60 


- 

- 

1.11 

1.09 

1.05 

1.04 

- 

- 

- 

918/559 


±18 

±0.04 

- 

±1.0 

+0.05 

-0.04 


~ 

- 

±0.07 

±0.08 

±0.08 

±0.08 

- 

- 

- 


C+BB+C 

432 

-0.91 

- 

37.6 

-1.50 


- 

21.130 

1.10 

1.07 

1.03 

1.03 

- 

- 

- 

909/558 


±20 

±0.06 

- 

±1.2 

+0.15 

-0.06 


~ 

+ 12.900 
-6.560 

±0.07 

±0.08 

±0.08 

±0.08 

- 

- 

- 





b0+bl+n3+r 

i6+n7 w: 

ith b0>756.48 keV and 

bl >609.51 keV + 

LLE + Front + Back 




Band 

280 

-0.73 

-2.29 

- 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 

1234/633 


±4 

±0.01 

±0.01 

- 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 


B + BB 

405 

-0.99 

-2.28 

40.9 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 

1100/631 


±18 

±0.02 

±0.01 

±0.8 

- 


- 

- 

- 

- 

- 

- 

- 

- 

- 


B + PL 

251 

-0.46 

-2.31 

- 

-1.95 


- 

- 

- 

- 

- 

- 

- 

- 

- 

1121/631 


±4 

±0.03 

±0.02 

- 

+0.03 

-0.05 



- 

- 

- 

- 

- 

- 

- 

- 


B+C 

254 

-0.46 

-2.38 

- 

-1.85 


- 

1319 

- 

- 

- 

- 

- 

- 

- 

1110/630 


±5 

±0.03 

±0.03 

- 

±0.02 


~ 

+948 

-503 

- 

- 

- 

- 

- 

- 

- 


B+BB+PL 

370 

-0.84 

-2.34 

38.0 

-1.88 


- 

- 

- 

- 

- 

- 

- 

- 

- 

1091/629 


±23 

±0.06 

±0.03 

±1.5 

±0.04 


- 

- 

- 

- 

- 

- 

- 

- 

- 


B+BB+C 

415 

-0.87 

-2.68 

36.9 

-1.73 


- 

789.600 

- 

- 

- 

- 

- 

- 

- 

1056/628 


±18 

±0.04 

+0.10 

-0.13 

±1.1 

±0.03 



+207.000 

-145.000 

- 

- 

- 

- 

- 

- 

- 


B+BB+B2 

430 

-0.90 

<-4 

37.6 

-1.53 


-2.29 

33.140 

- 

- 

- 

- 

- 

- 

- 

1016/627 


±14 

±0.05 

- 

±1.1 

+0.09 

-0.06 


±0.05 

+22.400 

-9.600 

- 

- 

- 

- 

- 

- 

- 


C+BB+C 

451 

-0.86 

- 

36.0 

-1.73 


- 

617.7 

- 

- 

- 

- 

- 

- 

- 

1084/629 


±11 

±0.03 

- 

±0.9 

±0.02 


- 

+ 104.0 
-82.9 

- 

- 

- 

- 

- 

- 

- 


C+BB+B2 

430 

-0.90 

- 

37.6 

-1.53 


-2.29 

33.770 

- 

- 

- 

- 

- 



1016/628 


±14 

±0.05 

- 

±1.1 

+0.09 

-0.05 


±0.05 

+ 15.600 
-10.200 

- 

- 

- 

- 

- 






b0+bl + 

n3+n6+r 

i7 with b0>756.48 

keV and bl>609.51 keV + LLE + Front 

+ Back - 

f EAC 




Band 

293 

-0.75 

-2.42 

- 

- 


- 

- 

- 

- 

- 

- 

0.47 

0.39 

0.11 

1189/630 


±4 

±0.01 

±0.03 

- 

- 


- 

- 

- 

- 

- 

- 

±0.06 

±0.07 

±0.07 


B + BB 

439 

-1.00 

-2.49 

39.7 

- 


- 

- 

- 

- 

- 

- 

0.38 

0.28 

0.30 

1056/628 


±17 

±0.02 

±0.04 

±0.8 

- 


- 

- 

- 

- 

- 

- 

±0.07 

±0.08 

±0.08 


B + PL 

259 

-0.52 

-2.36 

- 

-2.05 


- 

- 

- 

- 

- 

- 

0.64 

0.62 

0.64 

1105/628 


±5 

±0.04 

±0.03 

- 

+0.08 

-0.16 


- 

- 

- 

- 

- 

- 

±0.08 

+0.15 

-0.11 

+0.15 

-0.11 


B+C 

258 

-0.49 

-2.39 


-1.90 


- 

2564 

- 

- 

- 

- 

0.74 

0.84 

0.87 

1102/627 


±5 

±0.03 

±0.03 


+0.05 

-0.08 


- 

+4330 

-1140 

- 

- 

- 

- 

+0.14 

-0.11 

+0.29 

-0.21 

+0.29 

-0.21 


B+BB+PL 

425 

-0.93 

-2.55 

37.8 

-1.94 


- 

- 

- 

- 

- 

- 

0.35 

0.31 

0.33 

1049/626 


±19 

±0.04 

±0.06 

±1.1 

±0.05 


- 

- 

- 

- 

- 

- 

±0.07 

±0.10 

±0.10 


B+BB+C 

430 

-0.88 

-3.20 

36.8 

-1.66 


- 

1098 

- 

- 

- 

- 

1.46 

2.86 

2.91 

1029/625 


±13 

±0.04 

+0.27 

-0.47 

±1.0 

±0.04 


- 

+301 

-212 

- 

- 

- 

- 

±0.38 

+0.96 

-0.87 

+0.94 

-0.89 


B+BB+B2 

425 

-0.90 

<-10 

37.6 

-1.51 


-2.27 

20.030 

- 

- 

- 

- 

0.83 

0.88 

0.91 

1014/624 


±15 

±0.06 

- 

±1.2 

±0.07 


±0.08 

+ 19.600 
-8.040 

- 

- 

- 

- 

+0.21 

-0.17 

+0.32 

-0.25 

+0.32 

-0.25 


C+BB+C 

439 

-0.87 

- 

36.7 

-1.63 


- 

1041 

- 

- 

- 

- 

1.92 

3.88 

3.93 

1032/626 


±12 

±0.04 

- 

±1.0 

±0.03 


- 

+273 

-194 

- 

- 

- 

- 

±0.28 

+0.88 

-0.73 

+0.84 

-0.71 


C+BB+B2 

425 

-0.90 

- 

37.6 

-1.52 


-2.27 

20.580 

- 

- 

- 

- 

0.82 

0.86 

0.90 

1014/625 


±15 

±0.05 

- 

±1.2 

+0.09 

-0.05 


±0.08 

+ 17.100 
-8.410 

- 

- 

- 

- 

+0.21 

-0.15 

+0.32 

-0.23 

+0.32 

-0.23 





44 


S.Guiriec 


TABLE A3 

Parameters of the various tested models resulting from the coarse-time spectral analysis of GRB 080916C with their 

ASYMMETRICAL 1-0 UNCERTAINTIES. 


Time from T n 

Models 

Standard Model 

Additional Moi 

d. 

el 

EAC 

Cstat/dof 



C or Band 

BB 

PL . 

d: 

r C 



^start 

^stop 

Parameters 

^peak 

a 

P 

kT 

a 

r 

E 0 

nO/bO 

n3/b0 

n4/b0 

L/bO 


(s) 

(s) 


(keV) 



(keV) 


L 

(MeV) 










2 25+ 0 12 

-2.25 — o 16 

bO-f 

-n0+n3+n4 








- 0.1 

+4.3 

Band 

454+38 

-0.69+0.04 

- 

- 


- 

- 

- 

- 

- 

505/473 



B + BB 


-0.97+0.04 

<-5 

39.8+2.6 

- 


- 

- 

- 

- 

- 

477/471 



C+BB 


-0.97+0.04 

- 

39.8+2.7 

- 


- 

- 

- 

- 

- 

477/472 





-2 26~*" 012 
^ -0.16 

b 0 +n 0 -| 

-n3+n4 + EAC 









Band 

459 + 42 

-0.69+0.04 

- 

- 


- 

1.03+0.06 

0.98+0.05 

0.98+0.05 

- 

502/470 



B + BB 

1183 ±?82 

-0.97+0.04 

<-5 

40.0 + 2.6 

- 


- 

1.08+0.06 

1.03 + 0.05 

1.03+0.05 

- 

573/468 



C+BB 

+ 212 
USO-lgO 

-0.97+0.04 

- 

40.0 + 2.6 

- 


- 

1.08+0.06 

1.03 + 0.05 

1.03+0.05 

- 

573/469 






bO+nO+n3+n4+LLE 










Band 

496+33 

-0.72+0.03 

-2.67+0.08 

- 

- 


- 

- 

- 

- 

- 

524/484 



B + BB 

1069”*” 22 ^ 

iuoy_i 76 

-0.95+0.05 

2 q *+°- 12 
-2.95_q 15 

39.4+2.8 

- 


- 

- 

- 

- 

- 

493/482 



C+BB 

1193+ 217 

1 Ayd —186 

-0.97+0.04 

- 

39.8+2.7 

- 


- 

- 

- 

- 

- 

493/483 





9 90 +O.II 
-2.25_o i4 

bO + nO+n3 + n4+LLE + EAC 




7 8 + 6-6 
7 -°_3 9 




Band 

451+38 

-0.69+0.04 

- 

- 


- 

- 

- 

- 

517/483 



B + BB 

10S1 ±?93 

-0.94+0.05 

2 on + 6- 32 
-2.80 — 0 48 

39.4+2.8 

- 


- 

- 

- 

- 

1 58 + 2 ' 97 
1,5 ° —1.17 

493/481 



C+BB 


-0 97 + 0 - 04 
u ‘ y '— 0.12 


39.8+2.7 

- 


- 

- 

- 

- 

1 kc + 3-00 
i.58_1 10 

493/483 



1291 + 250 


-2 oi - * -0,10 

^ UA —0.13 

bO-f 

-n0+n3+n4 








+4.3 

+7.5 

Band 

-1.15+0.03 

- 

- 


- 

- 

- 

- 

- 

519/473 



B + BB 

2464+0™ 

-1.23+0.03 

-2 10+ 013 

AU —0.18 

44 4 + 6 ' 7 
-5.9 

- 


- 

- 

- 

- 

- 

508/471 



B + PL 

967+ 247 

4,0 -178 

-0.97+0.13 

<-5 

- 

-1 47 + 0 - 10 
1,4 -0.05 


- 

- 

- 

- 

- 

517/471 



C+BB+PL 

1554+331 

1554 — 268 

-0.70+0.10 

- 

39.4+3.8 

-1.56 + 0.04 


- 

- 

- 

- 

- 

495/470 



1299”*” 337 

_266 


-2 ()1 + 010 
^ UA —0.13 

b 0 +n 0 -| 

-n3+n4 + EAC 









Band 

-1.15+0.03 

- 

- 


- 

1.06+0.07 

0.96 + 0.06 

1.02+0.06 

- 

511/470 



B + BB 

oofio+647 

z383_549 

-1.24+0.03 

9 1 o+6-13 
'^ A -0.18 

44 - i ±*:l 

- 


- 

1.13+0.08 

1.04+0.07 

1.10+0.07 

- 

500/468 



B + PL 

»90±?SS 

-0.97+0.13 

<-5 

- 

+ 0.10 
14 -0.05 


- 

1.05+0.07 

0.95+0.06 

1.01+0.06 

- 

510/468 



C+BB+PL 

1520+2og 

-0.70+0.10 

- 

39.5 + 3.5 

-1.58 + 0.04 


- 

1.14+0.08 

1.04+0.07 

1.10+0.07 

- 

487/467 



1400”*” 3 '*' 2 

i4tuu_243 



bO+nO+n3+n4+LLE 










Band 

-1.16+0.03 

-2.08+0.04 

- 

- 


- 

- 

- 

- 

- 

546/484 



B + BB 

2474+ 790 

-459 

-1.23+0.03 

9 1O+0-04 
" Z A -0.06 

44 1 + 6 ' 7 
**-5.7 

- 


- 

- 

- 

- 

- 

534/482 



B+C 

1092 + 5|« 

-1.03+0.15 

<-2.5 

- 

-1.48 + 0.05 


104.3+3|-8 

- 

- 

- 

- 

538/481 



B+BB+C 

looolis? 

-0.70+0.10 

<-2.5 

39.3+3.6 

! 58 + 0.02 
A ' 58 —0.06 


i 2 6.8±“i; 

- 

- 

- 

- 

515/479 



C+C2 

10«l?il 

-0.99+0.13 

- 

- 

1 48 + 0.12 
-l-48_ 0 05 


98.5 4 if 7 ;i 

- 

- 

- 

- 

539/482 



C+BB+PL 

1786+|“| 

-0.70+0.10 

- 

39.1 + 3.3 

-1.67+0.01 


- 

- 

- 

- 

- 

543/481 



C+BB+C 

I 6 O 8+275 

-0.70+0.10 

- 

39.3+3.7 

-1.58 + 0.02 


1 17.3^ 3 x g 

- 

- 

- 

- 

516/480 



1418+247 



bO+nO+n3 + n4+LLE + EAC 




n 88+ 9 ■ 34 
°- 88 — 0.26 




Band 

-1.16+0.03 

-2.10+0.10 

- 

- 


- 

- 

- 

- 

545/483 



B + BB 

2615 + gQ 8 

-1.24+0.03 

-2.17+0.11 

—±I: 4 9 

- 


- 

- 

- 

- 

0 91+6-83 

uyi — 0.25 

534/481 



B + C 

lOH—228 

-1 03 - *" 0-18 

A U -0.16 

<-5 


-1.42 + 0.05 


101 o - * -26-9 
iui.U -20 o 

- 

- 

- 

1.68+0.38 

535/480 



B+BB+C 

i K«n+334 
1569L 2 g5 

-0.70+0.10 

<-3 

39 5iJ;° 

-1 54 + 0 - 04 
A ‘°* — 0.06 


191 q+69-7 
121 -8_ 22.9 

- 

- 

- 

1 34 + 0 ' 44 

A t3 -0.42 

515/478 



c+c 

981+200 

-1.00+0.17 

- 

- 

-1.43 + 0.05 


102.6 + 26.4 

- 

- 

- 

1 69 - * -0 ' 37 
iD — 0.31 

535/481 



C+BB+PL 

1739+390 

i/5y_3H 

-0.70+0.10 

- 

39.3+3.5 

-1 64 +0 06 
1-0 -0.14 


- 

- 

- 

- 

1 97+0.79 

A -0.89 

543/480 



C+BB+C 

1526+335 

1525 — 268 

-0.70+0.10 


39.5+3.8 

-1.54 + 0.04 


120.0tl 3 2 4 4 

- 


- 

1 41 + 041 

1 ' 4i — 0.36 

515/479 
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Time from Tg 

Models 

Standard Model 


Additional Me 

)del 

EAC 

Cstat/dof 



C or Band 

BB 

PL 

or C 



Tstart 

Tstop 

Parameters 

■^peak 

a 

P 

kT 

a 

E 0 

nO/bO 

n3/b0 

n4/b0 

L/bO 


(s) 

(s) 


(keV) 



(keV) 


(MeV) 










-2 29"*" 0,11 
z -0.17 

b0+n0+n3+n4 







+7.5 

+71.0 

Band 

413 + 26 

-1.07+0.02 

- 

- 

- 

- 

- 

- 

- 

618/473 



B + BB 

4,0 -210 

-1.34+0.04 

9 OO+0.17 
-^. 4 Z _0 29 

45.5 + 2.0 

- 

- 

- 

- 

- 

- 

569/471 



B + PL 

263+25 

-0.40 + 0.19 

-2 16 + 009 

Z i -0.14 

- 

-1 84 +0-09 
10 -0.16 

- 

- 

- 

- 

- 

579/471 



C+BB 

l«3tSI 

-1.39+0.03 

- 

44.8+1.7 

- 

- 

- 

- 

- 

- 

573/472 



c+c 


-1.42+0.02 

- 

- 

+0.95+0.22 

64 + 7 keV 

- 

- 

- 

- 

573/471 



C+BB+PL 

681+f™ 

—155 

-1 io +°- 20 
±1U —0.27 

- 

42.9 + 3.4 

-1 69"*” 9-2 3 
1-0 -0.07 

- 

- 

- 

- 

- 

570/470 



C+BB+C 


-0.72 + 0.42 

- 

4 i.5+J;S 

-1.53+0.08 

7 247 + 11-600 
'•^ 4 '_3.220 

- 

- 

- 

- 

567/469 





9 ofi + 0.11 

-4.48_o 17 

b0+n0+n3+n4 + EAC 








Band 

-1.07+0.02 

- 

- 

- 

0.98 + 0.04 

0.98+0.04 

0.98+0.04 

- 

618/470 



B + BB 

96 + 330 

-1.35+0.04 

-2 34 + 0-18 
—0.32 

45.2+1.9 

- 

- 

1.06 + 0, 05 

1.06+0.05 

1.06+0.05 

- 

568/468 



B + PL 

255+25 

-0 36+ 0 12 
UcSO -0.18 

-2 16 + 0- 09 
Zit> -0.14 

- 

1 84 + 0-09 
-1 ‘ 84 —0.14 

- 

1.04+0.04 

1.04+0.05 

1.04+0.04 

- 

578/468 



C+BB 

1595+ 499 

±oyo_44 8 

-1.40+0.04 

- 

44.6+1.6 

- 

- 

1.06+0.05 

1.07+0.05 

1.06+0.05 

- 

571/469 



C+C 

1773+ 389 

1 ' ' J — 329 

-1.42+0.01 

- 

- 

+ 1.17+0.23 

57+6 keV 

1.05+0.05 

1.06+0.05 

1.05+0.05 

- 

571/468 



C+BB+PL 

714+ 473 

' 14 — 174 

1 15+ 0 - 22 
A * 10 — 0.16 

- 

43.3+3.1 

-1 69+ 0 - 35 
loy —0.07 

- 

1.05+0.05 

1.06+0.05 

1.05+0.05 

- 

569/467 



C+BB+C 

472 + 242 
z —148 

n cn +°- 81 
-0.80_o.4i 

- 

42.1+3;® 

-1 54 + 0 " 21 
±- ° 4 —0.09 

7343 ™° 

1.05+0.05 

1.06+0.05 

1.05+0.05 

- 

566/466 






bO+r 

iO+n3+n4+LLE 








Band 

402 + 23 

-1.07+0.02 

-2.17+0.01 

- 

- 

- 

- 

- 

- 

- 

632/484 



B + BB 

866+j33 

-1.33+0.04 

-2.19+0.03 

46.1 + 2.1 

- 

- 

- 

- 

- 

- 

583/482 



B + PL 

283+3| 

-0.54+0.11 

9 on + 0.05 

-4.4U_o 08 

- 

-1 93+0-0® 

± y -0.15 

- 

- 

- 

- 

- 

592/482 



B+BB+PL 

+294 

003 — 147 

+0.16 

±i4 -0.13 

9 OQ + 0.11 

-4.48_o 13 

43.1 + 3.2 

, o 4 +0.14 
-i.84_o sg 

- 

- 

- 

- 

- 

581/480 



B+BB+C 

691+308 

Dy 1 —180 

-1.14+0.19 

-2.56+0- 4 ® 

43.1+3.5 

-1.74 + 0.02 

362 «+ 888 -0 
804.o_ 1210 

- 

- 

- 

- 

580/479 



C+C 

336+20 

-0.66+0.09 

- 

- 

-1.75+0.02 

q + 121.0 

247 - 8 _ 62 .5 

- 

- 

- 

- 

608/482 



C+BB+PL 

A/1R +193 

6 46 _h 2 

-1.02+0.13 

- 

41.0 + 2.8 

-1.79+0.02 

- 

- 

- 

- 

- 

593/481 



C+BB+C 

7 3S+535 
73 ° — 176 

-1.15+0.20 

- 

43.2 + 3.2 

-1.71+0.05 

237.9+g 21 {° 

- 

- 

- 

- 

581/480 






bO+nO+i 

i3+n4+LLE + 

EAC 




0 96+ 9-3 ® 
uyo —0.27 




Band 

403+24 

-1.07+0.02 

-2.18+0.07 

- 

- 

- 

- 

- 

- 

632/483 



B + BB 

070+246 
873 _176 

-1.33+0.04 

-2.15+0.08 

46.0 + 2.1 

- 

- 

- 

- 

- 

1 17 +0.42 
— 0.33 

583/481 



B + PL 

271 + 22 

0 48 + 0-09 

-0.48_o.l4 

-2 14+0- 08 
- 0.10 

- 

-1 9o + 0 - 07 
Ayu —0.14 

- 

- 

- 

- 

o 7 +0.48 

-0.38 

591/481 



B+BB+PL 

002 + 3^2 

-1.14+0.15 

9 oc+0.17 

-4.48_o 27 

43.1+3.3 

-1 84+ 0 - 15 
1CS4 —0.48 

- 

- 

- 

- 

0 99 + 0 - 47 

uyy — 0.42 

581/479 



B+BB+C 

077 ±”g 

- 1.12 + 0.22 

9 eo + 0.38 
-4°8 _ 1 83 

43.1 + 3.5 

-1.73+0.02 

320 7+ 223 0 

84U. t _ 94 7 

- 

- 

- 

1 15+ 0 - 74 

—0.51 

580/478 



C+C 

306+22 

-0.52 + 0.12 

- 

- 

-1.67+0.03 

214 6+ 815 
Z1 ^' D —47.9 

- 

- 

- 

9 97 +O .73 
- 0.66 

604/481 



C+BB+PL 

073^129 

-1 07+ 0-11 
1U -0.27 

- 

41.4 + 3.2 

1 85+0-0® 

-i.85_o 25 

- 

- 

- 

- 

1.12 + 0.50 

592/480 



C+BB+C 

068^52 

-1.09+0.23 

- 

42.9 + 3.4 

1 «8 + 0 - ^® 
-l.D 8_ 0 Q6 

226 1 + 102.0 
44O.i _ 60 0 

- 

- 

- 

1 52+ 0 - 71 
1,0 -0.65 

581/479 
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TABLE A4 

Parameters of the various tested models resulting from the coarse-time spectral analysis of GRB 090926A with their 

ASYMMETRICAL 1-0 UNCERTAINTIES. 
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TABLE A5 

Parameters of the various tested models resulting from the fine-time spectral analysis of GRB 080916C with their 

ASYMMETRICAL 1 -a UNCERTAINTIES. 


Time from T 0 

Models 

Standard Model 

Additior 

al Model 

Cstat/dof 



C or Band 

BB 

PL 


O-'start 

^stop 

Parameters 

Epeak 

a 

0 

kT 

a 


(s) 

(s) 


(keV) 



(keV) 



- 0.1 

+0.7 

Band 

4f)8 + 162 

408 -110 

-0 77 +°' 17 
u ' -0.13 

9 n 7 +0.29 

L -1.14 



492/473 



C + BB 

2!07«g° 

i 91 +0.13 

_A ^ A -o.n 

- 

49.2 + 8.4 


488/472 



C + BB + PL 5params 

632+gf 

-0.7 (fix) 

- 

44 8+ 20 ’ 5 
44 -11.4 

15 (fix) 

492/472 

+0.7 

+ 1.2 

Band 

400+62 

438 -89 

0 58+ 0 - 14 
-0.5S_o.Q9 

<-2.5 

- 

- 

483/473 



C + BB 

541+ 317 

041 —141 

-0.70 + 0.21 

- 

43 0 + 12 5 

- 

481/472 



C + BB + PL5 p arams 

546±g 7 

-0.7 (fix) 

- 

43.1 + 14.0 

-1.5 (fix) 

482/472 

+ 1.2 

+ 1.7 

Band 

i 473 i 2 l° 

-0.89+0.05 

<-5 

- 

- 

467/473 



C + BB 

I 999+ 394 

i»a»_ 3 5 0 

-0.88+0.05 

- 

35- 611:1 

- 

456/472 



C+BB+PL 

1896+366 

-0 79 + 0 ’ 51 
u ' y —0.25 

- 

35.5 +^ 8 

40 + 0.15 
-1 ' 48 —0.23 

455/470 



C + BB + PL 5params 

1836+| 4 ° 

-0.7 (fix) 

- 

35-3 +il 

-1.5 (fix) 

456/472 

+ 1.7 

+2.3 

Band 

267+51 

n qq+ 0.19 
-0.3 8 _o 14 

2 1 5 + 0-20 
0.39 

- 

- 

533/473 



C + BB 

S =+?42 

-0.73+0.19 

- 

34.3+6.8 

- 

533/472 



C+BB+PL 

4 s 4 A§5 

-0.49+0.14 

- 

32.3 + 3.5 

1 55 + 0-14 
a ' oo _ 0.20 

532/470 



C + BB + PL5 p arams 

S 3 S + H 1 

-0.7 (fix) 

- 

04 c + 5 • 4 
34 ’ 8 — 4.8 

-1.5 (fix) 

533/472 

+2.3 

+ 2.8 

Band 

590+ 110 
oyu — 84 

-0.84+0.07 

<-5 

- 

- 

472/473 



C + BB 

1103^” 3 ^ 3 

iiUO-294 

-1.00+0.08 

- 

38.3+6.3 

- 

463/472 



C+BB+PL 

953+329 
yo -211 

-0 54 + 1 ’ 65 
u o -0.56 

- 

33 0 + 5.2 
33 - 8 _4.o 

-2 14+0-20 
ZA -0.26 

460/470 



C + BB + PL5 p arams 

766+ 242 
/DD — 162 

-0.7 (fix) 

- 

32.6 + 5.4 

-1.5 (fix) 

569/472 

+ 2.8 

+3.4 

Band 

492 + 130 

* y -89 

-0.75+0.09 

<-2.09 

- 

- 

482/473 



C + BB 

qk «+411 

856 -245 

-1.00+0.13 

- 

53.5+10.7 

- 

478/472 



C+BB+PL 

694+ 194 

Dy -135 

-0.63+0.13 

- 

4 °‘5 — 5.9 

1 57 +0.09 
1-57 —0.15 

476/470 



C + BB + PL5 p arams 

64 s) 1 ?!o 

-0.7 (fix) 

- 

47 .5±1 ° 8 8 

-1.5 (fix) 

477/472 

+3.4 

+3.8 

Band 

356±ig 

0 55+0-37 
-U.55 _ 0 15 

+0.19 
^ A -0.38 

- 

- 

478/473 



C + BB 

789+ 4 ® 2 
' oy _225 

-0.75 + 0.16 

- 

28 '9 ±f° 3 

- 

472/472 



C+BB+PL5params 

81 n+ 33 7 

810 -195 

-0.7 (fix) 

- 

2°' 7 — 4.6 

-1.5 (fix) 

471/472 

+3.8 

+4.3 

Band 

eistg? 

-0.80+0.09 

9 1 q + 0.19 
-^ A -0.36 

- 

- 

517/473 



C + BB 

is 47 ±!|| 

-0.99+0.07 

- 

3 4 - 9 ±5:J 

- 

512/472 



C+BB+PL 


-0.87+0.08 

- 

34 5+ 4 ’ 8 
34 ' 5 —4.0 

1 53 +°' 14 
A ’ t>cS — 0.21 

511/470 



C + BB + PL5 p arams 

1379+315 

' y — 408 

-0.7 (fix) 

- 

33 4 + 5 ’ 4 

^ 13 - 5.0 

-1.5 (fix) 

512/472 

+4.3 

+4.8 

Band 

1.570+774 

i&7U _549 

-0.96+0.07 

-1 97+ 016 
A ' y ^ — 0.24 

- 

- 

475/473 



C + BB 

5326 + ^300 
00 Z 0 — 1100 

-1.09+0.05 

- 

49 3+10 - 3 
y ’ - 8.1 

- 

471/472 



C+BB+PL 

2848+1040 

284 «_ 955 

-0.70+0.10 

- 


-1 45 + 0-07 
A ' 4 ° — 0.12 

469/470 



C + BB + PL5 p arams 

3212^” 782 

-703 

-0.7 (fix) 

- 

4 °' 3 — 6.0 

-1.5 (fix) 

470/472 

+4.8 

+5.5 

Band 

«23±??S 

- 0 . 86 + 0.11 

-1 89"*" 013 
Aoy — 0.21 

- 

- 

543/473 



C + BB 

1902 + “°° 

-1.07+0.08 

- 

44 '° —8.4 

- 

544/472 



C+BB+PL 

1325 ttll 

-0 60"*"°' 16 
uou -0.13 

- 


-1.48+0.05 

536/470 



C + BB + PL 5params 

1436+ 4 ^® 

i4oo_ 38 9 

-0.7 (fix) 

- 

4 2.4 +°; 8 

-1.5 (fix) 

537/472 

+5.5 

+ 6.2 

Band 

oopo + 1700 
oouz _ 1230 

-1.34+0.03 

2 1 5 + 0-28 
215 -0.77 

- 

- 

470/473 



C+BB+PL 

1983 _ 662 

-0.70+0.10 

- 

34 +^ 9 

-1.68+0.06 

462/471 



C + BB + PL 5params 

1322+ 7 ° 4 

-0.7 (fix) 

- 

q, 7 + 21.8 
* A -10.0 

-1.5 (fix) 

470/472 

+6.2 

+7.0 

Band 

73 «±SK 

-1.22+0.06 

-1 98+ 0,18 
Ayo —0.38 

- 

- 

494/473 



C + BB 

4726 ± 3 il§S 

-1.40+0.05 

- 

39.6t?;° 

- 

491/472 



C+BB+PL 

1155 +H 7 

-0.70+0.10 

- 

34 '° — 6.1 

1 ««+0.08 
- 1(36 -0.15 

489/471 



C + BB + PL 5params 

623+ 73 | 

-0.7 (fix) 


31- 5 —9.4 

-1.5 (fix) 

495/472 
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Time from T 0 

Models 

Standard Model 

Additior 

al Model 

Cstat/dof 



C or Band 

BB 

PL 


Tstart 

Tstop 

Parameters 

Epeak 

0 

0 

kT 

a 


(s) 

(s) 


(keV) 



(keV) 



+7.0 

+7.5 

Band 

509+™ 

-1.07+0.14 

<-1.75 



467/473 



C + BB 

1788^1900 

-1.28 + 0.12 

- 

39 '0ii5.i 

- 

465/472 



C+BB+PL 

X274+- 0 

-0 94 + 0 - 19 
uy4 —0.14 

- 

34 +^°5 3 

-1 67+ 012 
± D -0.24 

465/470 



C + BB + PL5 p arams 

549 1? 46 

-0.7 (fix) 

- 

27.8 ±l 2 2 3 

-1.5 (fix) 

468/472 

+7.5 

+7.9 

Band 

5382 t“tS 

-1.34+0.06 

<-1.70 

- 

- 

534/473 



C + BB + PL5 p arams 

44 + 12 
* -9 

-0.7 (fix) 

- 

91 3±1| § 

-1.5 (fix) 

527/472 

+7.9 

+8.5 

Band 

535+305 

03 -141 

- 1 . 11 + 0.10 

<-1.91 

- 

- 

492/473 



C + BB + PL5params 

35±& 

-0.7 (fix) 

- 

61.9±1° 9 3 

-1.5 (fix) 

487/472 

+8.5 

+9.5 

Band 


, 1 o+0.13 
-±.io_o 12 

1 09+013 

-1.S4 —0 32 

- 

- 

530/473 



C + BB 

5496~*”® 339 

o^yo_2740 

-1.41+0.06 

- 

37 6”*" 9 ' 4 

- 

526/472 



C+BB+PL 

i oqo+2350 

1382 _746 

-0.60+0.10 

- 

34 1 + 7 ’ 7 
-6.9 

-1 64+ 0 ' 29 
1D - 0.20 

526/471 



C + BB + PL 5params 

798+ 1050 
/yo —395 

-0.7 (fix) 

- 

32 -5±l;l 

-1.5 (fix) 

530/472 

+9.5 

+ 11.5 

Band 

ooi +206 

3 -105 

-1.09+0.11 

-1 69"*" 007 
1 Dy - 0.10 

- 

- 

534/473 



C + BB 

79n o +6300 
7208_273 q 

-1.38+0.05 


34 -5 ±11 

- 

533/472 



C+BB+PL 

679+ 341 

o ' y -200 

-0.63+0.15 

- 

27 ..±»;3 

-1.52+0.03 

533/471 



C + BB + PL5 p arams 

,-24+500 

524 -204 

-0.7 (fix) 

- 

28.4+7 

-1.5 (fix) 

534/472 

+ 11.5 

+ 14.0 

Band 

618+™ 

-1 04 - *" 0-09 
1U -0.06 

<-2.06 

- 

- 

527/473 



C + BB 

1289^^ 

- 1 . 20 + 0.10 

- 

40.8+9.5 

- 

524/472 



C+BB+PL 

831+36 

-0.66+0.17 

- 

34.61 3 ’ 8 

1 ficj+0-09 

-i.o5 —q 17 

525/471 



C + BB + PL5 p arams 

51 + 2 138 

-0.7 (fix) 

- 

31.8+ 3 ’ 4 
— 6.5 

-1.5 (fix) 

527/472 

+ 14.0 

+ 18.0 

Band 

429+ 109 
— 81 

-1.06+0.06 

1 07 +O.O 9 
- 1 ’ 87 -0.14 

- 

- 

490/473 



C + BB 

(?9qo+6690 

o2oo_292o 

-1.42+0.05 

- 

47.9 + 5.7 

- 

490/472 



C+BB+PL 

471 + 190 

4 ' 1 — 109 

-0.91+0.12 

- 

48 ' 8 — 7.5 

-1.52+0.08 

487/470 



C + BB + PL5 p arams 

coo+ 28 9 

333 — 142 

-0.7 (fix) 

- 

33 

-1.5 (fix) 

488/472 

+ 18.0 

+ 22.0 

Band 

3 "±S 

-0.97+0.05 

<-2.27 

- 

- 

482/473 



C + BB 

oooitfi 

-1.23+0.12 

- 

48 .3+9 ; 3 

- 

477/472 



C+BB+PL 


-0.79+0.09 

- 

42 ' 4 ±4.9 

-1.71+0.06 

476/470 



C + BB + PL5 p arams 

391 till 

-0.7 (fix) 

- 

40.5± 34 4 9 

71 +o.io 

L -' - 0.20 

476/471 

+ 22.0 

+26.0 

Band 

« 4 tSS 

-0.95+0.05 

2 90 + O. 17 
z z -0.34 

- 

- 

486/473 



C + BB 

•5QKQ + 1560 
3053_ H40 

-1.41+0.05 

- 

52.3 + 3.1 

- 

465/472 



C+BB+PL 

<™±5IS 

-0.71+0.06 

- 

47 - 8 + 4 -2 
— 3.6 

-1.68+0.07 

468/470 



C + BB + PL 5params 

685 ±K 2 7 

-0.7 (fix) 

- 

i7S ±tl 

-1.68+0.07 

468/471 

+ 26.0 

+32.0 

Band 

489+ I5 

— 65 

-1.10+0.04 

<-2.23 

- 

- 

456/473 



C + BB 

so, +877 

331 — 250 

-1.29+0.10 

- 

4 6-6 ±11 

- 

451/472 



C+BB+PL 


-0.76+0.09 

- 

37 2 + 4 ’ 1 

4 ' ' z — 3.7 

-1.77+0.05 

448/470 



C + BB + PL5params 

527 1 1 " 

-0.7 (fix) 

- 

3 6.4 + 5;7 

-1 76+ 0,09 
A- ' D -0.17 

448/471 

+32.0 

+40.0 

Band 

4 56+t 

-1.08+0.04 

<-5 

- 

- 

535/473 



C + BB 

.+424 

664 -201 

-1.33+0.09 

- 

55.+?;? 

- 

525/472 



C + BB + PL 5params 

502 j: 33 ^ 

-0.7 (fix) 

- 

43 - 2 ±? 2 .I 

-1.73+0.13 

524/471 

+40.0 

+53.0 

Band 

482 + 8 5 

-1.08+0.05 

<-5 

_ 


518/473 



C + BB 

1 4ns+ 754 
1408 _597 

-1 39 + 0 - 09 
lc5y - 0.06 

- 

48 ' 3 —1. 1 


506/472 



C+BB+PL 

cl q + 178 

618 -117 

-0.66+0.13 

- 

36.5 + 3.2 

- 1.90+0.12 

504/470 



C + BB + PL5 p arams 

636+ 354 

03 -169 

-0.7 (fix) 

- 

36.8 + 5.2 

-1 91+0-16 
y -0.26 

504/471 

+53.0 

+71.0 

Band 

204“*" 40 
u -26 

-1.05+0.07 

-2 29 + 0-21 
^ -1.40 

- 

- 

492/473 



C + BB 

oi -+388 
317 —93 

-1.31+0.18 


33 K +10 0 

- 

491/472 



C + BB + PL 5params 

65+ 90 

b5 -29 

-0.7 (fix) 

_ 

41 0 + 12 - 6 
- 6.8 

-1.5 (fix) 

489/472 
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TABLE A6 

Parameters of the various tested models resulting from the fine-time spectral analysis of GRB 090926A with their 

ASYMMETRICAL 1-a UNCERTAINTIES. 


Time from Tq 

Models 

Standard Model 

Additional Model 

Cstat/dof 



C or Band 

BB 

PL 01 

C 


^start 

^stop 

Parameters 

E peak 

a 

0 

kT 

a 

E 0 


(s) 

(s) 


(keV) 



(keV) 


(MeV) 


+ 0.0 

+ 1.5 

Band 

326+ 112 

-0.71+0.13 

9 pio+0- 1 4 

-2.08_o.32 




609/566 



C + BB 

551 + 1320 

±t>Oi_79o 

+ 0.17 
1 Z 1 - 0.10 


44.0 + 5.5 

- 

- 

605/565 



C+BB+PL 5params 

561±f 2 | 

-0.7 (fix) 

- 

37.6 + 7.0 

-1.5 (fix) 

- 

607/565 

+ 1.5 

+ 1.8 

Band 

302+« 

-0.43+0.12 

-1.90+0.08 

- 

- 

- 

597/566 



C + BB 


-1.08+0.06 

- 

43.9+4.0 

- 

- 

595/565 



C+BB+PL5params 

*»i±SS 

-0.7 (fix) 

- 

42.4 + 3.7 

-1.5 (fix) 

- 

599/565 

+ 1.8 

+ 2.0 

Band 

2S5 ±it 

-0.29+0.13 

-2.24 + 0.16 

- 

- 

- 

570/566 



C + BB 

667 + 151 

DD —105 

-0.72+0.11 

- 

36.6+4.7 

- 

- 

562/565 



C + BB 

fjCI + 83 

° &1 _69 

-0.7 (fix) 

- 

36.0+3.3 

- 

- 

562/566 



C+BB+PL 5params 

6 S 2 + 1° 8 

-0.7 (fix) 

- 

37.3+3.5 

-1.5 (fix) 

- 

561/565 

+ 2.0 

+ 2.3 

Band 

305±f 8 

-0.25+0.09 

_2 47 + 0 ' 14 
- 0.21 

- 

- 

- 

599/566 



C + BB 


-0.56+0.09 

- 

36.2 + 4.2 

- 

- 

585/565 



C + BB 

«38±“ 

-0.7 (fix) 

- 

40.9 + 2.8 

- 

- 

587/566 

+2.3 

+ 2.4 

Band 

324+6| 

-0.47+0.13 

-2 30 - *" 0-18 
z ‘ au -0 .35 

- 

- 

- 

568/566 



C + BB 

500^67 

0 55 + 0-14 
-0.55_o ii 


23 + 4.3 

- 

- 

564/565 



C + BB 

660±1| 

-0.7 (fix) 

- 

29.sjl6; 2 

- 

- 

566/566 



C+BB+PL 5params 

660 tH 

-0.7 (fix) 

- 

32.818.9 

-1.5 (fix) 

- 

564/565 

+ 2.4 

+ 2.5 

Band 

270 ±t 3 

-0.24 + 0.13 

9 97 + 0 .13 
-0.20 

- 

- 

- 

567/566 



C + BB 

‘«±S 

-0 44 + 0 15 
0,4 -0.13 


27.4±5 : * 

- 

- 

565/565 



C + BB 

568t'l 

-0.7 (fix) 

~ 

36.6+31 

- 

- 

567/566 



C+BB+PL5params 

590+88 
oyu —68 

-0.7 (fix) 

- 

38.4lt;8 2 

-1.5 (fix) 

- 

567/565 

+ 2.5 

+ 2.6 

Band 

338+54 

-0 46+ 013 
U4D —0.17 

- 2 - 51 -2.12 

- 

- 

- 

539/566 



C + BB 

*^+124 

566_si 

-0.69+0.12 

~ 

34.6l?;I 

- 

- 

533/565 



C + BB 


-0.7 (fix) 

- 

35.018.0 

- 

- 

533/566 



C+BB+PL5params 

666 i™ 

-0.7 (fix) 

- 

34.6l®;“ 

-1.5 (fix) 

- 

533/565 

+ 2.6 

+ 2.7 

Band 

268+66 

-0.37+0.13 

- 2 * 2 ±°o.ll 

- 

- 

- 

595/566 



C + BB 

662+ 147 

00 -104 

-0.78+0.11 

~ 

35.6+4.2 

- 

- 

581/565 



C + BB 


-0.7 (fix) 

~ 

33 2 ^” 3-3 
dd "^_2.9 

- 

- 

576/566 

+ 2.7 

+ 2.9 

Band 

304 ttt 

-0.49+0.09 

-2 48+ 0 ' 16 
—0.28 

- 

- 

- 

581/566 



C + BB 

668 ± 1 ° 

-0.76+0.09 

- 

33.8+4.7 

- 

- 

573/565 



C + BB 

493tg 

-0.7 (fix) 

- 

33.41 3 2 

- 

- 

573/566 

+ 2.9 

+3.0 

Band 


-0.51+0.06 

<-5 

_ 

_ 

_ 

545/566 



C + BB 

„„+52 

660-45 

n qq+ 0-20 

-0' 38 -0.32 

- 

2 °- 4 1 !;o 

- 

- 

542/565 



C + BB 

678 + g8 

-0.7 (fix) 

- 

88 . 0 I 42 8 

- 

- 

547/566 

+3.0 

+ 3.1 

Band 

520+q6 

-0.57+0.07 

-3 09+ 0 52 
Juy —2.41 

- 

- 

- 

659/566 



C + BB 

r + 164 
645 -105 

-0.72+0.14 

- 

89’ 4 1 40 7 

- 

- 

658/565 



C + BB 

620 tit 

-0.7 (fix) 

- 

68.0143 4 

- 

- 

659/566 



C+BB+PL5 p arams 

652 + ™ 

-0.7 (fix) 

- 

eo q+16.3 
oo-o_ IQ s 

-1.5 (fix) 

- 

657/565 

+3.1 

+3.2 

Band 

407+45 

-0.38+0.08 

_ 2 70+0-25 

'-0.40 

- 

- 

- 

536/566 



C + BB 

579+^7 

579 — 60 

-0.53+0.10 

- 

36.0+7.5 

- 

- 

532/565 



C + BB 

687+ 84 
— 56 

-0.7 (fix) 

- 

44.916.8 

- 

- 

535/566 



C+BB+PL5 p arams 

696 + 1? 

— 56 

-0.7 (fix) 

- 

48.6l5;6 

-1.5 (fix) 

- 

534/565 

+3.2 

+3.3 

Band 

362+31 

-0.32+0.08 

-3 Q 1+ 0,33 
— 0.62 

- 

- 

- 

536/566 



C + BB 

490 _199 

-0.63+0.18 

- 

57.7+11.0 

- 

- 

536/565 



C + BB 

623+6g 

-0.7 (fix) 

- 

68.7l 7 6 ;I 

- 

- 

537/566 



C+BB+PL5 p arams 

632±J 9 2 

-0.7 (fix) 

- 

68 . + 5'.8 

-1.5 (fix) 

- 

536/565 

+3.3 

+ 3.4 

Band 

371+33 

-0.44+0.07 

3 1 ej+0-36 
J ‘ LO —0.82 

- 

- 

- 

511/566 



C + BB 

484+ 89 
— 65 

-0.63+0.14 

- 

44 0 - * -9,8 
^ -13.2 

- 

- 

510/564 



C + BB 

518+ 44 

518 -40 

-0.7 (fix) 

- 

47 4 + 7 ' 4 

- 

- 

510/566 

+3.4 

+3.5 

Band 

333+54 

-0.50+0.11 

-2 69~*” 9-39 
^ ,oy —0.73 

- 

- 

- 

490/566 



C + BB 

M O + 80 

5A>5_60 

-0.69+0.10 

- 

3i.9l|;| 

- 

- 

482/565 



C + BB 

519+ 50 

0_Ly — 43 

-0.7 (fix) 

- 

09 o+4.0 
^'•’-3.4 

- 

- 

482/566 



C+BB+PL5 p arams 

521 + 53 
521 -44 

-0.7 (fix) 

- 

3 2 -6i 4 3 : 2 8 

-1.5 (fix) 


482/565 
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S.Guiriec 


Time from T 0 

Models 

Standard Mod 

b! 

Additional Model 

Cstat/dof 



C or Band 

BB 

PL or 

C 


'I'start 

Tstop 

Parameters 

E peak 

a 

P 

kT 

a 

E 0 


( 8 ) 

(s) 


(keV) 



(keV) 


(MeV) 


+3.5 

+3.6 

Band 


-0.39+0.09 

-2 43 + 0 - 14 
— 0.20 

_ 

_ 

_ 

577/566 



C+BB 

536±™ 

-0.58+0.10 

- 

33.3+I ' 8 

- 

- 

575/565 



C+BB 

598 t|o 

-0.7 (fix) 

- 

35.5+i't 
— 3.6 

- 

- 

576/566 



C + BB + PL5params 

603”^”®^ 

*> Ud _50 

-0.7 (fix) 

- 

36.2±fl 

-1.5 (fix) 

- 

576/565 

+3.6 

+3.7 

Band 

323+“ 

-0.52+0.10 

<-2.64 

- 

- 

- 

570/566 



C + BB 

420^“ 

-0.60+0.10 

- 

26.2+|;4 

- 

- 

563/565 



C + BB 

4531“ 

-0.7 (fix) 

- 

30.7±S;« 

- 

- 

564/566 

+3.7 

+3.8 

Band 


-0.39+0.10 

9 97 + 0.11 
z z -0.15 

- 

- 

- 

569/566 



C + BB 


-0.87+0.12 

- 

49 S + 4 ' 7 
42 ’ 8 —5.0 

- 

- 

569/565 



C + BB 

5531“ 

-0.7 (fix) 

- 

37.2 + 4 4 ; 8 

- 

- 

572/566 



C + BB + PL5params 

6001“ 

-0.7 (fix) 

- 

46-4^3 § 

-1.5 (fix) 

- 

567/565 

+3.8 

+3.9 

Band 

271 ±25 

-0.38+0.09 

o ^4+0.17 

-2.54 — 0 26 

- 

- 

- 

551/566 



C + BB 

626 -143 

-1.00+0.15 

- 


- 

- 

547/565 



C + BB 


-0.7 (fix) 

- 

41.2 + 8 ; 7 

- 

- 

550/566 



C + BB + PL 5params 

4381“ 

-0.7 (fix) 

- 

43.0+11 

-1.5 (fix) 

- 

545/565 

+3.9 

+4.0 

Band 

265+“ 

-0.39+0.10 

9 qc + 0-13 
-2.3 8 _o 19 

- 

- 

- 

622/566 



C + BB 

««±ss 

-0.74+0.11 

- 

356il“ 

- 

- 

621/565 



C + BB 

4641“ 

-0.7 (fix) 

- 

34.oi 4 ; 2 5 

- 

- 

621/566 



C + BB + PL5 p arams 

48 °it! 

-0.70 (fix) 

- 

36.3+i; 9 

-1.5 (fix) 

- 

617/565 

+4.0 

+4.1 

Band 

220 l 24 

-0.25+0.10 

-2.16+0.08 

_ 

_ 

_ 

560/566 



C + BB 

QQO + 1 94 

893_i 34 

-0 96"*" 012 
uyo — 0.20 

- 

35.9±!;! 

- 

- 

561/565 



C+BB+PL 

746lf 4 8 

0 OO + 0.31 

- 0 - 82 _o i4 

- 

34.7+3.2 

-i ck+0.34 
1 ’ £>£) —0.19 

- 

560/563 



C + BB 

878 —53 

-0.7 (fix) 

- 

30-0±2.J 

- 

- 

566/566 



C + BB + PL5params 

626"*" 88 

-0.7 (fix) 

- 

32.6 + 2.4 

-1.5 (fix) 

- 

560/565 

+4.1 

+4.2 

Band 

268+“ 

-0.44+0.09 

-2.25 + 0.11 

- 

- 

- 

576/566 



C + BB 

772 t^ 

-0 93+ 0 11 
u,y -0.18 

- 

35.4+ 4 -? 

- 

- 

601/565 



C+BB+PL 

60911“ 

-0 47 + 0 52 
0,4 -0.32 

- 

31 7 + 3 ' 8 
d± '— 3.1 

1 65 + 0 - 13 
-i. 66-0 09 

- 

595/563 



C + BB 

544 t« 

-0.7 (fix) 

- 

28.5 + 2.4 

- 

- 

576/566 



C + BB + PL5p arams 

5731“ 

-0.7 (fix) 

- 

31 2 "*" 2-8 
^ - 2.6 

-1.5 (fix) 

- 

571/565 

+4.2 

+4.3 

Band 

2221 “ 

-0.25+0.11 

-2.11+0.07 

- 

- 

- 

545/566 



C + BB 

-"liss 

-0.83+0.11 

- 

33.5 + 3.3 

- 

- 

560/565 



C+BB+PL 

573l“6 

-0 63 - * -0 ’ 29 
u.Od-o is 

- 

32.2 + 3.3 

-1 49 + 0 - 19 

1 ^ -0.09 

- 

549/563 



C + BB 

5701“ 

-0.7 (fix) 

- 

30.0+2.1 

- 

- 

561/566 



C + BB + PL 5params 

620 + 11 
— 63 

-0.7 (fix) 

- 

33.5 + 2.4 

-1.5 (fix) 

- 

550/565 

+4.3 

+4.4 

Band 

211 + 20 

-0.25+0.10 

-2.15+0.08 

- 

- 

- 

563/566 



C + BB 

619+ 180 
oiy —106 

-0.87+0.12 

- 

34.0 + 3.8 

- 

- 

577/565 



C+BB+PL 

49811“ 

0.65-0 2 i 

- 

32.0+4.2 

-1 45 + 0 ' 26 
— 0.12 

- 

565/563 



C + BB 

5041“ 

-0.7 (fix) 

- 

29.5 + 2.3 

- 

- 

579/566 



C + BB + PL5params 

847 ±S 

-0.7 (fix) 

- 

33.6+ - 2.5 

-1.5 (fix) 

- 

566/565 

+4.4 

+4.5 

Band 

293+35 

-0.55+0.08 

-2.35+0.15 

- 

- 

- 

590/566 



C + BB 

55211“ 

- 0 . 88 + 0.11 

- 

37.0 ± 5 6 ; 9 7 

- 

- 

596/565 



C+BB+PL 

5031113 

0 «q+ 0*20 
-0.S3-0 12 

- 

36-361 

-1 34 + 0 - 73 
^ - 0.20 

- 

593/563 



C + BB 

44011® 

-0.7 (fix) 

- 

97 8^~ 8 ’ 9 
27 ' 8 —4.8 

- 

- 

599/566 



C + BB + PL5 p arams 

4551“ 

-0.7 (fix) 

- 

32 7 + 5 ' 5 
-4.7 

-1.5 (fix) 

- 

593/565 

+4.5 

+4.7 

Band 

314+ 27 

■ 514 -24 

-0.59+0.06 

-2 59 + 0 15 
2 .t ,y _o 21 

- 

- 

- 

571/566 



C + BB 

450+ 88 

-0.80+0.09 

- 

37. 1 +6 7 ; 7 

- 

- 

574/565 



C+BB+PL 

432+ 62 

4 -37 

-0 74 + 0 - 20 
-0.28 

- 

35.61!;! 

1 55 + 0 ’ 14 

i. 66-0 22 

- 

573/563 



C + BB 

405+23 

-0.7 (fix) 

- 

29 7 + 7 ' 2 
zy ' -5.3 

- 

- 

576/566 



C + BB + PL5params 

411 + 25 

-0.7 (fix) 

- 

32 - 9 l!;! 

-1.5 (fix) 

- 

574/565 

+4.7 

+4.9 

Band 

295+27 

-0.60+0.06 

-2.17+0.08 

- 

- 

- 

612/566 



C + BB 

61711“ 

-0.98+0.10 

- 

42 . 6+!;6 

- 

- 

634/565 



C+BB+PL 

502+® 7 

n S4+0-17 
-0.S4-0 is 

- 

41.5+!! 

— 6.3 

-1 43 + 0 ' 21 
1,4 - 0.10 

- 

614/563 



C + BB 

429+30 

-0.7 (fix) 

- 

2 4.1+ 8 ;6 

- 

- 

643/566 



C + BB + PL5params 

462+40 

-0.7 (fix) 

- 

38.1+6;! 

-1.5 (fix) 

- 

615/565 

+4.9 

+5.2 

Band 

310+30 

-0.63+0.06 

-2.06+0.06 

- 

- 

- 

625/566 



C + BB 

2072"^ 94 

U ' -464 

-1.24+0.05 

- 

45.3 + 2.5 

- 

- 

638/565 



C+BB+PL 

071 + 282 
871 -178 

-1 04+ 0 - 15 
± u - 0.10 

- 

42.9 + 3.2 

-1 40+ 0 - 23 
1,4U — 0.12 

- 

622/563 



C + BB 

467+33 

-0.7 (fix) 

- 

26.6+5.5 

- 

- 

686/566 



C + BB + PL5 p arams 

590^” 72 
oyu —58 

-0.7 (fix) 


38.ol!;! 

-1.5 (fix) 


628/565 
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Time fr 

om Tq 

Models 


Standard Mod 

el 

Additional Model 

Cstat / dof 



C or Band 

BB 

PL or 

C 


^start 

Tstop 

Parameters 

E peak 

a 

P 

kT 

<* 

E 0 


(8) 

(s) 


(keV) 



(keV) 


(MeV) 


+5.2 

+5.6 

Band 

390+30 

-0.68+0.04 

-2 39+®’ 11 

z ‘ av -0.14 

_ 

_ 

_ 

609/566 



C + BB 

709+ 154 
' uy —109 

-1.01+0.08 

- 

51.5 + 24.5 

- 

- 

615/565 



C+BB+PL 

592+ 125 
oyz —95 

-0 90+ 0 - 25 
u y -0.13 

- 

er. « + 5.7 

50.6 g 3 

44 +0.55 
± ’ 44 —0.16 

- 

603/563 



C + BB 

446+21 

-0.70 (fix) 

- 

84 q+0-5 
34 ’ 3 — 11.4 

- 

- 

630/566 



C + BB + PL5 p arams 

487+31 

-0.7 (fix) 

- 

44 1 + 7 - 2 
' -5.9 

-1.5 (fix) 

- 

605/565 

+5.6 

+6.0 

Band 

338+30 

-0.68+0.05 

-2.25+0.10 

- 

- 

- 

594/566 



C + BB 

574+122 

-0.97+0.09 

- 

47.0+5.5 

- 

- 

611/565 



C+BB+PL 

420"*" 80 
4 ^ u —50 

-0.60+0.30 

- 

39.0+10.5 

1 Kq + 0.16 

± ' £>d —0.07 

- 

587/563 



C + BB + PL5 params 

435 + 31 

-0.7 (fix) 

- 


-1.5 (fix) 

- 

587/565 

+6.0 

+6.4 

Band 

433+42 

-0.87+0.04 

9 oc+0.13 
-Z.oo — o 20 

- 

- 

- 

628/566 



C + BB 

«i±Jg 3 

-1.01+0.09 


59.3+”;S 

- 

- 

636/565 



C + PL 

«e±ff 

-0.70+0.14 

- 

- 

1 55+ 015 
± - ot> —0.07 

- 

621/565 



C+BB+PL 

44 «i! 4 7 

-0.81+0.22 

- 

k R q+26.4 
58-9 —14.4 

-1.52 + 0.04 

- 

621/563 



C+PL 

402 + 20 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

622/567 



C + BB + PL5 p arams 

402 + 20 

-0.7 (fix) 

- 

40 O"*" 7-4 
uu -5.9 

-1.5 (fix) 

- 

622/565 

+6.4 

+6.7 

Band 

582 + 44 

-0.76+0.04 

9 78 + 0.24 
- 2 78 -0.47 

- 

- 

- 

607/566 



C + BB 

eootH 1 

-0.87+0.08 

- 

63.0±^|| 

- 

- 

608/565 



C+PL 

532+H 

-0.50+0.11 

- 

- 

-1.62 + 0.07 

- 

597/565 



C+BB+PL 


-0.7 (fix) 

- 

52.2 + 31 , 0 

-1.61 + 0.09 

- 

596/564 



C + PL 

580+27 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

601/567 



C + BB + PL5 p arams 

606+40 

-0.7 (fix) 

- 

311 12.6 

-1.5 (fix) 

- 

599/565 

+6.7 

+ 7.2 

Band 

501 + 24 

-0.68+0.03 

<-3.05 

- 

- 

- 

719/566 



C + BB 

527 -48 

-0.86+0.07 

- 

84.9ll 2 8 7 

- 

- 

678/565 



C + PL 

443 + 20 

-0.42+0.09 

- 

- 

71 +0.07 

±ri - 0.11 

- 

679/565 



C+BB+PL 

50+f 4 

-0.79+0.08 

- 

85 7+ 7 ■ 7 
85-7 — 10.8 

-1.54+0.16 

- 

676/563 



C + BB 


-0.7 (fix) 

- 

“3 tit 

- 

- 

685/566 



C + PL 

5OO+1 0 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

691/567 



C + BB + PL5 p arams 

473+ 32 
* -36 

-0.7 (fix) 

- 

86-5±i?:? 

-1.5 (fix) 

- 

675/565 

+7.2 

+ 7.9 

Band 

243+13 

-0.59+0.04 

9 71+0.16 
" z - '-0.22 

- 

- 

- 

634/566 



C + BB 

39S±“ 

-0.98+0.08 

- 

38.2 + 2.2 

- 

- 

627/565 



C+PL 

244+14 

-0.41+0.13 

- 

- 

-1 61+ 009 
1-01 — 0.06 

- 

637/565 



C+BB+PL 

287 

— 56 

-0.58+0.40 

- 

36.4+5.0 

-1.41 + 0.05 

- 

624/563 



C + BB 

292+11 

-0.7 (fix) 

- 

33i±l;i 

- 

- 

639/566 



C+PL 

288+ 8 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

648/567 



C + BB + PL5 p arams 

296+14 

-0.7 (fix) 

- 

»»-»±i.5 

-1.5 (fix) 

- 

627/565 

+7.9 

+8.0 

Band 

H>7±“ 

n 4R+0-23 
-0.48 _ q 19 

+0.12 

Zi4 -0.17 

- 

- 

- 

531/566 



C + BB 

3-4+5? 

1 O1+0.26 

1Z1 —0.22 

- 

20.0+|^ 

- 

- 

534/565 



C+PL 

128 + 20 

-0 47+ 0 ' 49 
u ’ 4 -0.28 

- 

- 

-1 54+ 0-22 
10 -0.12 

- 

534/565 



C+BB+PL 

186+J? 

-0.73+0.20 

- 

iT.ajl 3 ; 3 

-1.58 + 0.11 

- 

532/563 



C + PL 

142+11 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

535/567 



C + BB + PL5 p arams 

l®7 + 27 

-0.7 (fix) 

- 

16 7+3-3 
‘-3.4 

-1.5 (fix) 

- 

532/565 

+8.0 

+8.2 

Band 

343 + 34 

-0 72 + 0 07 
U- ' -0.04 

<-4.03 

- 

- 

- 

568/566 



C + BB 

„ s1 +85 

381 -58 

-1.00 + 0.14 

- 

67.7+®-8 

— 6.3 

- 

- 

562/565 



C + BB + PL5 p arams 

329+60 

-0.7 (fix) 

- 

«+14.6 

566 -10.0 

-1.5 (fix) 

- 

563/565 

+8.2 

+8.8 

Band 

288+13 

-0.60+0.03 

-3 60~*~°' 44 
4.ou_i 20 

- 

- 

- 

628/566 



C + BB 

358^44 

-0.93+0.09 

- 

49 - 5 ±l:l 

- 

- 

616/565 



C + PL 

261 + 13 

-0.21+0.13 

- 

- 

+0.07 
— 0.10 

- 

614/565 



C+BB+PL 

302t 3 4 

-0.32+0.25 

- 

36.8l 3 -^ 

, o-i +0.08 

_1-81 — 0.14 

- 

610/563 



C + BB 

302 + 16 

-0.7 (fix) 

- 

51 o + 1 0 • 3 
51 - 2 -7.8 

- 

- 

622/566 



C + BB + PL5 p arams 

298+19 

-0.7 (fix) 

- 

51 -°±6.3 

-1.5 (fix) 

- 

617/565 

+8.8 

+9.5 

Band 

237+12 

-0.67+0.04 

-2 94+ 0,20 
^ y4 -0.30 

- 

- 

- 

643/566 



C + BB 

386+J 3 

-1.07+0.08 

- 

37.0+2.3 

- 

- 

628/565 



C+BB+PL 

359+19 

-0.70+0.10 

- 

33.9+3.1 

-1.55+0.08 

- 

626/563 



C + BB 

256 + 9.1 

-0.7 (fix) 

- 

29 6+ 11 - 9 
y- -6.4 

- 

- 

648/566 



C + BB + PL5 p arams 

401+60 

-0.7 (fix) 

- 

34 0+ 7 - 5 
d4-u -6.0 

-1.5 (fix) 

- 

630/565 

+9.5 

+9.7 

Band 

401 + 16 

-0.44 + 0.04 

<-3.58 

- 

- 

- 

565/566 



C + BB 

418+48 

-0.63+0.11 

- 

oq +19.2 
83 ‘ —10.3 

- 

- 

559/565 



C + BB 

453 + 30.0 

-0.7 (fix) 

- 

63 3"*" 7-9 
D4.4_6 4 

- 

- 

559/566 



C + BB + PL 5params 

474+® 1 
4 ' 4 -46 

-0.7 (fix) 


& i +8.1 

61,6 -5.3 

-1.5 (fix) 


559/564 
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S.Guiriec 


Time from T 0 

Models 

Standard Model 


Additional M 

odel 

Cstat/dof 



C or Band 

BB 

PL or C 


^start 

Tstop 

Parameters 

■^peak 

a 

P 

kT 

a 

E 0 


(s) 

(8) 


(keV) 



(keV) 


(MeV) 


+9.7 

+9.9 

Band 

285 + 13 

-0.51 + 0.05 

o ns+ 9 -28 

-3.08_o 47 


_ 

_ 

635/566 



C + BB 

409l|2 

-1.00 + 0.12 

- 

54.9 + 3.0 

- 

- 

624/565 



C+BB+PL 

355 ±Is 

n 8R+ 9-49 

-0.88_ 0 is 

- 

55.4+4 4 3 

1 51+ 0 - 44 
- 1 - 51 —0.17 

- 

607/589 



C + BB 

322 + 22 

-0.7 (fix) 

- 


- 

- 

625/566 



C + BB + PL5 p arams 

315 + 31 

-0.7 (fix) 

- 

48.2++ + 

-1.5 (fix) 

- 

611/565 

+9.90 

+9.95 

Band 

293^ 

-1.21+0.07 

-1.88+0.06 

- 

- 

- 

531/566 



BB + PL 

- 

- 

- 

42.3+3'.) 

-1.63+0.02 

- 

530/566 



BB + C 

- 

- 

- 

40.0+; 3 

-1.58+0.03 

27 640+ 38- ^ 99 

' ,D u _n. 6 o 0 

526/565 



BB + C 

- 

- 

- 

35.0+-5 

-1.5 (fix) 

q fioo+3.950 

9.823_ 2 680 

566/563 



C + BB + PL5 p arams 

2^8 2 

-0.7 (fix) 

- 

+ 2 +>±l 2 0 

-1.5 (fix) 

- 

526/565 



C + BB + C5 params 

22 + 1 2 

-0.7 (fix) 

- 

39.3+;5 

-1.5 (fix) 

23 660 + 18 ’ 900 
jo.ddu _ 8 5 8Q 

519/564 

+9.95 

+ 10.00 

Band 

4SS+153 

488 -101 

-1.34+0.05 

_ 2 07+ 0 - 11 
z v -0.17 

- 

- 

- 

572/566 



BB + PL 

- 

- 

- 

«o.7±S;i 

-1.71 + 0.03 

- 

572/566 



BB + C 

- 

- 

- 

55.61-,“ 

-1.65 + 0.05 

16 870+ 8 -*"® 99 
10.0 iu — g.g 0 0 

569/565 



C+BB+PL 

390+ 121 

oyu _98 

-0.80+0.10 

- 

54.4+-13 

-1.74+0.05 

— 

566/563 



BB + C 


- 

- 

37.0+;! 

-1.5 (fix) 

0 148 + 1-060 
^.I48_ 0 661 

544/566 



C + BB + PL 5params 

30 + 1° 

-0.7 (fix) 

- 

44.3+6.6 

-1.5 (fix) 

- 

546/565 



C + BB+C5 p arams 

21±5 

-0.7 (fix) 

- 

38.i+;! 

-1.5 (fix) 

9 242+ 9 ' 980 
y.z4^_ 4 330 

537/564 

+ 10.0 

+ 10.1 

Band 

329l« 

-1.18+0.05 

46+0.19 

z -^-0.35 

- 

- 

- 

494/566 



C + BB 


-1.31+0.10 


35.6+g 9 ^ 8 

- 

- 

495/565 



C+PL 

295t 33 

n 84+ 0 - 25 

-0.84-0 2i 

- 

- 

-1.79 + 0.09 

- 

492/565 



C + C 


n R k + 0.41 
-0.65_ q 33 

- 

- 

-1.65+0.09 

9 379 + 15 - 700 
y.o' y _4.75o 

490/564 



C+BB+PL 

3 ™+ 

-0.76+0.14 

- 

27.5+4.8 

-l.86_o.l2 

- 

473/564 



C+BB+C 

339 -l 2 6 

-0.76+0.14 

- 

28.8+4.8 

.1 71+0-10 
' — 0.14 

10 340+ 37-999 

1U.04U — 4.930 

471/563 



BB+C 


- 

- 

35.8+;6 

-1.5 (fix) 

•1 977 + 0.502 

L Z ‘ ‘ —0.330 

475/566 



C+C 

171 + 17 

' -14 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

4 953+ 2 - 860 
4.yo4_ 2 060 

479/566 



C + BB + PL5 p arams 

405 + 22® 
uo -128 

-0.7 (fix) 

- 

32.6+3.5 

-1.85 + 0.08 

- 

473/564 



C + BB + C5 params 

26 + 1 6 

-0.7 (fix) 

- 

36 -9±f* 

-1.5 (fix) 

4.878± 2 2 90 7 ° 0 

471/564 

+ 10.1 

+ 10.2 

Band 

197+20 

-0.92+0.07 

o co + 0.20 

-4.oo — o 32 

- 

- 

- 

548/566 



C + BB 

291+8° 

y -46 

-1.22+0.13 

- 

34.7±5.5 

- 

- 

548/565 



C+PL 

188 + 20 

0 c ;q+ 9 - 49 
-u.5y_ 0 32 

- 

- 

-1 70+°-20 
±- ' u —0.08 

- 

546/565 



C+BB+PL 

2S8+H 

-1.10+0.08 

- 

35.5+-“ 

-1.60 + 0.09 

- 

545/563 



BB+C 

- 

- 

- 

35.8+-.1 

-1.5 (fix) 

°- 850 +.Im 

550/566 



C+C 

170 + 11 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

10 010+ 20 - 800 
iu.uiu_ 6 420 

546/566 



C + BB + PL5 p arams 

S2 til 8 

-0.7 (fix) 

- 

5 6.2+ + 5 

-1.5 (fix) 

- 

548/565 



C + BB + C5 params 


-0.7 (fix) 

- 

51.4+8 5 ° 

-1.5 (fix) 

9 099+ 21 500 
y.uyy_s Q20 

544/564 

+ 10.2 

+ 10.6 

Band 

185 + 7 

-0.95+0.04 

<-5 

- 

- 

- 

526/566 



C + BB 

203± 22 

-1.14+0.10 

- 

35-3l|;l 

- 

- 

522/565 



C+PL 

176+7 

-0.93+0.04 

- 

- 

-1.24 + 0.18 

- 

517/565 



C+BB+PL 


-1.11+0.11 

- 

3-5+;) 

! 9q +0.42 
-i.zo — o 36 

- 

513/563 



BB + C 

- 

- 

- 

34.0+1.8 

-1.5 (fix) 

0 502+ 0-093 
U.5U^ — 0.079 

536/566 



C+PL 

146+5 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

533/567 



C+C2 

174+17 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

0 297+ 9-978 
u ^ y -0.057 

531/566 



C + BB + PL5 p arams 

139+ 28 
— 56 

-0.7 (fix) 

- 

43.0+10.0 

-1.75+0.08 

- 

520/564 



C + BB + C5 params 

69 + 18 

-0.7 (fix) 

- 

48.1 + 5.5 

-1.5 (fix) 

39 70o+ 1 2 1000 
jy.' uu _2 3 .600 

523/564 

+ 10.6 

+ 10.9 

Band 

233+8 

-0.62+0.04 

<-5 

- 

- 

- 

513/566 



C + BB 

247+ 3 ° 

-0.77+0.10 

- 

44.7+° 9 3 

- 

- 

503/565 



C+PL 

222 + 9 

-0.42+0.14 

- 

- 

-1 74+ 0 ’ 17 
±- ' 4 —0.20 

- 

509/565 



C+BB+PL 

232 + 9 

-0.55+0.05 

- 

42.0+8-2 

, 71 +0.09 
-0.16 

- 

498/563 



C+BB 

199+ 19 

±yy_ 30 

-0.7 (fix) 

- 

52.4+7.8 

- 

- 

505/566 



C + BB + PL 5params 


-0.7 (fix) 

- 

52.3+6.5 

-1.5 (fix) 

- 

501/565 

+ 10.9 

+ 11.3 

Band 

303 + 15 

-0.71+0.04 

-3 14+0-33 
-0.84 

- 

- 

- 

511/566 



C + BB 

o««+29 

3®®_24 

-0.84+0.06 

- 

35.6+5.0 

- 

- 

505/565 



C+PL 

302 + 13 

-0.65+0.09 

- 

- 

1 58+0-54 

— 0.13 

- 

508/565 



C+BB+PL 

349 + 24 

-0.51+0.20 

- 

29.5+;3 

-1 79+0-09 

'-0.13 

- 

498/563 



C + BB 

293 + 10 

-0.7 (fix) 

- 

319+ 90 
dl-y -4.8 

- 

- 

515/566 



C+PL 

289+8 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

515/567 



C + BB + PL5 p arams 

294+12 

-0.7 (fix) 

- 

34 9+ 6 - 0 

^-4.3 

-1.5 (fix) 

- 

507/565 

+ 11.3 

+ 11.8 

Band 

221 + 13 

-0.76+0.04 

-2.39+0.11 

- 

- 

- 

666/566 



C + BB 

3 5 0±f 4 

-1.00+0.05 

- 

26.7+2.4 

- 

- 

662/565 



C+PL 

250 + 11 

-0.80+0.06 

- 

- 

44+0.19 

i44 -0.10 

- 

670/565 



C+BB+PL 


-0 93+ 9- 21 
uy -0.09 

- 

26.2 + 2.8 

-1.48 + 0.07 

- 

650/563 



C+PL 

218+6 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

657/567 



C + BB + PL5 p arams 

232 + 12 

-0.7 (fix) 


23.0+;! 

-1.5 (fix) 


650/565 
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1 Time fr 

om Tq 

Models 

Standard Model 

Additional Mo> 

del 

Cstat / dof 



C or Band 

BB 

PL 

or C 


^start 

Tstop 

Parameters 

^peak 


/3 

kT 

a 

E 0 


(s) 

(s) 


(keV) 



(keV) 


(MeV) 


+ 11.8 

+ 12.2 

Band 

230 + 12 

-0.89 + 0.04 

3 15+0-35 
- 5,ao -0.83 

_ 

_ 

_ 

516/566 



C + BB 

300+3| 

-1.09 + 0.07 

- 

31.0+3.4 

- 

- 

509/565 



C + PL 

227+14 

-0.75+0.14 

- 

- 

-1 79+ 0 - 20 
— 0.30 

- 

516/565 



C+BB+PL 

29l"^ 88 
^ yA —28 

70+0.27 

-0.78_o 24 

- 

26.s±|;X 

-2 00"*" 0 ' 11 
2UU -0.15 

- 

491/563 



C+PL 

183+6 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

516/567 



C + BB+PL 5params 

315 + 100 

-0.7 (fix) 

- 

31.4+2-6 

-1.5 (fix) 

- 

495/565 

+ 12.2 

+ 12.6 

Band 

148 + 17 

-0.96 + 0.08 

_ 2 39+0-15 
dy -0.26 

_ 

_ 

_ 

579/566 



C + BB 

212"*" 45 

21 —31 

-1 16+ 019 
aad -0.13 


20 - 5 ±t.5 

- 

- 

583/565 



C + PL 

158 + 16 

-0 92"*" 0-21 
u,y2 - 0.14 

- 

- 

-1 62+ 0 ' 25 
AO2 -0.11 

- 

581/565 



C+BB+PL 

164+11 

-0.79 + 0.08 

- 

9 9+3-5 
yy — 2.4 

-1 66+ 025 

1 DD —0.13 

- 

580/563 



C+PL 

124+6 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

587/567 



C+BB+PL5 p arams 

142 + 11 

-0.7 (fix) 

- 

5 9+ 21 
°' y —1.7 

-1.5 (fix) 

- 

581/565 

+ 12.6 

+ 12.9 

Band 

92 + 11 

-0.88 + 0.15 

.2 41+ 0-16 
Z - 4A —0.22 

- 

- 

- 

569/566 



C + BB 

33<A“ 8 

-1.63 + 0.03 

- 

49.5+2Q 

- 

- 

566/565 



C + PL 

S2±8 

-0 03^~ 9- ^ 4 
u u -0.58 

- 

- 

1 70+O.I3 
' A ' -0.08 

- 

566/565 



C + BB 

99±12 

-0.7 (fix) 

- 

4 8^” 2 ' 8 

4 -1.1 

- 

- 

560/563 



C + PL 

85±5 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

563/567 



C+BB+PL5 p arams 

70±7 

-0.7 (fix) 

- 

25.2++J 
— 3.6 

-1.66+0.06 

- 

555/565 

+ 12.9 

+ 13.1 

Band 


-1.05+0.08 

<-2.26 

- 

- 

- 

569/566 



C + BB 

691+^2 

-1.28+0.14 

- 

47.7+12.5 

- 

- 

566/565 



C + PL 

352+H 

-0 64+ 0-38 

U D -0.43 

- 

- 

-1.85+0.12 

- 

566/565 



C+BB+PL 


-0.63+0.12 

- 

27 9+ 9-0 

2 ‘ y —5.0 

2 4Q+0.36 
-2.49 — 0 44 

- 

560/563 



C + BB 

310”^ 33 
* 1U -29 

-0.7 (fix) 

- 

2.8+0.5 

- 

- 

578/566 



C+PL 

272+30 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

588/567 



C + BB+PL 5params 

797+629 

12 -349 

-0.7 (fix) 

- 

35.2+ 8 ;! 

-1.5 (fix) 

- 

575/565 

+ 13.1 

+ 13.3 

Band 

301+21 

-0.70+0.06 

<-3.38 




593/566 



C + BB 

3»7±i“ 

i 17+0-51 
11 -0.22 

- 

56.5± 8 ;J 

- 

- 

581/565 



C + PL 

268 + 22 

-0 16+ 0-29 
uad -0.41 

- 

- 

-1.93+0.12 

- 

583/565 



C+BB + PL5 p arams 

413 + 228 

4A,i —103 

-0.7 (fix) 

- 

45.5 + 8 

-1.5 (fix) 

- 

581/564 

+ 13.3 

+ 13.9 

Band 

130+9 

-0.62 + 0.09 

-2 60"*"°' 18 
^.bU_o 26 

- 

- 

- 

562/566 



C + BB 

971+108 

2 ' A — 62 

-1.27+0.14 

- 

25 2“*~-*" 8 
25.2-2 o 

- 

- 

558/565 



C+BB+PL5 p arams 

144+10 

-0.7 (fix) 

- 

23 2”*" 8 ' 8 
-5.9 

-1.5 (fix) 

- 

558/565 

+ 13.9 

+ 15.7 

Band 

136+ 28 

-1.16 + 0.08 

24+0.14 
^ -0.44 

- 

- 

- 

656/566 



C + BB 

228 + |^ 

-1.41 + 0.18 

- 

91 9 + 5.6 

21 -7.9 

- 

- 

656/565 



C + PL 

124 + 11 

-0.77+0.11 

- 

- 

-1.73+0.04 

- 

650/565 



C+BB+PL 

121 + 16 

-0.77+0.11 

- 

12.4+7-0 

-1.71 + 0.05 

- 

650/563 



C + BB 

136+8 

-0.7 (fix) 

- 

4.4 + 0.7 

- 

- 

672/566 



C+PL 

96+5 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

674/567 



BB+C 

- 

- 

- 

22-3ti;S 

-1.5 (fix) 

587+ 155 

5» 7 _116 

660/566 

+ 15.7 

+ 16.4 

Band 

122 + 8 

-0.65+0.08 

9 qq+ 0-27 
-z.oO — o 54 

- 

- 

- 

599/566 



C + BB 

136 + 12 

-0.92 + 0.09 

- 

28 +tit 

- 

- 

601/565 



C + PL 

128+5 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

609/567 



BB + C 

- 

- 

- 

26.7+1.6 

-1.5 (fix) 

,qq + 192 
-137 

615/566 



C+PL 

128+5 

-0.7 (fix) 

- 

- 

-1.5 (fix) 

- 

609/567 



C + BB+PLgpaj-ajng 

109 + 10 

-0.7 (fix) 

- 


-1.5 (fix) 

- 

600/565 

+ 16.4 

+20.0 

Band 

102"^” 38 
AU2 —25 

-1.52+0.10 

-2 07+ 0 ’ 11 
2U -0.17 

- 

- 

- 

679/566 



BB + PL 

- 

- 

- 

«.i+l:S 

-1.86+0.05 

- 

678/566 



C + BB + C5 params 

15 + 4 

-0.7 (fix) 

- 

18.8+3.0 

-1.5 (fix) 

1 fiftn +5.470 
ll380 _0.877 

684/564 




